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 mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray-tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass $M_{\nu} =$ 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of $M_{\nu}$. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

Context: The protostellar envelopes, outflow and large-scale chemistry of Class~0 and Class~I objects have been well-studied, but while previous works have hinted at or found a few Keplerian disks at the Class~0 stage, it remains to be seen if their presence in this early stage is the norm. Likewise, while complex organics have been detected toward some Class~0 objects, their distribution is unknown as they could reside in the hottest parts of the envelope, in the emerging disk itself or in other components of the protostellar system, such as shocked regions related to outflows. Aims: In this work, we aim to address two related issues regarding protostars: when rotationally supported disks form around deeply embedded protostars and where complex organic molecules reside in such objects. Methods: We observed the deeply embedded protostar, L483, using Atacama Large Millimeter/submillimeter Array (ALMA) Band~7 data from Cycles~1 and 3 with a high angular resolution down to $\sim$~0.1$^{\prime\prime}$ (20~au) scales. Results: We find that the kinematics of CS~$J=7$--$6$ and H$^{13}$CN~$J=4$--$3$ are best fitted by the velocity profile from infall under conservation of angular momentum and not by a Keplerian profile. The spatial extents of the observed complex organics are consistent with an estimated ice sublimation radius of the envelope at $\sim$~50~au, suggesting that the complex organics exist in the hot corino of L483. Conclusions: We find that L483 does not harbor a Keplerian disk down to at least $15$~au in radius. Instead, the innermost regions of L483 are undergoing a rotating collapse. This result highlights that some Class~0 objects contain only very small disks, or none at all, with the complex organic chemistry taking place on scales inside the hot corino of the envelope, in a region larger than the emerging disk.

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 report the first high-resolution (sub-arcminute) large-scale mapping $^{12}$CO and $^{13}$CO observations of the molecular clouds associated with the giant outer Galaxy HII region CTB~102 (KR 1). These observations were made using a newly commissioned receiver system on the 13.7-m radio telescope at the Taeduk Radio Astronomy Observatory. Our observations show that the molecular clouds have a spatial extent of $60 \times 35$ pc and a total mass of $10^{4.8} - 10^{5.0}$ M$_\odot$. Infrared data from WISE and 2MASS were used to identify and classify the YSO population associated with ongoing star formation activity within the molecular clouds. We directly detect 18 class I/class II YSOs and six transition disk objects. Moving away from the HII region, there is an age/class gradient consistent with sequential star formation. The infrared and molecular-line data were combined to estimate the star formation efficiency (SFE) of the entire cloud as well as the SFE for various sub-regions of the cloud. We find that the overall SFE is between $\sim5 - 10$%, consistent with previous observations of giant molecular clouds. One of the sub-regions, region 1a, is a clear outlier, with a SFE of 17 $-$ 35% on a 5 pc spatial scale. This high SFE is more typical for much smaller (sub-pc scale) star-forming cores, and we think region 1a is likely an embedded massive protocluster.

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 present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

The HIgh-Resolution Mid-infrarEd Spectrometer (HIRMES) is the 3rd Generation Instrument for the Stratospheric Observatory For Infrared Astronomy (SOFIA), currently in development at the NASA Goddard Space Flight Center (GSFC), and due for commissioning in 2019. By combining direct-detection Transition Edge Sensor (TES) bolometer arrays, grating-dispersive spectroscopy, and a host of Fabry-Perot tunable filters, HIRMES will provide the ability for High Resolution (R~100,000), Mid-Resolution (R~10,000), and Low-Resolution (R~600) slit-spectroscopy, and 2D Spectral Imaging (R~2000 at selected wavelengths) over the 25 - 122 {\mu}m mid-far infrared waveband. The driving science application is the evolution of proto-planetary systems via measurements of water-vapor, water-ice, deuterated hydrogen (HD), and neutral oxygen lines. However, HIRMES has been designed to be as flexible as possible to cover a wide range of science cases that fall within its phase-space, all whilst reaching sensitivities and observing powers not yet seen thus far on SOFIA, providing unique observing capabilities which will remain unmatched for decades.

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 performed a molecular line survey of 82 high-mass protostellar objects in search for inflow signatures associated with high-mass star formation. Using the H$^{13}$CO$^+$ (1$-$0) line as an optically thin tracer, we detected a statistically significant excess of blue asymmetric line profiles in the HCO$^+$ (1$-$0) transition, but nonsignificant excesses in the HCO$^+$ (3$-$2) and H$_2$CO (2$_{12}$$-$1$_{11}$) transitions. The negative blue excess for the HCN (3$-$2) transition suggests that the line profiles are affected by dynamics other than inflow motion. The HCO$^+$ (1$-$0) transition thus seems to be the suitable tracer of inflow motions in high-mass star-forming regions, as previously suggested. We found 27 inflow candidates that have at least one blue asymmetric profile and no red asymmetric profile, and derived the inflow velocities to be 0.23$-$2.00 km s$^{-1}$ for 20 of them using a simple two-layer radiative transfer model. Our sample is divided into two groups in different evolutionary stages. The blue excess of the group in relatively earlier evolutionary stages was estimated to be slightly higher than that of the other in the HCO$^+$ (1$-$0) transition.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

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.

