We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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.

Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks.

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 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.

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.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present our analysis of the magnetic field structures from 6000 au to 100 au scales in the Class 0 protostar B335 inferred from our JCMT POL-2 observations and the ALMA archival polarimetric data. To interpret the observational results, we perform a series of (non-)ideal MHD simulations of the collapse of a rotating non-turbulent dense core, whose initial conditions are adopted to be the same as observed in B335, and generate synthetic polarization maps. The comparison of our JCMT and simulation results suggests that the magnetic field on a 6000 au scale in B335 is pinched and well aligned with the bipolar outflow along the east-west direction. Among all our simulations, the ALMA polarimetric results are best explained with weak magnetic field models having an initial mass-to-flux ratio of 9.6. However, we find that with the weak magnetic field, the rotational velocity on a 100 au scale and the disk size in our simulations are larger than the observational estimates by a factor of several. An independent comparison of our simulations and the gas kinematics in B335 observed with the SMA and ALMA favors strong magnetic field models with an initial mass-to-flux ratio smaller than 4.8. We discuss two possibilities resulting in the different magnetic field strengths inferred from the polarimetric and molecular-line observations, (1) overestimated rotational-to-gravitational energy in B335 and (2) additional contributions in the polarized intensity due to scattering on a 100 au scale.

Black-hole binary (BHB) systems comprise a stellar-mass black hole and a closely orbiting companion star. Matter is transferred from the companion to the black hole, forming an accretion disk, corona and jet structures. The resulting release of gravitational energy leads to emission of X-rays. The radiation is affected by special/general relativistic effects, and can serve as a probe of the properties of the black hole and surrounding environment, if the accretion geometry is properly identified. Two competing models describe the disk-corona geometry for the hard spectral state of BHBs, based on spectral and timing measurements. Measuring the polarization of hard X-rays reflected from the disk allows the geometry to be determined. The extent of the corona differs between the two models, affecting the strength of relativistic effects (e.g., enhancement of polarization fraction and rotation of polarization angle). Here, we report observational results on linear polarization of hard X-ray (19-181 keV) emission from a BHB, Cygnus X-1, in the hard state. The low polarization fraction, <8.6% (upper limit at 90% confidence level), and the alignment of the polarization angle with the jet axis show that the dominant emission is not influenced by strong gravity. When considered together with existing spectral and timing data, our result reveals that the accretion corona is either an extended structure, or is located far from the black hole in the hard state of Cygnus X-1.

We report $0.14"$ resolution observations of the dust continuum at band 7, and the CO(3--2) and HCO$^{+}$(4--3) line emissions toward the transitional disk around Sz 91 with Atacama Large Millimeter/submillimeter Array (ALMA). The dust disk appears to be an axisymmetric ring, peaking a radius of $\sim$95~au from a Gaussian fit. The Gaussian fit widths of the dust ring are 24.6 and 23.7~au for the major and the minor axes, respectively, indicating that the dust ring is not geometrically thin. The gas disk extends out to $\sim$320~au and is also detected in the inner hole of the dust ring. A twin-line pattern is found in the channel maps of CO, which can be interpreted as the emission from the front and rear of the flared gas disk. We perform the radiative transfer calculations using RADMC-3D, to check whether the twin-line pattern can be reproduced under the assumption that the flared gas disk has a power-law form for the column density and $T_\mathrm{gas}=T_\mathrm{dust}$. The thermal Monte Carlo calculation in RADMC-3D shows that the disk temperature has a gradient along the vertical direction beyond the dust ring, as it blocks the stellar radiation, and thus the twin-line pattern can be naturally explained by the flared gas disk in combination with the dust ring. In addition, no significant depletion of the CO molecules in the cold midplane achieves a reasonable agreement with the observed twin-line pattern. This result indicates that the CO emission from the rear surface must be heavily absorbed in the cold midplane.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal MHD effects proposed to enable the formation of Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare gas kinematics traced by the H13CO+ (3-2) and C18O (2-1) lines in the Class 0 protostar B335 obtained with our ALMA observations. A central compact (~1"-2") component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. We constructed kinematical models to fit the observed velocity structures and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. The infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85+/-0.2 km/s at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km/s at a radius of 100 au. The Hall parameter for H13CO+ is estimated to be >>1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between ions and neutrals could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0.1") are needed to examine this possibility and ambipolar diffusion on a smaller scale.

The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335.

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.

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.

Strong magnetic fields, synchrotron emission, and Compton scattering are omnipresent in compact celestial X-ray sources. Emissions in the X-ray energy band are consequently expected to be linearly polarized. X-ray polarimetry provides a unique diagnostic to study the location and fundamental mechanisms behind emission processes. The polarization of emissions from a bright celestial X-ray source, the Crab, is reported here for the first time in the hard X-ray band (~20-160 keV). The Crab is a complex system consisting of a central pulsar, a diffuse pulsar wind nebula, as well as structures in the inner nebula including a jet and torus. Measurements are made by a purpose-built and calibrated polarimeter, PoGO+. The polarization vector is found to be aligned with the spin axis of the pulsar for a polarization fraction, PF = (20.9 $\pm$ 5.0)%. This is higher than that of the optical diffuse nebula, implying a more compact emission site, though not as compact as, e.g., the synchrotron knot. Contrary to measurements at higher energies, no significant temporal evolution of phase-integrated polarisation parameters is observed. The polarization parameters for the pulsar itself are measured for the first time in the X-ray energy band and are consistent with observations at optical wavelengths.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.

