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 study the prospects of detecting magnetic helicity in galaxies by observing the dust polarization of the edge-on galaxy NGC 891. Our numerical results of mean-field dynamo calculations show that there should be a large-scale component of the rotationally invariant parity-odd B polarization that we predict to be negative in the first and third quadrants, and positive in the second and fourth quadrants. The large-scale parity-even E polarization is predicted to be negative near the axis and positive further away in the outskirts. These properties are shown to be mostly a consequence of the magnetic field being azimuthal and the polarized intensity being maximum at the center of the galaxy and are not a signature of magnetic helicity.

Magnetic fields are ubiquitous in our Universe, but remain poorly understood in many branches of astrophysics. A key tool for inferring astrophysical magnetic field properties is dust emission polarimetry. The James Clerk Maxwell Telescope (JCMT) is planning a new 850$\mu$m camera consisting of an array of 7272 paired Microwave Kinetic Inductance Detectors (MKIDs), which will inherently acquire linear polarization information. The camera will allow wide-area polarization mapping of dust emission at 14$^{\prime\prime}$-resolution, allowing magnetic field properties to be studied in a wide range of environments, including all stages of the star formation process, Asymptotic Giant Branch stellar envelopes and planetary nebula, external galaxies including starburst galaxies and analogues for the Milky Way, and the environments of active galactic nuclei (AGN). Time domain studies of AGN and protostellar polarization variability will also become practicable. Studies of the polarization properties of the interstellar medium will also allow detailed investigation of dust grain properties and physics. These investigations would benefit from a potential future upgrade adding 450$\mu$m capability to the camera, which would allow inference of spectral indices for polarized dust emission in a range of environments. The enhanced mapping speed and polarization capabilities of the new camera will transform the JCMT into a true submillimetre polarization survey instrument, offering the potential to revolutionize our understanding of magnetic fields in the cold Universe.

The dependence of polarization fraction $p$ on total intensity $I$ in polarized submillimeter emission measurements is typically parameterized as $p\propto I^{-\alpha}$ $(\alpha \leq 1)$, and used to infer dust grain alignment efficiency in star-forming regions, with an index $\alpha=1$ indicating near-total lack of alignment of grains with the magnetic field. In this work we demonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurements of $\alpha \sim 1$ even in data with significant signal-to-noise in Stokes $Q$, $U$ and $I$ emission, and so with robust measurements of polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model is applicable, $\alpha$ can be recovered by fitting the $p-I$ relationship with the mean of the Rice distribution, without statistical debiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850$\mu$m observations of three clumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, $\alpha\approx 0.34$, while in the more isolated Oph B and C, despite their differing star formation histories, $\alpha \sim 0.6-0.7$. Our results thus suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field than by star formation history. We further find that grains may remain aligned with the magnetic field at significantly higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within star-forming clumps and cores.

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.

To shed light on the early phase of a low-mass protostar formation process, we conducted interferometric observations towards a protostar GF9-2 using the CARMA and SMA. The observations have been carried out in the CO J=3-2 line and in the continuum emission at the wavelengths of 3 mm, 1 mm and 850 micron. All the continuum images detected a single point-like source with a radius of 250+/-80 AU at the center of the previously known ~3 Msun molecular cloud core. A compact emission is detected towards the object at the Spitzer MIPS and IRAC bands as well as the four bands at the WISE. Our spectroscopic imaging of the CO line revealed that the continuum source is driving a 1000 AU scale molecular outflow, including a pair of lobes where a collimated "higher" velocity red lobe exists inside a poorly collimated "lower" velocity red lobe. These lobes are rather young and the least powerful ones so far detected. A protostellar mass of M~<0.06 Msun was estimated using an upper limit of the protostellar age of (4+/-1)x10^3 yrs and an inferred non-spherical steady mass accretion rate of ~10^{-5} Msun/yr. Together with results from an SED analysis, we discuss that the outflow system is driven by a protostar whose surface temperature of~3,000K, and that the natal cloud core is being dispersed by the outflow.

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 the first high-resolution, submillimeter-wavelength polarimetric observations of -- and thus direct observations of the magnetic field morphology within -- the dense gas of the Pillars of Creation in M16. These 850$\,\mu$m observations, taken as part of the BISTRO (B-Fields in Star-forming Region Observations) Survey using the POL-2 polarimeter on the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT), show that the magnetic field runs along the length of the pillars, perpendicular to, and decoupled from, the field in the surrounding photoionized cloud. Using the Chandrasekhar-Fermi method we estimate a plane-of-sky magnetic field strength of $170-320\,\mu$G in the Pillars, consistent with their having been formed through compression of gas with initially weak magnetization. The observed magnetic field strength and morphology suggests that the magnetic field may be slowing the pillars' evolution into cometary globules. We thus hypothesize that the evolution and lifetime of the Pillars may be strongly influenced by the strength of the coupling of their magnetic field to that of their parent photoionized cloud -- i.e. that the Pillars' longevity results from magnetic support

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 determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We report 3 mm continuum, CH3CN(5-4) and 13CS(2-1) line observations with CARMA, in conjunction with 6 and 1.3 cm continuum VLA data, and 12 and 25 micron broadband data from the Subaru Telescope toward the massive proto-star IRAS18566+0408. The VLA data resolve the ionized jet into 4 components aligned in the E-W direction. Radio components A, C, and D have flat cm SEDs indicative of optically thin emission from ionized gas, and component B has a spectral index alpha = 1.0, and a decreasing size with frequency proportional to frequency to the -0.5 power. Emission from the CARMA 3 mm continuum, and from the 13CS(2-1), and CH3CN(5-4) spectral lines is compact (i.e. < 6700 AU), and peaks near the position of VLA cm source, component B. Analysis of these lines indicates hot, and dense molecular gas, typical for HMCs. Our Subaru telescope observations detect a single compact source, coincident with radio component B, demonstrating that most of the energy in IRAS18566+0408 originates from a region of size < 2400 AU. We also present UKIRT near-infrared archival data for IRAS18566+0408 which show extended K-band emission along the jet direction. We detect an E-W velocity shift of about 10 km/sec over the HMC in the CH3CN lines possibly tracing the interface of the ionized jet with the surrounding core gas. Our data demonstrate the presence of an ionized jet at the base of the molecular outflow, and support the hypothesis that massive protostars with O-type luminosity form with a mechanism similar to lower 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.