We have investigated the TUKH122 prestellar core in the Orion A cloud using ALMA 3 mm dust continuum, N$_2$H$^+$ ($J=1-0$), and CH$_3$OH ($J_K=2_K-1_K$) molecular line observations. Previous studies showed that TUKH122 is likely on the verge of star formation because the turbulence is almost dissipated and chemically evolved among other starless cores in the Orion A cloud. By combining ALMA 12-m and ACA data, we recover extended emission with a resolution of $\sim5"$ corresponding to 0.01 pc and identify 6 condensations with a mass range of $0.1-0.4$ $M_\odot$ and a radius of $\lesssim0.01$ pc. These condensations are gravitationally bound following a virial analysis and are embedded in the filament including the elongated core with a mass of $\sim29$ $M_\odot$ and a radial density profile of $r^{-1.6}$ derived by {\it Herschel}. The separation of these condensations is $\sim0.035$ pc, consistent with the thermal jeans length at a density of $4.4\times10^5$ cm$^{-3}$. This density is similar to the central part of the core. We also find a tendency that the N$_2$H$^+$ molecule seems to deplete at the dust peak condensation. This condensation may be beginning to collapse because the linewidth becomes broader. Therefore, the fragmentation still occurs in the prestellar core by thermal Jeans instability and multiple stars are formed within the TUKH122 prestellar core. The CH$_3$OH emission shows a large shell-like distribution and surrounds these condensations, suggesting that the CH$_3$OH molecule formed on dust grains is released into gas phase by non-thermal desorption such as photoevaporation caused by cosmic-ray induced UV radiation.

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

Understanding the formation of wide binary systems of very low mass stars (M $\le$ 0.1 Msun) is challenging. The most obvious route is via widely separated low-mass collapsing fragments produced through turbulent fragmentation of a molecular core. However, close binaries/multiples from disk fragmentation can also evolve to wide binaries over a few initial crossing times of the stellar cluster through tidal evolution. Finding an isolated low mass wide binary system in the earliest stage of formation, before tidal evolution could occur, would prove that turbulent fragmentation is a viable mechanism for (very) low mass wide binaries. Here we report high resolution ALMA observations of a known wide-separation protostellar binary, showing that each component has a circumstellar disk. The system is too young to have evolved from a close binary and the disk axes are misaligned, providing strong support for the turbulent fragmentation model. Masses of both stars are derived from the Keplerian rotation of the disks; both are very low mass stars.

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 observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

We report the detection of a wing component in NH$_3$ emission toward the starless core TUKH122 with subthermal turbulence in the Orion A cloud. This NH$_3$ core is suggested to be on the verge of star formation because the turbulence inside the NH$_3$ core is almost completely dissipated, and also because it is surrounded by CCS, which resembles the prestellar core L1544 in Taurus showing infall motions. Observations were carried out with the Nobeyama 45 m telescope at 0.05 km s$^{-1}$ velocity resolution. We find that the NH$_3$ line profile consists of two components. The quiescent main component has a small linewidth of 0.3 km s$^{-1}$ dominated by thermal motions, and the red-shifted wing component has a large linewidth of 1.36 km s$^{-1}$ representing turbulent motions. These components show kinetic temperatures of 11 K and $<$ 30 K, respectively. Furthermore, there is a clear velocity offset between the NH$_3$ quiescent gas ($VLSR=3.7$ km s$^{-1}$) and the turbulent gas ($VLSR=4.4$ km s$^{-1}$). The centroid velocity of the turbulent gas corresponds to that of the surrounding gas traced by the $^{13}$CO ($J=1-0$) and CS ($J=2-1$) lines. LVG model calculations for CS and CO show that the turbulent gas has a temperature of $8-13$ K and an H$_2$ density of ~ $10^4$ cm$^{-3}$, suggesting that the temperature of the turbulent component is also ~ 10 K. The detections of both NH$_3$ quiescent and wing components may indicate a sharp transition from the turbulent parent cloud to the quiescent dense core.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

We present observations of HDCO and H2CO emission toward a sample of 15 Class 0 protostars in the Orion A and B clouds. Of these, eleven protostars are Herschel-identified PACS Bright Red Sources (PBRS) and four are previously identified protostars. Our observations revealed the chemical properties of the PBRS envelope for the first time. The column densities of HDCO and H2CO are derived from single dish observations at an angular resolution of ~20 arcsec (~8400 AU). The degree of deuteration in H2CO ([HDCO]/[H2CO]) was estimated to range from 0.03 to 0.31. The deuterium fractionation of most PBRS (70%) is similar to that of the non-PBRS sources. Three PBRS (30%) exhibit high deuterium fractionation, larger than 0.15. The large variation of the deuterium fractionation of H2CO in the whole PBRS sample may reflect the diversity in the initial conditions of star forming cores. There is no clear correlation between the [HDCO]/[H2CO] ratio and the evolutionary sequence of protostars.