Polarimetric observations of celestial sources in the hard X-ray band stand to provide new information on emission mechanisms and source geometries. PoGO+ is a Compton scattering polarimeter (20-150 keV) optimised for the observation of the Crab (pulsar and wind nebula) and Cygnus X-1 (black hole binary), from a stratospheric balloon-borne platform launched from the Esrange Space Centre in summer 2016. Prior to flight, the response of the polarimeter has been studied with polarised and unpolarised X-rays allowing a Geant4-based simulation model to be validated. The expected modulation factor for Crab observations is found to be $M_{\mathrm{Crab}}=(41.75\pm0.85)\%$, resulting in an expected Minimum Detectable Polarisation (MDP) of $7.3\%$ for a 7 day flight. This will allow a measurement of the Crab polarisation parameters with at least $5\sigma$ statistical significance assuming a polarisation fraction $\sim20\%$ $-$ a significant improvement over the PoGOLite Pathfinder mission which flew in 2013 and from which the PoGO+ design is developed.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and $^{12}$CO, $^{13}$CO, and C$^{18}$O J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ~0".25 (35 AU). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-$\mu$m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius $\sim$70 AU that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 AU, or an edge that transitions smoothly over an annular width of 10 AU near 60 AU. In gas, the data are consistent with a cavity that is smaller, about 15 AU in radius, and whose surface density at 15 AU is $10^{3\pm1}$ times smaller than the surface density at 70 AU; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 AU or a double-drop model as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities, e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation.

We present ALMA 1.3 mm continuum, 12CO, C18O, and SO data for the Class 0 protostars, Lupus 3 MMS, IRAS 15398-3559, and IRAS 16253-2429 at resolutions of ~100 AU. By measuring a rotational profile in C18O, a 100 AU Keplerian disk around a 0.3 Msun protostar is observed in Lupus 3 MMS. No 100 AU Keplerian disks are observed in IRAS 15398-3559 and IRAS 16253-2429. Nevertheless, embedded compact (<30 AU) continuum components are detected. The C18O emission in IRAS 15398-3559 shows signatures of infall with a constant angular momentum. IRAS 16253-2429 exhibits signatures of infall and rotation, but its rotational profile is unresolved. By fitting the C18O data with our kinematic models, the protostellar masses and the disk radii are inferred to be 0.01 Msun and 20 AU in IRAS 15398-3559, and 0.03 Msun and 6 AU in IRAS 16253-2429. By comparing the specific angular momentum profiles from 10,000 to 100 AU in 8 Class 0 and I protostars, we find that the evolution of envelope rotation can be described with conventional inside-out collapse models. In comparison with a sample of 18 protostars with known disk radii, our results reveal signs of disk growth, with the disk radius increasing as M*^{0.8+/-0.14} or t^{1.09+\-0.37} in the Class 0 stage, where M* is the protostellar mass and t is the age. The disk growth rate slows down in the Class I stage. Besides, we find a hint that the mass accretion rate declines as t^{-0.26+\-0.04} from the Class 0 to I stages.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

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.

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.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.

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.

The formation scenario of a gapped disk, i.e., transitional disk, and its asymmetry is still under debate. Proposed scenarios such as disk-planet interaction, photoevaporation, grain growth, anticyclonic vortex, eccentricity, and their combinations would result in different radial distributions of the gas and the small (sub-$\mu$m size) and large (millimeter size) dust grains as well as asymmetric structures in a disk. Optical/near-infrared (NIR) imaging observations and (sub-)millimeter interferometry can trace small and large dust grains, respectively; therefore multi-wavelength observations could help elucidate the origin of complicated structures of a disk. Here we report SMA observations of the dust continuum at 1.3~mm and $^{12}$CO~$J=2\rightarrow1$ line emission of the pre-transitional protoplanetary disk around the solar-mass star PDS~70. PDS~70, a weak-lined T Tauri star, exhibits a gap in the scattered light from its disk with a radius of $\sim$65~AU at NIR wavelengths. However, we found a larger gap in the disk with a radius of $\sim$80~AU at 1.3~mm. Emission from all three disk components (the gas and the small and large dust grains) in images exhibits a deficit in brightness in the central region of the disk, in particular, the dust-disk in small and large dust grains has asymmetric brightness. The contrast ratio of the flux density in the dust continuum between the peak position to the opposite side of the disk reaches 1.4. We suggest the asymmetries and different gap-radii of the disk around PDS~70 are potentially formed by several (unseen) accreting planets inducing dust filtration.

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.

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 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.

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

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

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 interferometric observations of the VeLLO L1521F-IRS, in 12CO (2-1) line emission and the 1.3 mm continuum emission, using the Submillimeter Array (SMA). With the 12CO high-resolution observations, we have spatially resolved a compact but poorly collimated molecular outflow associated with L1521F-IRS for the first time. The blueshifted and redshifted lobes are aligned along the east and west side of L1521F-IRS with a lobe size of ~1000 AU. The estimated outflow mass, maximum outflow velocity, and outflow force are (9.0-80)x10^{-4} M_{sol}, 7.2 km s^{-1}, and (7.4-66)x10^{-7} M_{sol} km s^{-1} yr^{-1}, respectively. The estimated outflow parameters such as size, mass, and momentum rate are similar to values derived for other VeLLOs, and are located at the lower end of values compared to previously studied outflows associated with low- to high-mass star forming regions. Low-velocity less collimated (1.5 km s^{-1}/1200 AU) and higher-velocity compact (4.0 km s^{-1}/920 AU) outflow components are suggested by the data. These velocity structures are not consistent with those expected in the jet driven or wind driven outflow models, perhaps suggesting a remnant outflow from the FHSC as well as an undeveloped outflow from the protostar. Detection of an infrared source and compact millimeter continuum emission suggest the presence of the protostar, while its low bolometric luminosity (0.034-0.07 L$_{\odot}$), and small outflow, suggests that L1521F is in the earliest protostellar stage ($<10^{4}$ yr) and contains a substellar mass object. The bolometric (or internal) luminosity of L1521F-IRS suggests that the current mass accretion rate is an order-of-magnitude lower than expected in the standard mass accretion model (~10^{-6} M_{sol} yr^{-1}), which may imply that L1521F-IRS is currently in a low activity phase.

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.