Using the ALMA archival data of both CO(6--5) line and 689 GHz continuum emission towards the archetypical Seyfert galaxy, NGC 1068, we identified a distinct continuum peak separated by 14 pc from the nuclear radio component S1 in projection. The continuum flux gives a gas mass of ~2x10^5 Msun and bolometric luminosity of ~10^8 Lsun, leading to a star formation rate of ~0.1 Msun/yr. Subsequent analysis on the line data suggest that the gas has a size of ~10 pc, yielding to mean H2 number density of ~10^5 cm^{-3}. We therefore refer to the gas as "massive dense gas cloud": the gas density is high enough to form a "proto starcluster" whose stellar mass of ~10^4 Msun. We found that the gas stands a unique position between galactic and extraglactic clouds in the diagrams of start formation rate (SFR) vs. gas mass proposed by Lada et al. and surface density of gas vs. SFR density by Krumholz and McKee. All the gaseous and star-formation properties may be understood in terms of the turbulence-regulated star formation scenario. Since there are two stellar populations with the ages of 300 Myr and 30 Myr in the 100 pc-scale circumnulear region, we discuss that NGC1068 has experienced at least three episodic star formation events with a tendency that the inner star-forming region is the younger. Together with several lines of evidence that the dynamics of the nuclear region is decoupled from that of the entire galactic disk, we discuss that the gas inflow towards the nuclear region of NGC 1068 may be driven by a past minor merger.

NGC 2264-C is a high-mass protocluster where several star-formation events are known to have occurred. To investigate whether past protostellar activity has left a chemical imprint in this region, we mapped it in SiO($J = 2-1$), a shock tracer, and several other molecular lines with the Nobeyama 45 m telescope. Our observations show the presence of a complex network of protostellar outflows. The strongest SiO emission lies beyond a radius of $\sim 0.1$ pc with respect to the center of the clump, and is characterized by broad ($> 10$ km s$^{-1}$) lines and abundances of $\sim 1.4 \times 10^{-8}$ with respect to H$_2$. Interestingly, SiO appears relatively depleted ($\chi_\mathrm{SiO} \sim 4 \times 10^{-9}$) within this radius, despite it being affected by molecular outflow activity. We attribute this to fast condensation of SiO back onto dust grains and/or rapid gas-phase destruction of SiO, favored by the high density present in this area ($> 10^6$ cm$^{-3}$). Finally, we identify a peripheral, narrow-line ($\sim 2$ km s$^{-1}$) component, where SiO has an abundance of a few times 10$^{-11}$. After considering different options, we conclude that this weak emission may be tracing protostellar shocks from the star formation episode that preceded the current one, which have decelerated over time and eventually resulted in SiO being largely depleted/destroyed. Alternatively, a population of unresolved low-mass protostars may be responsible for the narrow SiO emission. High-angular resolution observations are necessary to distinguish between these two possibilities and thus understand the role of SiO as a chemical tracer of past star-formation episodes in massive protoclusters.