Observations of the isolated globule B335 with ALMA have yielded absorption features against the continuum that are redshifted from the systemic velocity in both HCN and HCO$^+$ lines. These features provide unambiguous evidence for infall toward a central luminosity source. Previously developed models of inside-out collapse can match the observed line profiles of HCN and HCO$^+$ averaged over the central 50 AU. At the new distance of 100 pc, the inferred infall radius is 0.012 pc, the mass infall rate is $3 \times 10^{-6}$ solar masses per year, the age is 50,000 years, and the accumulated mass in the central zone is 0.15 solar masses, most of which must be in the star or in parts of a disk that are opaque at 0.8 mm. The continuum detection indicates an optically thin mass (gas and dust) of only $7.5\times 10^{-4}$ solar masses in the central region, consistent with only a very small disk mass.

We have observed the HC3N (J=10-9) and N2H+ (J=1-0) lines toward the Vela C molecular clouds with the Mopra 22 m telescope to study chemical characteristics of dense cores. The intensity distributions of these molecules are similar to each other at an angular resolution of 53", corresponding to 0.19 pc suggesting that these molecules trace the same dense cores. We identified 25 local peaks in the velocity-integrated intensity maps of the HC3N and/or N2H+ emission. Assuming LTE conditions, we calculated the column densities of these molecules and found a tendency that N2H+/HC3N abundance ratio seems to be low in starless regions while it seems to be high in star-forming regions, similar to the tendencies in the NH3/CCS, NH3/HC3N, and N2H+/CCS abundance ratios found in previous studies of dark clouds and the Orion A GMC. We suggest that carbon chain molecules, including HC3N, may trace chemically young molecular gas and N-bearing molecules, such as N2H+, may trace later stages of chemical evolution in the Vela C molecular clouds. It may be possible that the N2H+/HC3N abundance ratio of ~ 1.4 divides the star-forming and starless peaks in the Vela C, although it is not as clear as those in NH3/CCS, NH3/HC3N, and N2H+/CCS for the Orion A GMC. This less clear separation may be caused by our lower spatial resolution or the misclassification of star-forming and starless peaks due to the larger distance of the Vela C. It might be also possible that the HC3N (J=10-9) transition is not a good chemical evolution tracer compared with CCS (J=4-3 and 7-6) transitions.

Star formation activities in the NGC 2024 FIR 4 region were studied by imaging centimeter continuum sources and water maser sources using several archival data sets from the Very Large Array. The continuum source VLA 9 is elongated in the northwest-southeast direction, consistent with the FIR 4 bipolar outflow axis, and has a flat spectrum in the 6.2-3.6 cm interval. The three water maser spots associated with FIR 4 are also distributed along the outflow axis. One of the spots is located close to VLA 9, and another one is close to an X-ray source. Examinations of the positions of compact objects in this region suggest that the FIR 4 cloud core contains a single low-mass protostar. VLA 9 is the best indicator of the protostellar position. VLA 9 may be a radio thermal jet driven by this protostar, and it is unlikely that FIR 4 contains a high-mass young stellar object (YSO). A methanol 6.7 GHz maser source is located close to VLA 9, at a distance of about 100 AU. The FIR 4 protostar must be responsible for the methanol maser action, which suggests that methanol class II masers are not necessarily excited by high-mass YSOs. Also discussed are properties of other centimeter continuum sources in the field of view and the water masers associated with FIR 6n. Some of the continuum sources are radio thermal jets, and some are magnetically active young stars.

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

Ten very low luminosity objects were observed multiple times in the 8.5 GHz continuum in search of protostellar magnetic activities. A radio outburst of IRAM 04191+1522 IRS was detected, and the variability timescale was about 20 days or shorter. The results of this survey and archival observations suggest that IRAM 04191+1522 IRS is in active states about half the time. Archival data show that L1014 IRS and L1148 IRS were detectable previously and suggest that at least 20%-30% of very low luminosity protostars are radio variables. Considering the variability timescale and flux level of IRAM 04191+1522 IRS and the previous detection of the circular polarization of L1014 IRS, the radio outbursts of these protostars are probably caused by magnetic flares. However, IRAM 04191+1522 IRS is too young and small to develop an internal convective dynamo. If the detected radio emission is indeed coming from magnetic flares, the discovery implies that the flares may be caused by the fossil magnetic fields of interstellar origin.

We imaged two starless molecular cloud cores, TUKH083 and TUKH122, in the Orion A giant molecular cloud in the CCS and ammonia (NH$_3$) emission with the Very Large Array. TUKH122 contains one NH$_3$ core "TUKH122-n," which is elongated and has a smooth oval boundary. Where observed, the CCS emission surrounds the NH$_3$ core. This configuration resembles that of the N$_2$H$^+$ and CCS distribution in the Taurus starless core L1544, a well-studied example of a dense prestellar core exhibiting infall motions. The linewidth of TUKH122-n is narrow (0.20 km s$^{-1}$) in the NH$_3$ emission line and therefore dominated by thermal motions. The smooth oval shape of the core boundary and narrow linewidth in NH$_3$ seem to imply that TUKH122-n is dynamically relaxed and quiescent. TUKH122-n is similar to L1544 in the kinetic temperature (10 K), linear size (0.03 pc), and virial mass ($\sim$ 2 $M_{\odot}$). Our results strongly suggest that TUKH122-n is on the verge of star formation. TUKH122-n is embedded in the 0.2 pc massive (virial mass $\sim$ 30 $M_{\odot}$) turbulent parent core, while the L1544 NH$_3$ core is embedded in the 0.2 pc less-massive (virial mass $\sim$ 10 $M_{\odot}$) thermal parent core. TUKH083 shows complicated distribution in NH$_3$, but was not detected in CCS. The CCS emission toward TUKH083 appears to be extended, and is resolved out in our interferometric observations.

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