The circumstellar disk of AB Aurigae has garnered strong attention owing to the apparent existence of spirals at a relatively young stage and also the asymmetric disk traced in thermal dust emission. However, the physical conditions of the spirals are still not well understood. The origin of the asymmetric thermal emission is unclear. We observed the disk at 230 GHz (1.3 mm) in both the continuum and the spectral line ^12CO J=2-1 with IRAM 30-m, the Plateau de Bure interferometer, and the Submillimeter Array to sample all spatial scales from 0.37" to about 50". To combine the data obtained from these telescopes, several methods and calibration issues were checked and discussed. The 1.3 mm continuum (dust) emission is resolved into inner disk and outer ring. Molecular gas at high velocities traced by the CO line is detected next to the stellar location. The inclination angle of the disk is found to decrease toward the center. On a larger scale, based on the intensity weighted dispersion and the integrated intensity map of ^12CO J=2-1, four spirals are identified, where two of them are also detected in the near infrared. The total gas mass of the 4 spirals (M_spiral) is 10^-7 < M_spiral < 10^-5 M_sun, which is 3 orders of magnitude smaller than the mass of the gas ring. Surprisingly, the CO gas inside the spiral is apparently counter-rotating with respect to the CO disk, and it only exhibits small radial motion. The wide gap, the warped disk, and the asymmetric dust ring suggest that there is an undetected companion with a mass of 0.03 M_sun at a radius of 45 AU. Although an hypothetical fly-by cannot be ruled out, the most likely explanation of the AB Aurigae system may be inhomogeneous accretion well above or below the main disk plane from the remnant envelope, which can explain both the rotation and large-scale motions detected with the 30-m image.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

We present 1.3-mm subarcsecond SMA observations of the prototypical Class 0 protostar VLA1623. We report the detection of 1.3-mm continuum emission both from the central protostellar component VLA1623 and two additional sources, Knot-A and Knot-B, which have been already detected at longer wavelengths. Knot-A and Knot-B are both located along the western cavity wall opened by the protostellar outflow from VLA1623. Our SMA observations moreover show that these two continuum sources are associated with bright, high-velocity 12CO(2-1) emission, slightly shifted downstream of the outflow propagation direction with respect to the 1.3-mm continuum emission peaks. The alignment of Knot-A and Knot-B along the protostellar outflow cavity, the compactness of their 1.3-mm continuum emission and the properties of the associated CO emission suggest that these two sources trace outflow features due to shocks along the cavity wall, rather than protostellar objects. While it was considered as one of the best examples of a close protobinary system so far, the present analysis suggests that the prototypical Class 0, VLA1623, is single on the scales a>100 AU probed by our SMA observations. Moreover, we present here the second robust case of compact millimeter continuum emission produced by interactions between the protostellar jet and the envelope of a Class 0 protostar, which suggests a high occurrence of these outflow features during the embedded phase.

FS Tau B is one of the few T Tauri stars that possess a jet and a counterjet as well as an optically-visible cavity wall. We obtained images and spectra of its jet-cavity system in the near-infrared H and K bands using Subaru/IRCS and detected the jet and the counterjet in the [Fe II] 1.644 \mu m line for the first time. Within the inner 2" the blueshifted jet is brighter, whereas beyond ~ 5" the redshifted counterjet dominates the [Fe II] emission. The innermost blueshifted knot is spectrally resolved to have a large line width of ~ 110 km/s, while the innermost redshifted knot appears spectrally unresolved. The velocity ratio of the jet to the counterjet is ~ 1.34, which suggests that FS Tau B is driving an asymmetric jet, similar to those found in several T Tauri Stars. Combining with optical observations in the literature, we showed that the blueshifted jet has lower density and higher excitation than the redshifted counterjet. We suggest that the asymmetry in brightness and velocity is the manifestation of a bipolar outflow driving at different mass-loss rates, while maintaining balance of linear momentum. A full explanation to the asymmetry in the FS Tau B system awaits detail modeling and further investigation of the kinematic structure of the wind-associated cavity walls.

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.

We present \sim1" resolution 1.3 millimeter dust continuum and spectral line (12CO and 13CO J=2-1) observations of the transitional disk system HD 135344B obtained with the Submillimeter Array. The disk shows a Keplerian rotation pattern with an inclination of \sim11{\deg}, based on the spatially and spectrally resolved 12CO and 13CO emission. The data show clear evidence for both dust and gas surface density reductions in the inner region of the disk (radius \lesssim 50 AU) from the continuum and 13CO J=2-1 data, respectively. The presence of this inner cavity in both the dust and gas is more consistent with clearing by giant planet formation than by photoevaporation or by grain growth. There is also an indication of global CO gas depletion in the disk, as the mass estimated from 13CO emission (\sim3.8 \times 10^(-4) M\odot) is about two orders of magnitude lower than that derived from the 1.3 mm continuum (\sim2.8 \times 10^(-2) M\odot).