Spectral line survey observations are conducted toward the high-mass protostar candidate NGC 2264 CMM3 in the 4 mm, 3 mm, and 0.8 mm bands with the Nobeyama 45 m telescope and the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. In total, 265 emission lines are detected in the 4 mm and 3 mm bands, and 74 emission lines in the 0.8 mm band. As a result, 36 molecular species and 30 isotopologues are identified. In addition to the fundamental molecular species, many emission lines of carbon-chain molecules such as HC5N, C4H, CCS, and C3S are detected in the 4 mm and 3 mm bands. Deuterated molecular species are also detected with relatively strong intensities. On the other hand, emission lines of complex organic molecules such as HCOOCH3, and CH3OCH3 are found to be weak. For the molecules for which multiple transitions are detected, rotation temperatures are derived to be 7-33 K except for CH3OH. Emission lines with high upper-state energies (Eu > 150 K) are detected for CH3OH, indicating existence of a hot core. In comparison with the chemical composition of the Orion KL, carbon-chain molecules and deuterated molecules are found to be abundant in NGC 2264 CMM3, while sulfur-bearing species and complex organic molecules are deficient. These characteristics indicate chemical youth of NGC 2264 CMM3 in spite of its location at the center of the cluster forming core, NGC 2264 C.

To study physical properties of the natal filament gas around the cloud core harboring an exceptionally young low-mass protostar GF9-2, we carried out J=1-0 line observations of 12CO, 13CO, and C18O molecules using the Nobeyama 45m telescope. The mapping area covers ~1/5 of the whole filament. Our 13CO and C18O maps clearly demonstrate that the core formed at the local density maxima of the filament, and the internal motions of the filament gas are totally governed by turbulence with Mach number of ~2. We estimated the scale height of the filament to be H = 0.3 ~ 0.7 pc, yielding the central density of n_c = 700 ~4200 cm^-3. Our analysis adopting an isothermal cylinder model shows that the filament is supported by the turbulent and magnetic pressures against the radial and axial collapse due to self-gravity. Since both the dissipation time scales of the turbulence and the transverse magnetic fields can be comparable to the free-fall time of the filament gas of 10^6 years, we conclude that the local decay of the supersonic turbulence made the filament gas locally unstable, hence making the core collapse. Furthermore, we newly detected a gas condensation with velocity width enhancement to ~0.3 pc south-west of the GF9-2 core. The condensation has a radius of ~0.15 pc and an LTE mass of ~5 Msun. Its internal motion is turbulent with Mach number of ~3, suggestive of a gravitationally unbound state. Considering the uncertainties in our estimates, however, we propose that the condensation is a precursor of a cloud core which would have been produced by the collision of the two gas components identified in the filament.

We wish to investigate the structure, velocity field, and stellar content of the G16.59-0.05 high-mass star-forming region, where previous studies have established the presence of two almost perpendicular (NE-SW and SE-NW), massive outflows, and a rotating disk traced by methanol maser emission. We performed Very Large Array observations of the radio continuum and ammonia line emission, complemented by COMICS/Subaru and Hi-GAL/Herschel images in the mid- and far-infrared (IR). Our centimeter continuum maps reveal a collimated radio jet that is oriented E-W and centered on the methanol maser disk, placed at the SE border of a compact molecular core. The spectral index of the jet is negative, indicating non-thermal emission over most of the jet, except the peak close to the maser disk, where thermal free-free emission is observed. We find that the ammonia emission presents a bipolar structure consistent (on a smaller scale) in direction and velocity with that of the NE-SW bipolar outflow detected in previous CO observations. After analyzing our previous N2H+(1-0) observations again, we conclude that two scenarios are possible. In one case both the radio jet and the ammonia emission would trace the root of the large-scale CO bipolar outflow. The different orientation of the jet and the ammonia flow could be explained by precession and/or a non-isotropic density distribution around the star. In the other case, the N2H+(1-0) and ammonia bipolarity is interpreted as two overlapping clumps moving with different velocities along the line of sight. The ammonia gas also seems to undergo rotation consistent with the maser disk. Our IR images complemented by archival data allow us to derive a bolometric luminosity of about 10^4 L_sun and to conclude that most of the luminosity is due to the young stellar object associated with the maser disk.

We used the Very Long Baseline Array (VLBA) and the European VLBI Network (EVN) to perform phase-referenced VLBI observations of the three most powerful maser transitions associated with the high-mass star-forming region G28.87+0.07: the 22.2 GHz H$_{2}$O, 6.7 GHz CH$_{3}$OH, and 1.665 GHz OH lines. We also performed VLA observations of the radio continuum emission at 1.3 and 3.6 cm and Subaru observations of the continuum emission at 24.5 $\mu$m. Two centimeter continuum sources are detected and one of them (named "HMC") is compact and placed at the center of the observed distribution of H$_{2}$O, CH$_{3}$OH and OH masers. The bipolar distribution of line-of-sight (l.o.s) velocities and the pattern of the proper motions suggest that the water masers are driven by a (proto)stellar jet interacting with the dense circumstellar gas. The same jet could both excite the centimeter continuum source named "HMC" (interpreted as free-free emission from shocked gas) and power the molecular outflow observed at larger scales -- although one cannot exclude that the free-free continuum is rather originating from a hypercompact \ion{H}{2} region. At 24.5 $\mu$m, we identify two objects separated along the north-south direction, whose absolute positions agree with those of the two VLA continuum sources. We establish that $\sim$90% of the luminosity of the region ($\sim$\times10^{5} L_\sun$) is coming from the radio source "HMC", which confirms the existence of an embedded massive young stellar object (MYSO) exciting the masers and possibly still undergoing heavy accretion from the surrounding envelope.