The results of a maser survey toward ninety-nine protostars in the Orion molecular cloud complex are presented. The target sources are low-mass protostars identified from infrared observations. Single-dish observations were carried out in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. Most of the detected sources were mapped to determine the source positions. Five water maser sources were detected, and they are excited by HH 1-2 VLA 3, HH 1-2 VLA 1, L1641N MM1/3, NGC 2071 IRS 1/3, and an object in the OMC 3 region. The water masers showed significant variability in intensity and velocity with time scales of a month or shorter. Four methanol emission sources were detected, and those in the OMC 2 FIR 3/4 and L1641N MM1/3 regions are probably masers. The methanol emission from the other two sources in the NGC 2071 IRS 1-3 and V380 Ori NE regions are probably thermal. For the water masers, the number of detections per protostar in the survey region is about 2%, which suggests that the water masers of low-mass protostars are rarely detectable. The methanol class I maser of low-mass protostars is an even rarer phenomenon, with a detection rate much smaller than 1%.

We have mapped six molecular cloud cores in the Orion A giant molecular cloud (GMC), whose kinetic temperatures range from 10 to 30 K, in CCS and N2H+ with Nobeyama 45 m radio telescope to study their chemical characteristics. We identified 31 intensity peaks in the CCS and N2H+ emission in these molecular cloud cores. It is found for cores with temperatures lower than ~ 25 K that the column density ratio of N(N2H+)/N(CCS) is low toward starless core regions while it is high toward star-forming core regions, in case that we detected both of the CCS and N2H+ emission. This is very similar to the tendency found in dark clouds (kinetic temperature ~ 10 K). The criterion found in the Orion A GMC is N(N2H+)/N(CCS) ~ 2-3. In some cases, the CCS emission is detected toward protostars as well as the N2H+ emission. Secondary late-stage CCS peak in the chemical evolution caused by CO depletion may be a possible explanation for this. We found that the chemical variation of CCS and N2H+ can also be used as a tracer of evolution in warm (10-25 K) GMC cores. On the other hand, some protostars do not accompany N2H+ intensity peaks but are associated with dust continuum emitting regions, suggesting that the N2H+ abundance might be decreased due to CO evaporation in warmer star-forming sites.

We present Herschel/PACS observations of L1448-MM, a Class 0 protostar with a prominent outflow. Numerous emission lines are detected at 55 < lambda < 210 micrometer including CO, OH, H2O, and [O I]. We investigate the spatial distribution of each transition to find that lines from low energy levels tend to distribute along the outflow direction while lines from high energy levels peak at the central spatial pixel. Spatial maps reveal that OH emission lines are formed in a relatively small area, while [O I] emission is extended. According to the rotational diagram analysis, the CO emission can be fitted by two (warm and hot) temperature components. For H2O, the ortho-to-para ratio is close to 3. The non-LTE LVG calculations suggest that CO and H2O lines could instead be formed in a high kinetic temperature (T > 1000 K) environment, indicative of a shock origin. For OH, IR-pumping processes play an important role in the level population. The molecular emission in L1448-MM is better explained with a C-shock model, but the atomic emission of PACS [O I] and Spitzer/IRS [Si II] emission is not consistent with C-shocks, suggesting multiple shocks in this region. Water is the major line coolant of L1448-MM in the PACS wavelength range, and the best-fit LVG models predict that H2O and CO emit (50-80)% of their line luminosity in the PACS wavelength range.

The NGC 2024 FIR 6 region was observed in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. The water maser spectra displayed several velocity components and month-scale time variabilities. Most of the velocity components may be associated with FIR 6n, while one component was associated with FIR 4. A typical life time of the water-maser velocity-components is about 8 months. The components showed velocity fluctuations with a typical drift rate of about 0.01 km/s/day. The methanol class I masers were detected toward FIR 6. The methanol emission is confined within a narrow range around the systemic velocity of the FIR 6 cloud core. The methanol masers suggest the existence of shocks driven by either the expanding H II region of FIR 6c or the outflow of FIR 6n.

The NGC 2024 FIR 5/6 region was observed in the 6.9 mm continuum with an angular resolution of about 1.5 arcsec. The 6.9 mm continuum map shows four compact sources, FIR 5w, 5e, 6c, and 6n, as well as an extended structure of the ionization front associated with the optical nebulosity. FIR 6c has a source size of about 0.4 arcsec or 150 AU. The spectral energy distribution (SED) of FIR 6c is peculiar: rising steeply around 6.9 mm and flat around 1 mm. The possibility of a hypercompact H II region is explored. If the millimeter flux of FIR 6c comes from hot ionized gas heated by a single object at the center, the central object may be a B1 star of about 5800 solar luminosities and about 13 solar masses. The 6.9 mm continuum of FIR 6n may be a mixture of free-free emission and dust continuum emission. Archival data show that both FIR 6n and 6c exhibit water maser activity, suggesting the existence of shocked gas around them. The 6.9 mm continuum emission from FIR 5w has a size of about 1.8 arcsec or 760 AU. The SEDs suggest that the 6.9 mm emission of FIR 5w and 5e comes from dust, and the masses of the dense molecular gas are about 0.6 and 0.5 solar masses, respectively.