We report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

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

FN Tau is a rare example of very low-mass T Tauri stars that exhibits a spatially resolved nebulosity in near-infrared scattering light. To directly derive the parameters of a circumstellar disk around FN Tau, observations of dust continuum emission at 340 GHz are carried out with the Submillimeter Array (SMA). A point-like dust continuum emission was detected with a synthesized beam of ~ 0.7" in FWHM. From the analysis of the visibility plot, the radius of the emission is estimated to be <= 0.29", corresponding to 41 AU. This is much smaller than the radius of the nebulosity, 1.85" for its brighter part at 1.6 micron. The 340 GHz continuum emission observed with the SMA and the photometric data at lambda <= 70 micron are explained by a power-law disk model whose outer radius and mass are 41 AU and (0.24 - 5.9) x 10^{-3} M_{sun}, respectively, if the exponent of dust mass opacity (beta) is assumed to be 0-2. The disk model cannot fully reproduce the flux density at 230 GHz obtained with the IRAM 30-meter telescope, suggesting that there is another extended "halo" component that is missed in the SMA observations. By requiring the halo not to be detected with the SMA, the lower limit to the size of the halo is evaluated to be between 174 AU and 574 AU, depending on the assumed beta value. The physical nature of the halo is unclear, but it may be the periphery of a flared circumstellar disk that is not described well in terms of a power-law disk model, or a remnant of a protostellar envelope having flattened structure.

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.

We present Submillimeter Array observations of the GQ Lup system at 1.3 millimeters wavelength with $0\farcs4$ ($\sim$60 AU) resolution. Emission is detected from the position of the primary star, GQ Lup A, and is marginally resolved. No emission is detected from the substellar companion, GQ Lup B, $0\farcs7$ away. These data, together with models of the spectral energy distribution, suggest a compact disk around GQ Lup A with mass $\sim 3$ M$_{Jup}$, perhaps truncated by tidal forces. There is no evidence for a gap or hole in the disk that might be the signature of an additional inner companion body capable of scattering GQ Lup B out to $\sim100$ AU separation from GQ Lup A. For GQ Lup B to have formed {\it in situ}, the disk would have to have been much more massive and extended.

We present potential detections of H-band scattered light emission around four young star, selected from a total sample of 45 young stars observed with the CIAO coronagraph of the Subaru telescope. Two CTTS, CI Tau and DI Cep, and two WTTS, LkCa 14 and RXJ 0338.3+1020 were detected. In all four cases, the extended emission is within the area of the residual PSF halo, and is revealed only through careful data reduction. We compare the observed extended emission with simulations of the scattered light emission, to evaluate the plausibility and nature of the detected emission.

We made mid-infrared observations of the 10Msun Herbig Be star HD200775 with the Cooled Mid-Infrared Camera and Spectrometer (COMICS) on the 8.2m Subaru Telescope. We discovered diffuse emission of an elliptical shape extended in the north-south direction inabout 1000AU radius around unresolved excess emission. The diffuse emission is perpendicular to the cavity wall formed by the past outflow activity and is parallel to the projected major axis of the central close binary orbit. The centers of the ellipse contours of the diffuse emission are shifted from the stellar position and the amount of the shift increases as the contour brightness level decreases. The diffuse emission is well explained in all of geometry, size, and configuration by an inclined flared disk where only its surface emits the mid-infrared photons. Our results give the first well-resolved infrared disk images around a massive star and strongly support that HD200775 is formed through the disk accretion. The disk survives the main accretion phase and shows a structure similar to that around lower-mass stars with 'disk atmosphere'. At the same time, the disk also shows properties characteristic to massive stars such as photoevaporation traced by the 3.4mm free-free emission and unusual silicate emission with a peak at 9.2micron, which is shorter than that of many astronomical objects. It provides a good place to compare the disk properties between massive and lower-mass stars.

High mass star formation and the evolution of HII regions have a substantial impact on the morphology and star formation history of molecular clouds. The HII region Gum 48d, located in the Centaurus Arm at a distance of 3.5 kpc, is an old, well evolved HII region whose ionizing stars have moved off the main sequence. As such, it represents a phase in the evolution of HII regions that is less well studied than the earlier, more energetic, main sequence phase. In this paper we use multi-wavelength archive data from a variety of sources to perform a detailed study of this interesting region. Morphologically, Gum 48d displays a ring-like faint HII region associated with diffuse emission from the associated PDR, and is formed from part of a large, massive molecular cloud complex. There is extensive ongoing star formation in the region, at scales ranging from low to high mass, which is consistent with triggered star formation scenarios. We investigate the dynamical history and evolution of this region, and conclude that the original HII region was once larger and more energetic than the faint region currently seen. The proposed history of this molecular cloud complex is one of multiple, linked generations of star formation, over a period of 10 Myr. Gum 48d differs significantly in morphology and star formation that the other HII regions in the molecular cloud; these differences are likely the result of the advanced age of the region, and its different evolutionary status.