Aims: The main goal of this study is to perform a sub-arcsecond resolution analysis of the high-mass star formation region G19.61-0.23, both in the continuum and molecular line emission. While the centimeter continuum images will be discussed in detail in a forthcoming paper, here we focus on the (sub)mm emission, devoting special attention to the hot molecular core. Results: Our observations resolve the HMC into three cores whose masses are on the order of 10^1-10^3 Msun. No submm core presents detectable free-free emission in the centimeter regime, but they appear to be associated with masers and thermal line emission from complex organic molecules. Towards the most massive core, SMA1, the CH3CN (18_K-17_K) lines reveal hints of rotation about the axis of a jet/outflow traced by H2O maser and H13CO+ (1--0) line emission. Inverse P-Cygni profiles of the 13CO (3--2) and C18O (3--2) lines seen towards SMA1 indicate that the central high-mass (proto)star(s) is (are) still gaining mass with an accretion rate $ge 3 ~10^{-3}$ Msun/yr. Due to the linear scales and the large values of the accretion rate, we hypothesize that we are observing an accretion flow towards a cluster in the making, rather than towards a single massive star.

The present paper focuses on the high-mass star-forming region G23.01-0.41. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G23.01-0.41: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the HMC at 1.3 and 3.6 cm. Results: We have detected H2O, CH3OH, and OH maser emission clustered within 2000 AU from the center of a flattened HMC, oriented SE-NW, from which emerges a massive 12CO outflow, elongated NE-SW, extended up to the pc-scale. Although the three maser species show a clearly different spatial and velocity distribution and sample distinct environments around the massive YSO, the spatial symmetry and velocity field of each maser specie can be explained in terms of expansion from a common center, which possibly denotes the position of the YSO driving the maser motion. Water masers trace both a fast shock (up to 50 km/s) closer to the YSO, powered by a wide-angle wind, and a slower (20 km/s) bipolar jet, at the base of the large-scale outflow. Since the compact free-free emission is found offset from the putative location of the YSO along a direction consistent with that of the maser jet axis, we interpret the radio continuum in terms of a thermal jet. The velocity field of methanol masers can be explained in terms of a composition of slow (4 km/s in amplitude) motions of radial expansion and rotation about an axis approximately parallel to the maser jet. Finally, the distribution of line of sight velocities of the hydroxyl masers suggests that they can trace gas less dense (n(H2) < 10^6 cm^-3) and more distant from the YSO than that traced by the water and methanol masers, which is expanding toward the observer. (Abridged)

Among the tracers of the earliest phases in the massive star formation process, methanol masers have gained increasing importance. The phenomenological distinction between Class I and II methanol masers is based on their spatial association with objects such as jets, cores, and ultracompact HII regions, but is also believed to correspond to different pumping mechanisms: radiation for Class II masers, collisions for Class I masers. In this work, we have surveyed a large sample of massive star forming regions - 296 objects divided into two groups named 'High' and 'Low' according to their [25-12] and [60-12] IRAS colours - in Class I and II methanol masers. Previous studies indicate that the High sources are likely more evolved. Therefore, the sample can be used to assess the existence of a sequence for the occurrence of Class I and II methanol masers during the evolution of a massive star forming region. We observed the 6 GHz (Class II) CH3OH maser with the Effelsberg 100-m telescope, and the 44 GHz and 95 GHz (Class I) CH3OH masers with the Nobeyama 45-m telescope. We have detected: 55 sources in the Class II line (12 new detections); 27 sources in the 44 GHz Class I line (17 new detections); 11 sources in the 95 GHz Class I line (all except one are new detections). Our statistical analysis shows that the ratio between the detection rates of Class II and Class I methanol masers is basically the same in High and Low sources. Therefore, both masers are equally associated with each evolutionary phase. In contrast, all maser species have about 3 times higher detection rates in High than in Low sources. This might indicate that the phenomena that originate all masers become progressively more active with time, during the earliest evolutionary phases of a high-mass star forming region.