We obtained 350 and 850 micron continuum maps of the infrared dark cloud G049.40-00.01. Twenty-one dense clumps were identified within G049.40-00.01 based on the 350 micron continuum map with an angular resolution of about 9.6". We present submillimeter continuum maps and report physical properties of the clumps. The masses of clumps range from 50 to 600 M_sun. About 70% of the clumps are associated with bright 24 micron emission sources, and they may contain protostars. The most massive two clumps show extended, enhanced 4.5 micron emission indicating vigorous star-forming activity. The clump size-mass distribution suggests that many of them are forming high mass stars. G049.40-00.01 contains numerous objects in various evolutionary stages of star formation, from pre-protostellar clumps to HII regions.

The NGC 1333 IRAS 4B region was observed in the 6.9 mm and 1.3 cm continuum with an angular resolution of about 0.4 arcseconds. IRAS 4BI was detected in both bands, and BII was detected in the 6.9 mm continuum only. The 1.3 cm source of BI seems to be a disk-like flattened structure with a size of about 50 AU. IRAS 4BI does not show any sign of multiplicity. Examinations of archival infrared images show that the dominating emission feature in this region is a bright peak in the southern outflow driven by BI, corresponding to the molecular hydrogen emission source HL 9a. Both BI and BII are undetectable in the mid-IR bands. The upper limit on the far-IR flux of IRAS 4BII suggests that it may be a very low luminosity young stellar object.

The NGC 1333 IRAS 4A protobinary was observed in the 1.3 cm and 6.9 mm continuum and the ammonia and SiO lines, with an angular resolution of about 0.4 arcseconds. The continuum maps show the circumstellar structures of the two protostars, A1 and A2. The A1 system is brighter and more massive than the A2 system. The ratio of mass, including dense gas and protostar, is about 6. The properties of the circumstellar disks and outflows suggest that A1 may be younger than A2. The deflected part of the northeastern jet of A2 is bright in the SiO line, and the distance between the brightest peak and deflection point suggests that the enhancement of SiO takes about 100 yr after the collision with a dense core. The ammonia maps show a small structure that seems to be a part of the obstructing core. The outflow properties were studied by comparing interferometric maps of SiO, ammonia, formaldehyde, and HCN lines. Their overall structures agree well, suggesting that these species are excited by the same mechanism. However, the intensity distributions show that SiO is chemically unique. SiO may be directly linked to the primary jet while the other species may be tracing the entrained ambient gas.

The bipolar jet of the NGC 1333 IRAS 4A2 protostar shows a velocity gradient in the direction perpendicular to the jet axis. This lateral velocity gradient can be seen throughout the jet imaged in a silicon monoxide line, 2500-8700 AU from the driving source, and is consistent with the rotation of the accretion disk. If this gradient is caused by the rotation of the jet around its axis, the average specific angular momentum is about 1.5 x 10^21 cm^2 s^-1. Comparison of the kinematics between the jet and the disk suggests that the jet-launching region on the disk has a radius of about 2 AU, which supports the disk-wind models. The angular momentum transported away by the jet seems to be large enough for the protostar to accrete matter from the disk, confirming the crucial role of jets in the early phase of star formation process.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines with an angular resolution of 0.3 arcsec. The ammonia emission source of IRAS 4A2 is elongated in the direction perpendicular to the bipolar jet and has a size of 0.55 arcsec or 130 AU. This emission structure was interpreted as a circumstellar disk around the IRAS 4A2 protostar, and the rotation kinematics of the disk was investigated by making a position-velocity diagram along the major axis. Assuming a Keplerian rotation, the disk has a rotation velocity of 1.8 km/s at a radius of 20 AU, which implies a central object of about 0.08 solar masses. The collapse age of the protostar is about 50,000 yr. The mass, accretion rate, and age are consistent with what are expected from the standard theory of low-mass star formation. If IRAS 4A2 grows at this rate, it may become a star similar to the Sun.

We present the results of an extensive observational study of the active star-forming complex W51 that was observed in the J=2-1 transition of the 12CO and 13CO molecules over a 1.25 deg x 1.00 deg region with the University of Arizona Heinrich Hertz Submillimeter Telescope. We use a statistical equilibrium code to estimate physical properties of the molecular gas. We compare the molecular cloud morphology with the distribution of infrared (IR) and radio continuum sources, and find associations between molecular clouds and young stellar objects (YSOs) listed in Spitzer IR catalogs. The ratios of CO lines associated with HII regions are different from the ratios outside the active star-forming regions. We present evidence of star formation triggered by the expansion of the HII regions and by cloud-cloud collisions. We estimate that about 1% of the cloud mass is currently in YSOs.