Studying the earliest stages in the birth of stars is crucial for understanding how they form. Brown dwarfs with masses between that of stars and planets are not massive enough to maintain stable hydrogen-burning fusion reactions during most of their lifetime. Their origins are subject to much debate in recent literature because their masses are far below the typical mass where core collapse is expected to occur. We present the first confirmed evidence that brown dwarfs undergo a phase of molecular outflow that is typical of young stars. Using the Submillimeter Array, we have obtained a map of a bipolar molecular outflow from a young brown dwarf. We estimate an outflow mass of 1.6 x 10^-4 M_Sun and a mass-loss rate of 1.4 x 10^-9 M_Sun. These values are over two orders of magnitude smaller than the typical ones for T Tauri stars. From our millimiter continuum data and our own analysis of Spitzer infrared photometry, we estimate that the brown dwarf has a disk with a mass of 8 x 10^-3 M_Sun and an outer disk radius of 80 AU. Our results demonstrate that the bipolar molecular outflow operates down to planetary masses, occurring in brown dwarfs as a scaled-down version of the universal process seen in young stars.

Using the ASTE 10 m submillimeter telescope and the 1.4 m Infrared Survey Facility (IRSF), we performed an extensive outflow survey in the Orion Molecular Cloud -2 and -3 region. Our survey, which includes 41 potential star-forming sites, has been newly compiled using multi-wavelength data based on millimeter- and submillimeter-continuum observations as well as radio continuum observations. From the CO (3-2) observations performed with the ASTE 10 m telescope, we detected 14 CO molecular outflows, seven of which were newly identified. This higher detection rate, as compared to previous CO (1-0) results in the same region, suggests that CO (3-2) may be a better outflow tracer. Physical properties of these outflows and their possible driving sources were derived. Derived parameters were compared with those of CO outflows in low- and high-mass starforming regions. We show that the CO outflow momentum correlates with the bolometric luminosity of the driving source and with the envelope mass, regardless of the mass of the driving sources. In addition to these CO outflows, seven sources having NIR features suggestive of outflows were also identified.

We have mapped the proto-binary source IRAS 16293-2422 in CO 2-1, 13CO 2-1, and CO 3-2 with the Submillimeter Array (SMA). The maps with resolution of 1".5-5" reveal a single small scale (~3000 AU) bipolar molecular outflow along the east-west direction. We found that the blueshifted emission of this small scale outflow mainly extends to the east and the redshifted emission to the west from the position of IRAS 16293A. A comparison with the morphology of the large scale outflows previously observed by single-dish telescopes at millimeter wavelengths suggests that the small scale outflow may be the inner part of the large scale (~15000 AU) E-W outflow. On the other hand, there is no clear counterpart of the large scale NE-SW outflow in our SMA maps. Comparing analytical models to the data suggests that the morphology and kinematics of the small scale outflow can be explained by a wide-angle wind with an inclination angle of ~30-40 degrees with respect to the plane of the sky. The high resolution CO maps show that there are two compact, bright spots in the blueshifted velocity range. An LVG analysis shows that the one located 1" to the east of source A is extremely dense, n(H_2)~10^7 cm^-3, and warm, T_kin >55 K. The other one located 1" southeast of source B has a higher temperature of T_kin >65 K but slightly lower density of n(H_2)~10^6 cm^-3. It is likely that these bright spots are associated with the hot core-like emission observed toward IRAS 16293. Since both two bright spots are blueshifted from the systemic velocity and are offset from the protostellar positions, they are likely formed by shocks.

We present the results of the interferometric observations of the circumstellar disk surrounding MWC 297 in the continuum at 230 GHz (1.3 mm) and in the (J=2-1) rotational transitions of $^{12}$CO,$^{13}$CO and C$^{18}$O using the Submillimeter Array. At a distance of 250 pc, MWC 297 is one of the closest, young massive stars (M$_{\star}$ $\sim$10 M$_{\odot}$) to us. Compact continuum emission is detected towards MWC 297 from which we estimate a disk mass (gas+dust) of 0.07 M$_{\odot}$ and a disk radius of $\le$ 80 AU. Our result demonstrates that circumstellar disks can survive around massive stars well into their main sequence phase even after they have become optically visible. Complementing our observations with the data compiled from the literature, we find the submm dust opacity index $\beta$ to be between 0.1 and 0.3. If the emission is optically thin, the low value of $\beta$ indicates the presence of relatively large grains in the disk, possibly because of grain growth. We do not detect any CO emission associated with the continuum source. We argue that the $^{13}$CO emission from the disk is likely optically thin, in which case, we derive an upper limit to the gas mass which implies significant depletion of molecular gas in the disk. The mass of this disk and the evolutionary trends observed are similar to those found for intermediate mass Herbig Ae stars and low mass T Tauri stars.