The present paper focuses on the high-mass star-forming region G16.59-0.05. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G16.59-0.05: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the star-forming region at 1.3 and 3.6 cm. Results: This is the first work to report accurate measurements of the "relative" proper motions of the 6.7 GHz CH3OH masers. The different spatial and 3-D velocity distribution clearly indicate that the 22 GHz water and 6.7 GHz methanol masers are tracing different kinematic environments. The bipolar distribution of 6.7 GHz maser l.o.s. velocities and the regular pattern of observed proper motions suggest that these masers are tracing rotation around a central mass of about 35 solar masses. The flattened spatial distribution of the 6.7 GHz masers, oriented NW-SE, suggests that they can originate in a disk/toroid rotating around the massive YSO which drives the 12CO(2-1) outflow, oriented NE-SW, observed on arcsec scale. The extended, radio continuum source observed close to the 6.7 GHz masers could be excited by a wide-angle wind emitted from the YSO associated with the methanol masers, and such a wind is proven to be sufficiently energetic to drive the NE-SW 12CO(2-1) outflow. The H2O masers distribute across a region offset about 0".5 to the NW of the CH3OH masers, in the same area where emission of high-density molecular tracers, typical of HMCs, was detected. We postulate that a distinct YSO, possibly in an earlier evolutionary phase than that exciting the methanol masers, is responsible for the excitation of the water masers and the HMC molecular lines. (Abridged)

In order to investigate when and how the birth of a protostellar core occurs, we made survey observations of four well-studied dense cores in the Taurus molecular cloud using CO transitions in submillimeter bands. We report here the detection of unexpectedly warm (~ 30 - 70 K), extended (radius of ~ 2400 AU), dense (a few times 10^{5} cm^{-3}) gas at the heart of one of the dense cores, L1521F (MC27), within the cold dynamically collapsing components. We argue that the detected warm, extended, dense gas may originate from shock regions caused by collisions between the dynamically collapsing components and outflowing/rotating components within the dense core. We propose a new stage of star formation, "warm-in-cold core stage (WICCS)", i.e., the cold collapsing envelope encases the warm extended dense gas at the center due to the formation of a protostellar core. WICCS would constitutes a missing link in evolution between a cold quiescent starless core and a young protostar in class 0 stage that has a large-scale bipolar outflow.

Using the Caltech Submillimeter Observatory 10.4-meter telescope, we performed sensitive mapping observations of 12CO J=7-6 emission at 807 GHz towards Orion IRc2. The image has an angular resolution of 10", which is the highest angular resolution data toward the Orion Hot Core published for this transition. In addition, thanks to the on-the-fly mapping technique, the fidelity of the new image is rather high, particularly in comparison to previous images. We have succeeded in mapping the northwest-southeast high-velocity molecular outflow, whose terminal velocity is shifted by ~70-85 km/s with respect to the systemic velocity of the cloud. This yields an extremely short dynamical time scale of ~900 years. The estimated outflow mass loss rate shows an extraordinarily high value, on the order of 10^{-3} Msun/yr. Assuming that the outflow is driven by Orion IRc2, our result agrees with the picture so far obtained for a 20 Msun (proto)star in the process of formation.

We present spectroscopic evidence for infall motion of gas in the natal cloud core harboring an extremely young low-mass protostar GF9-2. We previously discussed that the ongoing collapse of the GF9-2 core has agreement with the Larson-Penston-Hunter (LPH) theoretical solution for the gravitational collapse of a core (Furuya et al.; paper I). To discuss the gas infall on firmer ground, we have carried out On-The-Fly mapping observations of the HCO+ (1--0) line using the Nobeyama 45m telescope equipped with the 25 Beam Array Receiver System. Furthermore, we observed the HCN (1--0) line with the 45m telescope, and the HCO+ (3--2) line with the Caltech Submillimeter Observatory 10.4 m telescope. The optically thick HCO+ and HCN lines show blueskewed profiles whose deepest absorptions are seen at the peak velocity of optically thin lines, i.e., the systemic velocity of the cloud (paper I), indicating the presence of gas infall toward the central protostar. We compared the observed HCO+ line profiles with model ones by solving the radiative transfer in the core under LTE assumption.We found that the core gas has a constant infall velocity of ~0.5 km/s in the central region, leading to a mass accretion rate of 2.5x10^{-5} Msun/yr. Consequently, we confirm that the gas infall in the GF9-2 core is consistent with the LPH solution.

We carried out an unbiased mapping survey of dense molecular cloud cores traced by the NH3 (1,1) and (2,2) inversion lines in the GF9 filament which contains an extremely young low-mass protostar GF9-2 (Furuya et al. 2006, ApJ, 653, 1369). The survey was conducted using the Nobeyama 45m telescope over a region of ~1.5 deg with an angular resolution of 73". The large-scale map revealed that the filament contains at least 7 dense cores, as well as 3 possible ones, located at regular intervals of ~0.9 pc. Our analysis shows that these cores have kinetic temperatures of $\lesssim$ 10 K and LTE-masses of 1.8 -- 8.2 Msun, making them typical sites of low-mass star formation. All the identified cores are likely to be gravitationally unstable because their LTE-masses are larger than their virial masses. Since the LTE-masses and separations of the cores are consistent with the Jeans masses and lengths, respectively, for the low-density ambient gas, we argue that the identified cores have formed via the gravitational fragmentation of the natal filamentary cloud.