The HH 1-2 region in the L1641 molecular cloud was observed in the near-IR J, H, and Ks bands, and imaging polarimetry was performed. Seventy six point-like sources were detected in all three bands. The near-IR polarizations of these sources seem to be caused mostly by the dichroic extinction. Using a color-color diagram, reddened sources with little infrared excess were selected to trace the magnetic field structure of the molecular cloud. The mean polarization position angle of these sources is about 111 deg, which is interpreted as the projected direction of the magnetic field in the observed region of the cloud. The distribution of the polarization angle has a dispersion of about 11 deg, which is smaller than what was measured in previous studies. This small dispersion gives a rough estimate of the strength of the magnetic field to be about 130 microG and suggests that the global magnetic field in this region is quite regular and straight. In contrast, the outflows driven by young stellar objects in this region seem to have no preferred orientation. This discrepancy suggests that the magnetic field in the L1641 molecular cloud does not dictate the orientation of the protostars forming inside.

The Serpens SMM 1 region was observed in the 6.9 mm continuum with an angular resolution of about 0.6 arcsec. Two sources were found to have steep positive spectra suggesting emission from dust. The stronger one, SMM 1a, is the driving source of the bipolar jet known previously, and the mass of the dense molecular gas traced by the millimeter continuum is about 8 solar mass. The newly found source, SMM 1b, positionally coincides with the brightest mid-IR source in this region, which implies that SMM 1b is yet another young stellar object. SMM 1b seems to be less deeply embedded than SMM 1a. SMM 1 is probably a protobinary system with a projected separation of 500 AU.

Circularly polarized 3.5 cm continuum emission was detected toward three radio sources in the R CrA region using the Very Large Array. The Class I protostar IRS 5b persistently showed polarized radio emission with a constant helicity over 8 yr, which suggests that its magnetosphere has a stable configuration. There is a good correlation between the Stokes I and Stokes V fluxes, and the fractional polarization is about 0.17. During active phases the fractional polarization is a weakly decreasing function of Stokes I flux, which suggests that IRS 5b is phenomenologically similar to other types of flare stars such as RS CVn binaries. The variability timescale of the polarized flux is about a month, and the magnetosphere of IRS 5b must be very large in size. The Class I protostar IRS 7A was detected once in circularly polarized radio emission, even though IRS 7A drives a thermal radio jet. This detection implies that the radio emission from the magnetosphere of a young protostar can escape the absorption by the partially ionized wind at least once in a while. The properties of IRS 7A and IRS 5b suggests that Class I protostars have organized peristellar magnetic fields of a few kilogauss and that the detectability of magnetospheric emission may depend on the evolutionary status of protostar. Also reported is the detection of circularly polarized radio emission toward the variable radio source B5.

Using a sub-pixel event repositioning technique, we spatially resolved X-ray emission from the infrared double system IRS 5 in the R Corona Australis molecular cloud with ~0.8" separation. As far as we know, this result - obtained from 8 Chandra archival observations between 2000 and 2005 - is the first X-ray study of individual sources in a Class I protostar binary system with a projected separation of less than 200 AU. We extracted light curves and spectra of the individual sources using a two-dimensional image fitting method. IRS 5a at the south, the source which was brighter in the near-infrared, showed three X-ray flares lasting >20 ksec, reminiscent of X-ray flares from pre-main sequence stars, while the northern source (IRS 5b) was quiescent in X-rays in all the observations except for a 2005 August 9 observation with a factor of ~2 flux enhancement. In quiescence, these sources showed almost identical X-ray spectra, with NH ~4e22 cm-2, kT ~2 keV, and log Lx ~30.2-3 ergs s-1. IRS 5a showed plasma at temperatures up to kT ~5-6 keV during flares, while the column density of IRS 5b increased by a factor of 2 during an observation on 2005 August 9. We discuss the evolutionary stages and variation of the X-ray activity of these sources.

The R CrA region was observed in the 3.5 and 6.2 cm continuum with high angular resolutions (0.6--1.7 arcseconds) using the Very Large Array. Archival data sets were also analyzed for comparison, which provided angular resolutions up to 0.3 arcseconds. A cluster of young stellar objects was detected, and a rich array of star forming activities was revealed. IRS 7A showed an enhanced outflow activity recently. The main peak of IRS 7A positionally coincides with an X-ray source, which suggests that the X-ray emission is directly related to the central protostar. The Class 0 source SMA 2 is associated with a double radio source, B 9a and 9b, and seems to be driving two outflows. The B 9 complex is probably a multiple-protostar system. Both B 9a and 9b are nonthermal radio sources with negative spectral indices. IRS 7B is a compact radio source surrounded by an extended structure. The compact source corresponds to the Class 0/I source SMA 1, and it is also closely associated with an X-ray source, suggesting that magnetic activities start early in the protostellar stage of evolution. The extended structure of IRS 7B may be a bipolar outflow. IRS 5 was resolved into two sources with a separation of 0.9 arcseconds. Both IRS 5a and 5b display radio flares and X-ray emission, suggesting that energetic magnetic processes are active in both members. The month-scale active phase of IRS 5b implies that the flare activity must involve large-scale magnetic fields. During the strong flare event of IRS 5b in 1998, IRS 5a also showed an enhanced level of radio emission. This concurrent activity suggests that IRS 5 may be an interacting young binary system, but the interaction mechanism is unknown. Alternatively, what was seen in the radio images could be a circumbinary halo.