(Abridged) With the SMA we have made high angular-resolution (~1" = 160 AU) observations of the protobinary system IRAS 16293-2422 in the J = 4-3 lines of HCN and HC^15N, and in the continuum at 354.5 GHz. The HCN (4-3) line was also observed using the JCMT to supply missing short spacing information. Submillimeter continuum emission is detected from the individual binary components with a separation of ~5". The HC^15N (4-3) emission has revealed a compact (~500 AU) flattened structure (P.A. = -16 degree) associated with Source A, and shows a velocity gradient along the projected minor axis, which can be interpreted as an infalling gas motion. Our HCN image including the short-spacing information shows an extended (~3000 AU) circumbinary envelope as well as the compact structure associated with Source A. A toy model consisting of a flattened structure with radial infall towards a 1 Msun central star reproduces the HCN/HC^15N position-velocity diagram along the minor axis of the HC^15N emission. In the extended envelope there is also a North-East (Blue) to South-West (Red) velocity gradient across the binary alignment, which is likely to reflect gas motion in the swept-up dense gas associated with the molecular outflow from Source A. At Source B, there is only a weak compact structure with much narrower line widths (~2 km/s) seen in the optically-thin HC^15N emission than that at Source A (>10 km/s), and there is no clearly defined bipolar molecular outflow associated with Source B. These results imply the different evolutionary stages between Source A and B in the common circumbinary envelope.

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

We report a dual-band observation at 223 and 654 GHz (460 micron) toward an ultracompact (UC) HII region, G240.31+0.07, with the Submillimeter Array. With a beam size of 1"5 X 0"8, the dust continuum emission is resolved into two clumps, with clump A coincident well with an H2O maser and the UC HII region. The newly discovered clump, B, about 1"3 (~8.3 X 10^3 AU) to the southwest of clump A, is also associated with H2O masers and may be a more recent star-forming site. The continuum flux densities imply an opacity spectral index of \beta = 1.5 +- 0.3, suggestive of a value lower than the canonical 2.0 found in the interstellar medium and in cold, massive cores. The presence of hot (~100 K) molecular gas is derived by the brightness ratio of two H2CO lines in the 223 GHz band. A radial velocity difference of 2.5 +- 0.4 km/s is found between the two clumps in C18O (6-5) emission. The total (nebular and stellar) mass of roughly 58 Msun in the central region is close to, but not by far larger than, the minimum mass required for the two clumps to be gravitationally bound for binary rotation. Our continuum data do not suggest a large amount of matter associated with the H2 knots that were previously proposed to arise from a massive disk or envelope.

The circumstellar dust disk of the Herbig Ae star AB Aur has been found to exhibit complex spiral-like structures in the near-IR image obtained with the Subaru Telescope. We present maps of the disk in both 12CO (3-2) and dust continuum at 345 GHz with the SMA at an angular resolution of 1.0"x0.7" (144AU x 100AU). The continuum emission traces a dust disk with a central depression and a maximum overall dimension of 450AU (FWHM). This dust disk exhibits several distinct peaks that appear to coincide with bright features in the near-IR image, in particular the brightest inner spiral arm. The CO emission traces a rotating gas disk of size 530AU x 330AU with a deprojected maximum velocity of 2.8km/s at 450AU. In contrast to the dust disk, the gas disk exhibits an intensity peak at the stellar position. Furthermore, the CO emission in several velocity channels traces the innermost spiral arm seen in the near-IR. We compare the observed spatial-kinematic structure of the CO emission to a simple model of a disk in Keplerian rotation, and find that only the emission tracing the main spiral arm clearly lies outside the confines of our model. This emission has a net outward radial motion compared with the radial velocity predicted by the model at the location of the main spiral arms. The disk of AB Aur is therefore quite different from the Keplerian disks seen around many Herbig Ae stars. The spiral-like structures of the disk with non-Keplerian motions we revealed in 12CO (3-2), together with the central depression of the dust disk, may be explained to be driven by the possible existence of a giant planet forming in the disk.

Arcsecond-resolution spectral observations toward the protobinary system IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillimeter Array. Complex organic molecules such as CH3OH and HCOOCH3 were detected. Together with the rich organic inventory revealed, it clearly indicates the existence of two, rather than one, compact hot molecular cores (smaller than or equal to 400 AU in radius) associated with each of the protobinary components identified by their dust continuum emission in the inner star-forming core.

We present the first interferometric CO(J=3-2) observations (beam size of 3.9"x1.6" or 160pc x 65pc) with the Submillimeter Array (SMA) toward the center of the Seyfert 2 galaxy M51. The image shows a strong concentration at the nucleus and weak emission from the spiral arm to the northwest. The integrated intensity of the central component in CO(J=3-2) is almost twice as high as that in CO(J=1-0), indicating that the molecular gas within an ~80 pc radius of the nucleus is warm (>~100 K) and dense (~10^4 cm^-3). Similar intensity ratios are seen in shocked regions in our Galaxy, suggesting that these gas properties may be related to AGN or starburst activity. The central component shows a linear velocity gradient (~1.4 km/s/pc) perpendicular to the radio continuum jet, similar to that seen in previous observations and interpreted as a circumnuclear molecular disk/torus around the Seyfert 2 nucleus. In addition, we identify a linear velocity gradient (~0.7 km/s/pc) along the jet. Judging from the energetics, the velocity gradient can be explained by supernova explosions or energy and momentum transfer from the jet to the molecular gas via interaction, which is consistent with the high intensity ratio.

Submillimeter Array (SMA) observations of the high-mass star-forming region IRAS 18089-1732 in the 1 mm and 850 $\mu$m band with 1 GHz bandwidth reveal a wealth of information. We present the observations of 34 lines from 16 different molecular species. Most molecular line maps show significant contributions from the outflow, and only few molecules are confined to the inner core. We present and discuss the molecular line observations and outline the unique capabilities of the SMA for future imaging line surveys at high spatial resolution.

We present initial results of an ongoing search for interferometric calibrators at submillimeter (sub-mm) wavelengths with the Submillimeter Array (SMA). Powerful radio galaxies are commonly used as calibrators at centimeter and millimeter wavelengths, but many are not strong enough to serve as calibrators at sub-mm wavelengths because of their rapidly declining flux densities toward shorter wavelengths. The inability to find a calibrator close to the target source may limit or even prevent us from imaging many interesting sources at sub-mm wavelengths. Here, we investigate whether high-mass protostellar objects and ultracompact HII regions can serve as useful calibrators for the SMA. The dust emission associated with these objects makes them among the brightest sub-mm sources in the sky. Our observations at 0.85 mm (345 GHz) with an angular resolution of ~3" reveal that although a large fraction of the dust emission originates from an extended ``halo'' component, a compact unresolved component often remains that when sufficiently strong may serve as a useful calibrator. These observations also provide a first glimpse at the small-scale distribution of dust around ultracompact HII regions and high-mass protostellar objects at sub-mm wavelengths. We discuss the origin of the core-halo structure seen in many sources, and conclude with suggestions for future searches for calibrators with the SMA.

We have imaged the circumstellar envelope around the binary protostar L1551 IRS 5 in CS (J=7-6) and 343 GHz continuum emission at ~ 3 arcsec resolution using the Submillimeter Array. The continuum emission shows an elongated structure (~ 220 x 100 AU) around the binary perpendicular to the axis of the associated radio jet. The CS emission extends over ~ 400 AU, appears approximately circularly symmetric, and shows a velocity gradient from southeast (blueshifted) to northwest (redshifted). The direction of the velocity gradient is different from that observed in C18O (J=1-0). This may be because rotation is more dominant in the CS envelope than the C18O envelope, in which both infall and rotation exist. The CS emission may be divided into two velocity components: (1) a "high" velocity disk-like structure surrounding the protostar, +-1.0 - 1.5 km s-1 from the systemic velocity, and (2) a "low" velocity structure, located southwest of the protostar, < 1.0 km s -1 from the systemic velocity. The high-velocity component traces warm and dense gas with kinematics consistent with rotation around the protostar. The low-velocity component may arise from dense gas entrained in the outflow. Alternatively, this component may trace infalling and rotating gas in an envelope with a vertical structure.

We have mapped the ultracompact HII region, G5.89-0.39, and its molecular surroundings with the Submillimeter Array at 2".8 x 1".8 angular resolution in 1.3 mm continuum, SiO(5-4), and eight other molecular lines. We have resolved for the first time the highly energetic molecular outflow in this region. At this resolution, the outflow is definitely bipolar and appears to originate in a 1.3 mm continuum source. The continuum source peaks in the center of the HII region. The axis of the outflow lines up with a recently discovered O5V star.

We present ~2"-4" aperture synthesis observations of the circumstellar disk surrounding the nearby young star TW Hya in the CO J=2--1 and J=3--2 lines and associated dust continuum obtained with the partially completed Submillimeter Array. The extent and peak flux of the 230 and 345 GHz dust emission follow closely the predictions of the irradiated accretion disk model of Calvet et al. (2002). The resolved molecular line emission extends to a radius of at least 200 AU, the full extent of the disk visible in scattered light, and shows a clear pattern of Keplerian rotation. Comparison of the images with 2D Monte Carlo models constrains the disk inclination angle to 7+/-1 degrees. The CO emission is optically thick in both lines, and the kinetic temperature in the line formation region is ~20K. Substantial CO depletion, by an order of magnitude or more from canonical dark cloud values, is required to explain the characteristics of the line emission.

SMA observations of the massive star-forming region IRAS 18089-1732 in the 1mm and 850mu band reveal outflow and disk signatures in different molecular lines. The SiO(5--4) data show a collimated outflow in the northern direction. In contrast, the HCOOCH3(20--19) line, which traces high-density gas, is confined to the very center of the region and shows a velocity gradient across the core. The HCOOCH3 velocity gradient is not exactly perpendicular to the outflow axis but between an assumed disk plane and the outflow axis. We interpret these HCOOCH3 features as originating from a rotating disk that is influenced by the outflow and infall. Based on the (sub-)mm continuum emission, the mass of the central core is estimated to be around 38M_sun. The dynamical mass derived from the HCOOCH3 data is 22Msun, of about the same order as the core mass. Thus, the mass of the protostar/disk/envelope system is dominated by its disk and envelope. The two frequency continuum data of the core indicate a low dust opacity index beta ~ 1.2 in the outer part, decreasing to beta ~ 0.5 on shorter spatial scales.

The molecular evolution that occurs in collapsing prestellar cores is investigated. To model the dynamics, we adopt the Larson-Penston (L-P) solution and analogues with slower rates of collapse. For the chemistry, we utilize the new standard model (NSM) with the addition of deuterium fractionation and grain-surface reactions treated via the modified rate approach. The use of surface reactions distinguishes the present work from our previous model. We find that these reactions efficiently produce H2O, H2CO, CH3OH, N2, and NH3 ices. In addition, the surface chemistry influences the gas-phase abundances in a variety of ways. The current reaction network along with the L-P solution allows us to reproduce satisfactorily most of the molecular column densities and their radial distributions observed in L1544. The agreement tends to worsen with models that include strongly delayed collapse rates. Inferred radial distributions in terms of fractional abundances are somewhat harder to reproduce. In addition to our standard chemical model, we have also run a model with the UMIST gas-phase chemical network. The abundances of gas-phase S-bearing molecules such as CS and CCS are significantly affected by uncertainties in the gas-phase chemical network. In all of our models, the column density of N2H+ monotonically increases as the central density of the core increases during collapse from 3 10^4 cm-3 to 3 10^7 cm-3. Thus, the abundance of this ion can be a probe of evolutionary stage. Molecular D/H ratios in assorted cores are best reproduced in the L-P picture with the conventional rate coefficients for fractionation reactions. If we adopt the newly measured and calculated rate coefficients, the D/H ratios, especially N2D+/N2H+, become significantly lower than the observed values.

We have carried out mapping observations of the molecular core associated with the young Class 0 protostar, IRAM 04191+1522, in the CH3OH (JK=2K-1K) and C34S (J=2-1) lines using the 45 m telescope at Nobeyama Radio Observatory. As well as an elongated envelope associated with the protostellar formation (size \~0.07 pc x 0.04 pc, mass ~ 2.3 Mo), two compact (~ 0.03 pc) condensations were found in the CH3OH line at the southern edge of the elongated envelope, where the blueshifted CO outflow emerging from the protostar is located. In contrast to the elongated envelope, those compact CH3OH condensations show much larger line width (up to 2.0 km s-1) with centroid velocities blueshifted by ~ 0.8 km s-1. The compact condensations have momenta (~ 0.06 Mo km s-1) comparable to that of the blueshifted molecular outflow. These results suggest that the compact condensations are probably formed in the course of interaction between the outflow and the ambient gas surrounding the protostar, and that such interaction may cause dissipation of a part of the ambient gas. No drastic, localized enhancement of the CH3OH abundance is, however, observed toward the compact condensations, implying that there seems to be no significant shock heating at the compact condensations. This may be because the CO outflow velocity (< 10 km s-1) is too low to cause effective heating to release CH3OH on dust grains into gas phase.

We have investigated the evolution and distribution of molecules in collapsing prestellar cores via numerical chemical models, adopting the Larson-Penston solution and its delayed analogues to study collapse. Molecular abundances and distributions in a collapsing core are determined by the balance among the dynamical, chemical and adsorption time scales. When the central density n_H of a prestellar core with the Larson-Penston flow rises to 3 10^6 cm^{-3}, the CCS and CO column densities are calculated to show central holes of radius 7000 AU and 4000 AU, respectively, while the column density of N2H+ is centrally peaked. These predictions are consistent with observations of L1544. If the dynamical time scale of the core is larger than that of the Larson-Penston solution owing to magnetic fields, rotation, or turbulence, the column densities of CO and CCS are smaller, and their holes are larger than in the Larson-Penston core with the same central gas density. On the other hand, N2H+ and NH3 are more abundant in the more slowly collapsing core. Therefore, molecular distributions can probe the collapse time scale of prestellar cores. Deuterium fractionation has also been studied via numerical calculations. The deuterium fraction in molecules increases as a core evolves and molecular depletion onto grains proceeds. When the central density of the core is n_H=3 10^6 cm^{-3}, the ratio DCO+/HCO+ at the center is in the range 0.06-0.27, depending on the collapse time scale and adsorption energy; this range is in reasonable agreement with the observed value in L1544.

Millimeter wavelength observations are presented of NGC 1333 IRAS 4, a group of highly-embedded young stellar objects in Perseus, that reveal motions of infall, outflow, rotation, and turbulence in the dense gas around its two brightest continuum objects, 4A and 4B. These data have finest angular resolution of approximately 2" (0.0034 pc) and finest velocity resolution of 0.13 km/s. Infall motions are seen from inverse P-Cygni profiles observed in H2CO 3_12-2_11 toward both objects, but also in CS 3-2 and N2H+ 1-0 toward 4A, providing the least ambiguous evidence for such motions toward low-mass protostellar objects. Outflow motions are probed by bright line wings of H2CO 3_12-2_11 and CS 3-2 observed at positions offset from 4A and 4B, likely tracing dense cavity walls. Rotational motions of dense gas are traced by a systematic variation of the N2H+ line velocities, and such variations are found around 4A but not around 4B. Turbulent motions appear reduced with scale, given N2H+ line widths around both 4A and 4B that are narrower by factors of 2 or 3 than those seen from single-dish observations. Minimum observed line widths of approximately 0.2 km/s provide a new low, upper bound to the velocity dispersion of the parent core to IRAS 4, and demonstrate that turbulence within regions of clustered star formation can be reduced significantly. A third continuum object in the region, 4B', shows no detectable line emission in any of the observed molecular species.

We present the aperture synthesis imaging of C3H2 2_{12}-1_{01} and HCO+ 1-0 lines and of continuum emission at 3.4 mm toward the class 0 young stellar object, IRAS 18148-0440 in the L483 molecular cloud. The continuum emission is detected at the IRAS position at a level of 16 mJy, indicating a compact source with a 0.13 solar mass. The C3H2 delineates an envelope near the IRAS position, with a size of 3000 AU x 2000 AU in R.A. and Dec. directions, respectively. A velocity gradient detected only along the axis of outflow can be explained in terms of free-fall motion of the envelope. The HCO+ line wing extends up to the velocity of +- 6 km s^{-1} relative to the systemic one, and the high velocity material shows a symmetric bipolarity and extends over 2' or 0.1 pc along the east-west direction. It is found that the outflow material is clumpy and the opening angle is widest for the slowest moving component. The core component of the HCO+ 1-0 line exhibits an anti-infall asymmetry not only in interferometric but also in single dish observations. It is attributable to the slow isotropic outward motion of gas between the flattened envelope and the collimated outflow.

We have carried out observations of the starless core L1544 in the CCS (J_N=3_2-2_1) line at 9 millimeters wavelength using the BIMA array. The maps show an elongated condensation, 0.15 x 0.045 pc in size, with stronger emission at the edges. The appearance is consistent with a flattened, ringlike structure viewed at high inclination to the line of sight. The CCS molecule is likely heavily depleted in the inner part of the core. The position velocity diagram along the major axis shows a remarkable pattern, a "tilted ellipse", that can be reproduced by a simple model ring with motions of both infall and rotation. The models suggest comparable velocities for infall and rotation, ~0.1 km/s, in the outermost envelope, at radius 15000 AU.