We used the CSO 10.4 meter telescope to image the 350 micron and 450 micron continuum and CO J=6-5 line emission of the IRAS 20126+4104 clump. The continuum and line observations show that the clump is isolated over a 4 pc region and has a radius of ~ 0.5 pc. Our analysis shows that the clump has a radial density profile propto r ^{-1.2} for r <~ 0.1 pc and has propto r^{-2.3} for r >~ 0.1 pc which suggests the inner region is infalling, while the infall wave has not yet reached the outer region. Assuming temperature gradient of r^{-0.35}, the power law indices become propto r ^{-0.9} for r < ~0.1 pc and propto r^{-2.0} for r >~ 0.1 pc. Based on a map of the flux ratio of 350micron/450micron, we identify three distinct regions: a bipolar feature that coincides with the large scale CO bipolar outflow; a cocoon-like region that encases the bipolar feature and has a warm surface; and a cold layer outside of the cocoon region. The complex patterns of the flux ratio map indicates that the clump is no longer uniform in terms of temperature as well as dust properties. The CO emission near the systemic velocity traces the dense clump and the outer layer of the clump shows narrow line widths (< ~3 km/s). The clump has a velocity gradient of ~ 2 km/s pc^{-1}, which we interpret as due to rotation of the clump, as the equilibrium mass (~ 200 Msun) is comparable to the LTE mass obtained from the CO line. Over a scale of ~ 1 pc, the clump rotates in the opposite sense with respect to the >~ 0.03 pc disk associated with the (proto)star. This is one of four objects in high-mass and low-mass star forming regions for which a discrepancy between the rotation sense of the envelope and the core has been found, suggesting that such a complex kinematics may not be unusual in star forming regions.