High-resolution SiO observations of the NGC 1333 IRAS 4A star-forming region showed a highly collimated outflow with a substantial deflection. The deflection was suggested to be caused by the interactions of the outflow and a dense cloud core. To investigate the deflection process of protostellar outflows, we have carried out three-dimensional hydrodynamic simulations of the collision of an outflow with a dense cloud. Assuming a power-law type density distribution of the obstructing cloud, the numerical experiments show that the deflection angle is mainly determined by the impact parameter and the density contrast between the outflow and the cloud. The deflection angle is, however, relatively insensitive to the velocity of the outflow. Using a numerical model with physical conditions that are particularly suitable for the IRAS 4A system, we produce a column-density image and a position-velocity diagram along the outflow, and they are consistent with the observations. Based on our numerical simulations, if we assume that the initial density and the velocity of the outflow are $\sim 10 \cm3$ and $\sim 100 \kms$, the densities of the dense core and ambient medium in the IRAS 4A system are most likely to be $\sim 10^5 \cm3$ and $\sim 10^2 \cm3$, respectively. We therefore demonstrate through numerical simulations that the directional variability of the IRAS 4A outflow can be explained reasonably well using the collision model.

The HH 24 MMS protostellar system was observed in the 6.9 mm continuum with a high angular resolution (0.5"). HH 24 MMS was resolved into two sources. The separation between sources 1 and 2 is ~0.9" or 360 AU. The spectral energy distribution suggests that the 6.9 mm flux is almost entirely from dust. The 6.9 mm image and the spectrum suggest that HH 24 MMS may be a protostellar binary system. Total mass including the accretion disks and the inner protostellar envelope is ~1.4 Msun. Disk masses of sources 1 and 2 are 0.6 and 0.3 Msun, respectively. Both sources are highly elongated. The difference in the position angle of the two disks is ~45 degree, which means that HH 24 MMS is a highly misaligned protobinary system. The misalignment suggests that turbulent fragmentation may be the formation mechanism relevant to the binary systems with a separation of a few hundreds of AU, such as the HH 24 MMS system.

We have obtained velocity-resolved spectra of the H_2 v=1-0 S(1) (lambda = 2.1218 micron) emission line at 2 arcsec angular resolution (or 0.08 pc spatial resolution) in four regions within the central 10 pc of the Galaxy where the supernova-like remnant Sgr A East is colliding with molecular clouds. To investigate the kinematic, physical, and positional relationships between the important gaseous components in the center, we compared the H_2 data cube with previously published NH_3 data. The projected interaction-boundary of Sgr A East is determined to be an ellipse with its center offset 1.5 pc from Sgr A* and dimensions of 10.8 pc X 7.6 pc. This H_2 boundary is larger than the synchrotron emission shell but consistent with the dust ring which is believed to trace the shock front of Sgr A East. Since Sgr A East is driving shocks into its nearby molecular clouds, we can determine their positional relationships using the shock directions as indicators. As a result, we suggest a revised model for the three-dimensional structure of the central 10 pc. The actual contact between Sgr A East and all of the surrounding molecular material, including the circum-nuclear disk and the southern streamer, makes the hypothesis of infall into the nucleus and feeding of Sgr A* very likely.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines and in the 1.3 cm continuum with a high resolution (about 1.0 arcsec). The ammonia maps show two compact sources, one for each protostar, and they are probably protostellar accretion disks. The disk associated with IRAS 4A2 is seen nearly edge-on and shows an indication of rotation. The A2 disk is brighter in the ammonia lines but dimmer in the dust continuum than its sibling disk, with the ammonia-to-dust flux ratios different by about an order of magnitude. This difference suggests that the twin disks have surprisingly dissimilar characters, one gas-rich and the other dusty. The A2 disk may be unusually active or hot, as indicated by its association with water vapor masers. The existence of two very dissimilar disks in a binary system suggests that the formation process of multiple systems has a controlling agent lacking in the isolated star formation process and that stars belonging to a multiple system do not necessarily evolve in phase with each other.

The protostar B335 was observed in the 1.3 mm continuum and in the H2CO 312 - 211 line with an angular resolution of about 8 arcsec. The mass of the inner envelope detected by the dust continuum emission is about 0.02 Msun. The H2CO spectrum at the protostellar position shows a blue-skewed double peak profile, suggesting that the kinematics of the inner envelope is dominated by infall motion. When the blueshifted and the redshifted peaks were imaged separately, however, there is a small east-west displacement between the maximum positions. This displacement suggests that some part of the H2CO emission might come from the outflowing gas. A combined effect of the infalling envelope and the outflow on the radiative transfer is discussed. This effect can make the line profile asymmetry severer than what is expected from infall-only models.

The NGC 1333 IRAS 4 region was observed in the 22 GHz H2O maser line with an angular resolution of about 0.08 arcseconds. Two groups of masers were detected, one near IRAS 4A and the other near BI. Among the eight maser spots detected near IRAS 4A, six spots are located close to A2, within 100 AU, and the maser velocities are also near the systemic velocity of the cloud core. These masers are probably related with the circumstellar disk. Since there is no maser spot detected around A1, the star forming process is relatively more active in A2 than in A1. Four maser spots were detected near IRAS 4BI. Since most of them are distributed along a straight line in the direction of the outflow, BI masers are most likely related with the jet or outflow. The disk-outflow dichotomy of H2O masers is discussed briefly. No maser was detected near H2O(C), another maser source reported previously.

The NGC 1333 region was observed in the H2 1-0 S(1) line. The H2 images cover a 5' x 7' region around IRAS 4. Numerous H2 emission features were detected. The northeast-southwest bipolar outflow driven by IRAS 4A was studied by combining the H2 images with SiO maps published previously. The SiO-H2 outflows are continuous on the southwestern side but show a gap on the northeastern side. The southwestern outflow lobe curves smoothly, and the position angle increases with the distance from the driving source. The base and the outer tip of the northeastern outflow lobe are located at positions opposite to the corresponding parts of the southwestern lobe. This point-symmetry suggests that the outflow axis may be drifting or precessing clockwise in the plane of the sky and that the cause of the axis drift may be intrinsic to the outflow engine. The axis drift model is supported by the asymmetric lateral intensity profile of the SiO outflow. The axis drift rate is about 0.011 deg yr-1. The middle part of the northeastern outflow does not exactly follow the point symmetry because of the superposition of two different kinds of directional variability: the axis drift of the driving source and the deflection by a dense core. The axis drift model provides a good explanation for the large deflection angle of the northeastern outflow. Other H2 emission features around the IRAS 4 region are discussed briefly. Some of them are newly found outflows, and some are associated with outflows already known before.

We present observations of 25 transitions of 17 isotopologues of 9 molecules toward B335. With a goal of constraining chemical models of collapsing clouds, we compare our observations, along with data from the literature, to models of chemical abundances. The observed lines are simulated with a Monte Carlo code, which uses various physical models of density and velocity as a function of radius. The dust temperature as a function of radius is calculated self-consistently by a radiative transfer code. The gas temperature is then calculated at each radius, including gas-dust collisions, cosmic rays, photoelectric heating, and molecular cooling. The results provide the input to the Monte Carlo code. We consider both {\it ad hoc} step function models for chemical abundances and abundances taken from a self-consistent modeling of the evolution of a star-forming core. The step function models can match the observed lines reasonably well, but they require very unlikely combinations of radial variations in chemical abundances. Among the self-consistent chemical models, the observed lines are matched best by models with somewhat enhanced cosmic-ray ionization rates and sulfur abundances. We discuss briefly the steps needed to close the loop on the modeling of dust and gas, including off-center spectra of molecular lines.

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

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 observed 23 Class 0 sources in the HCO+ J=4-3 and 3-2 lines. The mean bolometric temperature of the 16 sources with well-determined values is 44 K and the mean luminosity is 5.7 L_sun, excluding two sources of considerably higher luminosity. Nine sources, including three sources previously suggested to be collapsing, have the correct (blue) spectral line asymmetry for infall in both lines. Three sources have the opposite (red) asymmetry in both lines, and one source, L1157, has a red asymmetry in HCO+ J=4-3 and a blue asymmetry in J=3-2. The rest have no significant or consistent asymmetry. The H13CO+ J=4-3 and 3-2 lines were also observed to find the velocity of the ambient gas, and sources with an interesting line asymmetry were mapped. A Monte Carlo code was used to produce an evolutionary sequence of collapsing cloud models of the HCO+ J=4-3 and 3-2 lines and to compare various diagnostics of the resulting line profiles. The same code was used to compare infall models to the observations in one source, L1527. The results were consistent with previous collapse models. Based on integrated intensity maps of the line peaks and wings, as well as the velocity of the H13CO+ line, we select six of the nine sources with a blue line asymmetry as good candidates for protostellar collapse. Further evidence is needed to establish that infall is taking place. The HCO+ spectra are not conclusive because bipolar outflows produce strong emission, which can confuse the issue in any individual source. However, the predominance of blue asymmetries over red asymmetries is not naturally explained in outflow models, whereas it is expected in collapse models.

Extremely high velocity (EHV) wings, with full widths of 72 to 140 km/s, are seen on the CO J=3-2 lines toward W3 IRS 5, GL 490, NGC 2071, W28 A2, GL 2591, S140, and Cepheus A. The results of our survey suggest that EHV wings are common around infrared sources of moderate to high luminosity (500 to 4x10^5 Lsun) in dense regions. Line ratios imply that the EHV gas is usually optically thin and warm. Characteristic velocities range from 20 to 40 km/s, yielding timescales of 1600-4200 yr. Since most sources in this study are producing some ionizing photons, these short timescales suggest that neutral winds coexist with ionizing photons. We examined two possible sources for the EHV CO emission: a neutral stellar wind; and swept-up or entrained molecular gas. Neither can be ruled out. If the high-velocity (HV) gas is swept up by a momentum-conserving stellar wind traced by the extremely high velocity CO emission, most of the C in the winds from luminous objects cannot be in CO. If the EHV and HV forces are equal, the fraction of C in a form other than CO increases with source luminosity and with the production rate of ionizing photons.