Using the OVRO, Nobeyama, and IRAM mm-arrays, we searched for ``disk''-outflow systems in three high-mass (proto)star forming regions: G16.59-0.05, G23.01-0.41, and G28.87+0.07. These were selected from a sample of NH3 cores associated with OH and H2O maser emission and with no or very faint continuum emission. Our imaging of molecular line (including rotational transitions of CH3CN and 3mm dust continuum emission revealed that these are compact, massive, and hot molecular cores (HMCs), that is likely sites of high-mass star formation prior to the appearance of UCHII regions. All three sources turn out to be associated with molecular outflows from CO and/or HCO+ J=1--0 line imaging. In addition, velocity gradients of 10 -- 100 km/s per pc in the innermost densest regions of the G23.01 and G28.87 HMCs are identified along directions roughly perpendicular to the axes of the corresponding outflows. All the results suggest that these cores might be rotating about the outflow axis, although the contribution of rotation to gravitational equilibrium of the HMCs appears to be negligible. Our analysis indicates that the 3 HMCs are close to virial equilibrium due to turbulent pressure support. Comparison with other similar objects where rotating toroids have been identified so far shows that in our case rotation appears to be much less prominent; this can be explained by the combined effect of unfavorable projection, large distance, and limited angular resolution with the current interferometers.

Theory predicts and observations confirm that low-mass stars (like the Sun) in their early life grow by accreting gas from the surrounding material. But for stars ~ 10 times more massive than the Sun (~10 M_sun), the powerful stellar radiation is expected to inhibit accretion and thus limit the growth of their mass. Clearly, stars with masses >10 M_sun exist, so there must be a way for them to form. The problem may be solved by non-spherical accretion, which allows some of the stellar photons to escape along the symmetry axis where the density is lower. The recent detection of rotating disks and toroids around very young massive stars has lent support to the idea that high-mass (> 8 M_sun) stars could form in this way. Here we report observations of an ammonia line towards a high-mass star forming region. We conclude from the data that the gas is falling inwards towards a very young star of ~20 M_sun, in line with theoretical predictions of non-spherical accretion.

We present a study of the natal core harboring the class 0 protostar GF9-2 in the filamentary dark cloud GF 9 (d = 200 pc). GF9-2 stands unique in the sense that it shows H2O maser emission, a clear signpost of protostar formation, whereas it does not have a high-velocity large-scale molecular outflow evidenced by our deep search for CO wing emission. These facts indicate that GF9-2 core is early enough after star formation so that it still retains some information of initial conditions for collapse. Our 350 um dust continuum emission image revealed the presence of a protostellar envelope in the center of a molecular core. The mass of the envelope is ~0.6 Msun from the 350 um flux density, while LTE mass of the core is ~3 Msun from moleuclar line observations. Combining visibility data from the OVRO mm-array and the 45m telescope, we found that the core has a radial density profile of $\rho(r)\propto r^{-2}$ for 0.003 < r/pc < 0.08 region. Molecular line data analysis revealed that the velocity width of the core gas increases inward,while the outermost region maintains a velocity dispersion of a few times of the ambient sound speed. The broadened velocity width can be interpreted as infall. Thus, the collapse in GF9-2 is likely to be described by an extension of the Larson-Penston solution for the period after formation of a central star. We derived the current mass accretion rate of ~3E-05 Msun/year from infall velocity of ~ 0.3 km/s at r~ 7000 AU. All results suggest that GF9-2 core has been undergoing gravitational collapse for ~ 5000 years since the formation of central protostar(s), and that the unstable state initiated the collapse ~2E+05 years (the free-fall time) ago.

This article reports first results of a long-term observational program aimed to study the earliest evolution of jet/disk systems in low-mass YSOs by means of VLBI observations of the 22.2 GHz water masers. We report here data for the cluster of low-mass YSOs in the Serpens molecular core and for the single object RNO~15-FIR. Towards Serpens SMM1, the most luminous sub-mm source of the Serpens cluster, the water maser emission comes from two small (< 5 AU in size) clusters of features separated by ~25 AU, having line of sight velocities strongly red-shifted (by more than 10 km/s) with respect to the LSR velocity of the molecular cloud. The two maser clusters are oriented on the sky along a direction that is approximately perpendicular to the axis of the radio continuum jet observed with the VLA towards SMM1. The spatial and velocity distribution of the maser features lead us to favor the interpretation that the maser emission is excited by interaction of the receding lobe of the jet with dense gas in the accretion disk surrounding the YSO in SMM1. Towards RNO~15-FIR, the few detected maser features have both positions and (absolute) velocities aligned along a direction that is parallel to the axis of the molecular outflow observed on much larger angular scales. In this case the maser emission likely emerges from dense, shocked molecular clumps displaced along the axis of the jet emerging from the YSO. The protostar in Serpens SMM1 is more massive than the one in RNO~15-FIR. We discuss the case where a high mass ejection rate can generate jets sufficiently powerful to sweep away from their course the densest portions of circumstellar gas. In this case, the excitation conditions for water masers might preferably occur at the interface between the jet and the accretion disk, rather than along the jet axis.

We conducted a 4 epoch 3 month VLBA proper motion study of H$_2$O masers toward an intermediate-mass class 0 source IRAS 20050+2720 MMS1 (d=700 pc). From milli-arcsecond (mas) resolution VLBA images, we found two groups of H2O maser spots at the center of the submillimeter core of MMS1. One group consists of more than $\sim 50$ intense maser spots; the other group consisting of several weaker maser spots is located at 18 AU south-west of the intense group. Distribution of the maser spots in the intense group shows an arc-shaped structure which includes the maser spots that showed a clear velocity gradient. The spatial and velocity structures of the maser spots in the arc-shape did not significantly change through the 4 epochs. Furthermore, we found a relative proper motion between the two groups. Their projected separation increased by 1.13+/-0.11 mas over the 4 epochs along a line connecting them. The spatial and velocity structures of the intense group and the relative proper motions strongly suggest that the maser emission is associated with a protostellar jet. Comparing the observed LSR velocities with calculated radial velocities from a simple biconical jet model, we conclude that the most of the maser emission are likely to be associated with an accelerating biconical jet which has large opening angle. The large opening angle of the jet traced by the masers would support the hypothesis that poor jet collimation is an inherent property of luminous (proto)stars.

We present the results of high angular resolution millimeter observations of gas and dust toward G31.41+0.31 and G24.78+0.08, two high-mass star forming regions where four rotating massive toroids have been previously detected by Beltran et al. (2004). The CH3CN (12-11) emission of the toroids in G31.41+0.31 and core A1 in G24.78+0.08 has been modeled assuming that it arises from a disk-like structure seen edge-on, with a radial velocity field. For G31.41+0.31 the model properly fits the data for a velocity v_rot~1.7 km/s at the outer radius R_out~13400 AU and an inner radius R_inn~1340 AU, while for core A1 in G24.78+0.08 the best fit is obtained for v_rot~2.0 km/s at R_out~7700 AU and R_inn~2300 AU. Unlike the rotating disks detected around less luminous stars, these toroids are not undergoing Keplerian rotation. From the modeling itself, however, it is not possible to distinguish between constant rotation or constant angular velocity, since both velocity fields suitably fit the data. The best fit models have been computed adopting a temperature gradient of the type T proportional R^{-3/4}, with a temperature at the outer radius T_out~100 K for both cores. The M_dyn needed for equilibrium derived from the models is much smaller than the mass of the cores, suggesting that such toroids are unstable and undergoing gravitational collapse. The collapse is also supported by the CH3^{13}CN or CH3CN line width measured in the cores, which increases toward the center of the toroids. The estimates of v_inf and \dot M_acc are ~2 km/s and 3x10^{-2} M_sun/yr for G31.41+0.31, and ~1.2 km/s and ~9x10^{-3} M_sun/yr for G24.78+0.08 A1. Such large accretion rates could weaken the effect of stellar winds and radiation pressure and allow further accretion on the star.

Using the Owens Valley and Nobeyama Radio Observatory interferometers, we carried out an unbiased search for hot molecular cores and ultracompact UC HII regions toward the high-mass star forming region G19.61--0.23. In addition, we performed 1.2 mm imaging with SIMBA, and retrieved 3.5 and 2 cm images from the VLA archive data base. The newly obtained 3 mm image brings information on a cluster of high-mass (proto)stars located in the innermost and densest part of the parsec scale clump detected in the 1.2 mm continuum. We identify a total of 10 high-mass young stellar objects: one hot core (HC) and 9 UC HII regions, whose physical parameters are obtained from model fits to their continuum spectra. The ratio between the current and expected final radii of the UC \HII regions ranges from 0.3 to 0.9, which leaves the possibility that all O-B stars formed simultaneously. Under the opposite assumption -- namely that star formation occurred randomly -- we estimate that HC lifetime is less than $\sim$1/3 of that of UCHII regions on the basis of the source number ratio between them.

We report on the detection of four rotating massive disks in two regions of high-mass star formation. The disks are perpendicular to known bipolar outflows and turn out to be unstable but long lived. We infer that accretion onto the embedded (proto)stars must proceed through the disks with rates of ~10E-2 Msun/yr.

We have conducted multi-epoch synthesis imaging of 2 and 3 millimeter (mm) continuum emission and near infrared K band 2.2 um imaging of a flare event in January 2003 that occurred on the young stellar object GMR-A which is suggested to be a weak-line T Tauri star in the Orion cluster. Our mm data showed that the flare activity lasted at least over 13 days, whereas the K-band magnitude did not change during this event. In addition, we have succeeded in detecting short time variations of flux on the time scales of 15 minutes. The total energy of the flare is estimated to be 10^{35-36} erg, which makes it one of the most energetic flares reported to date. Comparing the mm continuum luminosities with reported X-ray luminosities, we conclude that the mm flare was similar in nature to solar and other stellar flares. Our results will be a crucial step toward understanding magnetically induced stellar surface activities in T Tauri stars.

We present results from a study of molecular gas and dust in the vicinity of the Herbig Haro object HH2. Emission from the sub-mm continuum, 12CO and HCO+ was mapped with angular resolutions ranging from 14 arcsec to 5 arcsec (or 0.01pc at the distance of HH2). The continuum shows an extended dust clump of mass 3.8Msun and temperature 22K, located downstream of the bright optical HH knots. However, a compact emission peak lies within ~0.01pc of the low-excitation H2-prominent shocks, with a luminosity consistent with local heating by the outflow. The HCO+ emission shows two velocity components: firstly, ambient-velocity gas lying in a region roughly corresponding to the dust clump, with abundance enhanced by a factor of a few close to the H2-prominent knots. Secondly a component of high-velocity emission (20 km/s linewidth), found mainly in a collimated jet linking the low-excitation HH objects. In this high-velocity jet, the line wings show an abundance ratio HCO+/CO proportional to v^2, with an HCO+ enhancement compared with ambient gas of up to 10^3 at the most extreme velocities. Such high abundances are consistent with models of shock chemistry in turbulent mixing layers at the interaction boundaries of jets. Extrapolating this effect to low velocities suggests that the more modest HCO+ enhancement in the clump gas could be caused by low velocity shocks. A UV precursor may not, therefore be necessary to explain the elevated HCO+ abundance in this gas.

We conducted VLA observations at 0".06 resolution of the 22 GHz water masers toward the Class 0 source S106 FIR (d=600 pc; 15" west of S106-IRS4) on two epochs separated by ~3 months. Two compact clusters of the maser spots were found in the center of the submillimeter core of S106 FIR. The separation of the clusters was ~80 mas (48 AU) along P. A. = 70 degrees and the size of each cluster was ~20 mas x 10 mas. The western cluster, which had three maser components, was 7.0 km/s redshifted with respect to the ambient cloud velocity. Each component was composed of a few spatially localized maser spots and was aligned on a line connecting the clusters. We found relative proper motions of the components with ~30 mas/yr (18 AU/yr) along the line. In addition, a series of single-dish observations show that the maser components drifted with a radial acceleration of ~1 km/s/yr. These facts indicate that the masers could be excited by a 10 AU-scale jet-like accelerating flow ejected from an assumed protostar located between the two clusters. The outflow size traced by the masers was 50 AU x 5 AU after correction for an inclination angle of 10 degrees which was derived from the relative proper motions and radial velocities of the maser components. The three-dimensional outflow velocity ranged from 40 to 70 km/s assuming symmetric motions for the blue and red components. Since no distinct CO molecular outflows have been detected so far, we suggest that S106 FIR is an extremely young protostar observed just after the onset of outflowing activity.