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.

The mass growth of protostars is a central element to the determination of fundamental stellar population properties such as the initial mass function. Constraining the accretion history of individual protostars is therefore an important aspect of star formation research. The goal of the study presented here is to determine whether high-mass (proto)stars gain their mass from a compact (<0.1pc) fixed-mass reservoir of gas, often referred to as dense cores, in which they are embedded, or whether the mass growth of high-mass stars is governed by the dynamical evolution of the parsec-scale clump that typically surrounds them. To achieve this goal, we performed a 350micron continuum mapping of 11 infrared dark clouds, along side some of their neighbouring clumps, with the ArT\'eMiS camera on APEX. By identifying about 200 compact ArT\'eMiS sources, and matching them with Herschel Hi-GAL 70micron sources, we have been able to produce mass vs. temperature diagrams. We compare the nature (i.e. starless or protostellar) and location of the ArT\'eMiS sources in these diagrams with modelled evolutionary tracks of both core-fed and clump-fed accretion scenarios. We argue that the latter provide a better agreement with the observed distribution of high-mass star-forming cores. However, a robust and definitive conclusion on the question of the accretion history of high-mass stars requires larger number statistics.

Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArT\'eMiS camera on the APEX telescope and combined the ArT\'eMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArT\'eMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

We identify 45 protostellar outflows in CO maps of the Orion A giant molecular cloud from the CARMA-NRO Orion survey. Our sample includes 11 newly detected outflows. We measure the mass and energetics of the outflows, including material at low-velocities by correcting for cloud contributions. The total momentum and kinetic energy injection rates of outflows is comparable to the turbulent dissipation rate of the cloud. We also compare the outflow position angles to the orientation of C$^{18}$O filaments. We find that the full sample of outflows is consistent with being randomly oriented with respect to the filaments. A subsample of the most reliable measurements shows a moderately perpendicular outflow-filament alignment which may reflect accretion of mass across filaments and onto the protostellar cores.

Herschel observations of nearby clouds in the Gould Belt support a paradigm for low-mass star formation, starting with the generation of molecular filaments, followed by filament fragmentation, and the concentration of mass into self-gravitating prestellar cores. With the unique far-infrared and submillimeter continuum imaging capabilities of the Herschel Space observatory, the closeby (d = 139 pc) Ophiuchus cloud was mapped at five wavelengths from 70 ${\mu}m$ to 500 ${\mu}m$ with the aim of providing a complete census of dense cores in this region, including unbound starless cores, bound prestellar cores, and protostellar cores. Taking advantage of the high dynamic range and multi-wavelength nature of the Herschel data, we used the multi-scale decomposition algorithms getsources and getfilaments to identify a complete sample of dense cores and filaments in the cloud and study their properties. The densest clouds of the Ophiuchus complex, L1688 and L1689, which thus far are only indirectly described as filamentary regions owing to the spatial distribution of their young stellar objects (YSOs), are confirmed to be dominated by filamentary structures. The tight correlation observed between prestellar cores and filamentary structures in L1688 and L1689 supports the view that solar-type star formation occurs primarily in dense filaments. While the sub clouds of the complex show disparities, L1689 being less efficient than L1688 at forming stars when considering their total mass budgets, both sub clouds share almost the same prestellar core formation efficiency in dense molecular gas. We also find evidence in the Herschel data for a remarkable concentric geometrical configuration in L1688 which is dominated by up to three arc-like compression fronts and presumably created by shockwave events emanating from the Sco OB2 association, including the neighboring massive (O9V) star ${\sigma}$ Sco.

We present a detailed study of the Orion B clouds (d~400 pc), imaged with the PACS/SPIRE cameras at 70-500 $\mu$m by the Herschel Gould Belt survey (HGBS). We release new high-res. maps of column density and dust temperature. In the filamentary sub-regions NGC2023/2024, NGC2068/2071, and L1622, 1768 starless dense cores were identified, ~28-45% of which are self-gravitating prestellar cores. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was found to be $t_{\rm pre}=1.7_{-0.6}^{+0.8}$ Myr. The prestellar core mass function (CMF) peaks at ~0.5 $M_\odot$ and is consistent with a power law with log slope -1.27$\pm$0.24 at the high-mass end, compared to the Salpeter slope of -1.35. In this region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A_V_bg~7 mag in background visual extinction, similar to the trend observed with Herschel in other clouds. This is not a sharp threshold, but a smooth transition between a regime with very low prestellar CFE at A_V_bg<5 and a regime with higher, roughly constant CFE at A_V_bg$\gtrsim$10. The total mass in the form of prestellar cores represents only ~20% of the dense molecular cloud gas at A_V_bg$\gtrsim$7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation between nearest core neighbors corresponds to the typical inner filament width of ~0.1 pc commonly observed in nearby molecular clouds. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser filaments.

The NIKA2 polarization channel at 260 GHz (1.15 mm) has been proposed primarily to observe galactic star-forming regions and probe the critical scales between 0.01-0.05 pc at which magnetic field lines may channel the matter of interstellar filaments into growing dense cores. The NIKA2 polarimeter consists of a room temperature continuously rotating multi-mesh HWP and a cold polarizer that separates the two orthogonal polarizations onto two 260 GHz KIDs arrays. We describe in this paper the preliminary results obtained during the most recent commissioning campaign performed in December 2018. We concentrate here on the analysis of the extended sources, while the observation of compact sources is presented in a companion paper [12]. We present preliminary NIKA2 polarization maps of the Crab nebula. We find that the integrated polarization intensity flux measured by NIKA2 is consistent with expectations.In terms of polarization angle, we are still limited by systematic uncertainties that will be further investigated in the forthcoming commissioning campaigns.

Herschel surveys of Galactic clouds support a paradigm for low-mass star formation in which dense filaments play a crucial role. The detailed fragmentation properties of star-forming filaments remain poorly understood, however, and the validity of the filament paradigm in the high-mass regime is still unclear. To investigate the density/velocity structure of the filament in the high-mass star-forming region NGC6334, we conducted ALMA observations in the 3mm continuum and the N2H+(1-0) line at ~3arcsec resolution. The filament was detected in both tracers. We identified 26 cores at 3mm and 5 velocity-coherent fiber-like features in N2H+ within the filament. The typical length of, and velocity difference between, the fiber-like features of the NGC6334 filament are reminiscent of the properties for the fibers of the low-mass star-forming filament B211/B213. Only 2 or 3 of the 5 velocity-coherent features are well aligned with the filament and may represent genuine, fiber sub-structures. The core mass distribution has a peak at ~10Msun. They can be divided into 7 groups of cores, closely associated with ArTeMiS clumps. The projected separation between cores and the projected spacing between clumps are roughly consistent with the effective Jeans length in the filament and a physical scale of about 4 times the filament width, respectively, suggesting a bimodal filament fragmentation process. Despite being one order of magnitude denser and more massive than the B211/B213 filament, the NGC6334 filament has a similar density/velocity structure. The difference is that the cores in NGC6334 appear to be an order of magnitude denser and more massive than the cores in Taurus. This suggests that dense filaments may evolve and fragment in a similar manner in low- and high-mass star-forming regions, and that the filament paradigm may hold in the intermediate-mass (if not high-mass) star formation regime.

We conducted an exploration of 12CO molecular outflows in the Orion A giant molecular cloud to investigate outflow feedback using 12CO (J = 1-0) and 13CO (J = 1-0) data obtained by the Nobeyama 45-m telescope. In the region excluding the center of OMC 1, we identified 44 12CO (including 17 newly detected) outflows based on the unbiased and systematic procedure of automatically determining the velocity range of the outflows and separating the cloud and outflow components. The optical depth of the 12CO emission in the detected outflows is estimated to be approximately 5. The total momentum and energy of the outflows, corrected for optical depth, are estimated to be 1.6 x 10 2 M km s-1 and 1.5 x 10 46 erg, respectively. The momentum and energy ejection rate of the outflows are estimated to be 36% and 235% of the momentum and energy dissipation rates of the cloud turbulence, respectively. Furthermore, the ejection rates of the outflows are comparable to those of the expanding molecular shells estimated by Feddersen et al. (2018, ApJ, 862, 121). Cloud turbulence cannot be sustained by the outflows and shells unless the energy conversion efficiency is as high as 20%.

We carried out mapping observations toward three nearby molecular clouds, Orion A, Aquila Rift, and M17, using a new 100 GHz receiver, FOREST, on the Nobeyama 45-m telescope. In the present paper, we describe the details of the data obtained such as intensity calibration, data sensitivity, angular resolution, and velocity resolution. Each target contains at least one high-mass star-forming region. The target molecular lines were $^{12}$CO ($J = 1 - 0$), $^{13}$CO ($J = 1 - 0$), C$^{18}$O ($J = 1 - 0$), N$_2$H$^+$ ($J=1-0$), and CCS ($J_N=8_7-7_6$), with which we covered the density range of 10$^2$ cm$^{-3}$ to 10$^6$ cm$^{-3}$ with an angular resolution of $\sim 20\arcsec$ and a velocity resolution of $\sim$ 0.1 km s$^{-1}$. Assuming the representative distances of 414 pc, 436 pc, and 2.1 kpc, the maps of Orion A, Aquila Rift, and M17 cover most of the densest parts with areas of about 7 pc $\times$ 15 pc, 7 pc $\times$ 7 pc, and 36 pc $\times$ 18 pc, respectively. On the basis of the $^{13}$CO column density distribution, the total molecular masses are derived to be $3.86 \times 10^4 M_\odot$, $2.67 \times 10^4 M_\odot$, and $8.1\times 10^5 M_\odot$ for Orion A, Aquila Rift, and M17, respectively. For all the clouds, the H$_2$ column density exceeds the theoretical threshold for high-mass star formation of $\gtrsim$ 1 g cm$^{-2}$, only toward the regions which contain current high-mass star-forming sites. For other areas, further mass accretion or dynamical compression would be necessary for future high-mass star formation. This is consistent with the current star formation activity. Using the $^{12}$CO data, we demonstrate that our data have enough capability to identify molecular outflows, and for Aquila Rift, we identify 4 new outflow candidates. The scientific results will be discussed in details in separate papers.

We present the results of mapping observations toward an outflow-shocked region, OMC-2 FIR 4 using the Nobeyama 45-m telescope. We observed the area in $^{13}$CO ($J=1-0$), C$^{18}$O ($J=1-0$), N$_2$H$^+$ ($J=1-0$), CCS ($J_N=8_7-7_6$), HCO$^+$ ($J=1-0$), H$^{13}$CO$^+$ ($J=1-0$), HN$^{13}$C ($J=1-0$), H$^{13}$CN ($J=1-0$), DNC ($J=1-0$), N$_2$D$^+$ ($J=1-0$), and DC$_3$N ($J=9-8$). We detected a dense molecular clump that contains FIR 4/5. We also detected in $^{13}$CO blueshifted and redshifted components driven presumably by protostellar outflows in this region. The axes of the FIR 3 and VLA 13 outflows, projected on the plane of the sky, appear to point toward the FIR 4 clump, suggesting that the clump may be compressed by protostellar outflows from Class I sources, FIR 3 and VLA 13. Applying the hyperfine fit of N$_2$H$^+$ lines, we estimated the excitation temperature to be $\sim$ 20 K. The high excitation temperature is consistent with the fact that the clump contains protostars. The CCS emission was detected in this region for the first time. Its abundance is estimated to be a few $\times 10^{-12}$, indicating that the region is chemically evolved at $\sim 10^5$ years, which is comparable to the typical lifetime of the Class I protostars. This timescale is consistent with the scenario that star formation in FIR 4 is triggered by dynamical compression of the protostellar outflows. The [HNC]/[HCN] ratio is evaluated to be $\sim 0.5$ in the dense clump and the outflow lobes, whereas it is somewhat larger in the envelope of the dense clump. The small [HNC]/[HCN] ratio indicates that the HNC formation was prevented due to high temperatures. Such high temperatures seem to be consistent with the scenario that either protostellar radiation or outflow compression, or both, affected the thermal properties of this region.

We present results of the classification of cloud structures toward the Orion A Giant Molecular Cloud based on wide-field $^{12}$CO ($J$ = 1--0), $^{13}$CO ($J$ = 1-0), and C$^{18}$O ($J$ = 1-0) observations using the Nobeyama 45 m radio telescope. We identified 78 clouds toward Orion A by applying Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) to the data cube of the column density of $^{13}$CO. Well-known subregions such as OMC-1, OMC-2/3, OMC-4, OMC-5, NGC 1977, L1641-N, and the dark lane south filament (DLSF) are naturally identified as distinct structures in Orion A. These clouds can also be classified into three groups: the integral-shaped filament, the southern regions of Orion A, and the other filamentary structures in the outer parts of Orion A and the DLSF. These groups show differences in scaling relations between the physical properties of the clouds. We derived the abundance ratio between $^{13}$CO and C$^{18}$O, $X_{^{13}\mathrm{CO}}/X_{\mathrm{C}^{18}\mathrm{O}}$, which ranges from 5.6 to 17.4 on median over the individual clouds. The significant variation of $X_{^{13}\mathrm{CO}}/X_{\mathrm{C}^{18}\mathrm{O}}$ is also seen within a cloud in both of the spatial and velocity directions and the ratio tends to be high at the edge of the cloud. The values of $X_{^{13}\mathrm{CO}}/X_{\mathrm{C}^{18}\mathrm{O}}$ decrease from 17 to 10 with the median of the column densities of the clouds at the column density of $N_{\mathrm{C^{18}O}} \gtrsim 1 \times 10^{15}$ cm$^{-2}$ or visual extinction of $A_V \gtrsim 3$ mag under the strong far-ultraviolet (FUV) environment of $G_0 > 10^{3}$, whereas it is almost independent of the column density in the weak FUV radiation field. These results are explained if the selective photodissociation of C$^{18}$O is enhanced under a strong FUV environment and it is suppressed in the dense part of the clouds.

The origin of the stellar initial mass function (IMF) is one of the most debated issues in astrophysics. Here, we explore the possible link between the quasi-universal filamentary structure of star-forming molecular clouds and the origin of the IMF. Based on our recent comprehensive study of filament properties from Herschel Gould Belt survey observations (Arzoumanian et al.), we derive, for the first time, a good estimate of the filament mass function (FMF) and filament line mass function (FLMF) in nearby molecular clouds. We use the observed FLMF to propose a simple toy model for the origin of the prestellar core mass function (CMF), relying on gravitational fragmentation of thermally supercritical but virialized filaments. We find that the FMF and the FLMF have very similar shapes and are both consistent with a Salpeter-like power-law function (d$N$/dlog$M_{\rm line} \propto M_{\rm line}^{-1.5\pm0.1}$) in the regime of thermally supercritical filaments ($M_{\rm line} > 16\, M_\odot$/pc). This is a remarkable result since, in contrast, the mass distribution of molecular clouds and clumps is known to be significantly shallower than the Salpeter power-law IMF, with d$N$/dlog$M_{\rm cl} \propto M_{\rm cl}^{-0.7}$. Since the vast majority of prestellar cores appear to form in thermally transcritical or supercritical filaments, we suggest that the prestellar CMF and by extension the stellar IMF are at least partly inherited from the FLMF through gravitational fragmentation of individual filaments.

A survey of molecular cores covering the infrared dark cloud known as the M17 southwest extension (M17 SWex) has been carried out with the 45 m Nobeyama Radio Telescope. Based on the N2H+ (J= 1-0) data obtained, we have identified 46 individual cores whose masses are in the range 43 to 3026 Mo. We examined the relationship between the physical parameters of the cores and those of young stellar objects (YSOs) associated with the cores found in the literature. The comparison of the virial mass and the core mass indicates that most of the cores can be gravitationally stable if we assume a large external pressure. Among the 46 cores, we found four massive cores with YSOs. They have large mass of >~ 1000Mo and line width of >~ 2.5 km s^-1 which are similar to those of clumps forming high mass stars. However, previous studies have shown that there is no active massive star formation in this region. Recent measurements of near-infrared polarization infer that the magnetic field around M17 SWex is likely to be strong enough to support the cores against self-gravity. We therefore suggest that the magnetic field may prevent the cores from collapsing, causing the low-level of massive star formation in M17 SWex.

We report on a possible cloud-cloud collision in the DR 21 region, which we found through molecular observations with the Nobeyama 45-m telescope. We mapped an area of 8'x12' around the region with twenty molecular lines including the 12CO(J=1-0) and 13CO(J=1-0) emission lines, and sixteen of them were significantly detected. Based on the 12CO and 13CO data, we found five distinct velocity components in the observed region, and we call molecular gas associated with these components -42, -22, -3, 9, and 17 km/s clouds taking after their typical radial velocities. The -3 km/s cloud is the main filamentary cloud (31,000 Mo) associated with young massive stars such as DR21 and DR21(OH), and the 9 km/s cloud is a smaller cloud (3,400 Mo) which may be an extension of the W75 region in the north. The other clouds are much smaller. We found a clear anticorrelation in the distributions of the -3 and 9 km/s clouds, and detected faint 12CO emission having intermediate velocities bridging the two clouds at their intersection. These facts strongly indicate that the two clouds are colliding against each other. In addition, we found that DR21 and DR21(OH) are located in the periphery of the densest part of the 9 km/s cloud, which is consistent with results of recent numerical simulations of cloud-cloud collisions. We therefore suggest that the -3 and 9 km/s clouds are colliding, and that the collision induced the massive star formation in the DR21 cloud. The interaction of the -3 and 9 km/s clouds was previously suggested by Dickel et al. (1978), and our results strongly support their hypothesis of the interaction.

We report possible interaction between multiple super nova remnants (SNRs) and Northern Coal Sack (NCS) which is a massive clump (~1000 Mo) in the Cyg OB 7 cloud complex and is forming a massive Class 0 object. We performed molecular observations of the 12CO(J=1-0), 13CO(J=1-0), and C18O(J=1-0) emission lines using the 45m telescope at the Nobeyama Radio Observatory, and we found that there are mainly four velocity components at Vlsr=-20, -6, -4, and 10 km/s. The -6 and -4 km/s components correspond to the systemic velocities of NCS and the Cygnus OB 7 complex, respectively, and the other velocity components originate from distinct smaller clouds. Interestingly, there are apparent correlations and anti-correlations among the spatial distributions of the four components, suggesting that they are physically interacting with one another. On a larger scale, we find that a group of small clouds belonging to the -20 and 10 km/s components are located along two different arcs around some SNRs including HB 21 which has been suggested to be interacting with the Cyg OB 7 cloud complex, and we also find that NCS is located right at the interface of the arcs. The small clouds are likely to be the gas swept up by the stellar wind of the massive stars which created the SNRs. We suggest that the small clouds alined along the two arcs recently encountered NCS and the massive star formation in NCS was triggered by the strong interaction with the small clouds.

We analyze $^{12}$CO($J$=2-1), $^{13}$CO($J$=2-1), C$^{18}$O ($J$=2-1), and 1.3 mm continuum maps of the $\rho$ Ophiuchi A photo-dissociation region (PDR) obtained with ALMA. Layered structures of the three CO isotopologues with an angular separation of 10 arcsec = 6.6$\times$10$^{-3}$ pc = 1400 au are clearly detected around the Be star, S1 (i.e., each front of emission shifts from the near to far side in order of $^{12}$CO, $^{13}$CO, and C$^{18}$O). We estimate the spatial variations of $X$($^{13}$CO)/$X$(C$^{18}$O) abundance ratios, and find that the abundance ratio is as high as 40 near the emission front, and decreases to the typical value in the solar system of 5.5 in a small angular scale of 4 arcsec = 2.6$\times$10$^{-3}$ pc = 560 au. We also find that the $I$($^{12}$CO(2-1))/$I$($^{13}$CO(2-1)) intensity ratio is very high ($>$21) in the flat-spectrum young stellar object, GY-51, located in the PDR. The enhancement of the ratios indicates that the UV radiation significantly affects the CO isotopologues via selective dissociation in the overall $\rho$ Ophiuchi A PDR, and that the $\rho$ Ophiuchi A PDR has a plane-parallel structure.

Herschel observations have emphasized the role of molecular filaments in star formation. However, the origin and evolution of these filaments are not yet well understood, partly because of the lack of kinematic information. To examine whether the B211/B213 filament is accreting background gas due to its gravitational potential, we produced a toy accretion model and compared its predictions to the 12CO(1--0) and 13CO(1--0) velocity patterns. We also examined the spatial distributions of Halpha, 857 GHz continuum, and HI emission to search for evidence of large-scale external effects. We estimated the depth of the cloud around the B211/B213 filament to be 0.3--0.7 pc under the assumption that the density of the gas is the same as the 13CO critical density. Compared to a linear extent of >10 pc in the plane of the sky, this suggests that the 3D morphology of the cloud is sheet-like. 12CO and 13CO PV diagrams perpendicular to the filament axis show that the emission from the gas surrounding B211/B213 is redshifted to the northeast of the filament and blueshifted to the southwest, respectively, and that the velocities of both components approach the filament velocity as the line of sight approaches the filament crest. The PV diagrams predicted by our accretion model are in good agreement with the observed 12CO and 13CO PV diagrams, supporting the scenario of mass accretion into the filament proposed by Palmeirim et al. Moreover, inspection of the distribution of the Halpha and 857 GHz emission in the Taurus-California-Perseus region suggests that the B211/B213 filament may have formed as a result of an expanding supershell generated by the Per OB2 association. Based on these results, we propose a scenario in which the B211/B213 filament was initially formed by large-scale compression of HI gas and then is now growing in mass due to the gravitational accretion of ambient cloud molecular gas.

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

We present results of the mapping observations covering a large area of 1 square degree around W40 and Serpens South carried out in the 12CO (J=1-0), 13CO (J=1-0), C18O (J=1-0), CCS (J_{N}=8_{7}-7_{6}), and N_2H+ (J=1-0) emission lines with the 45 m Nobeyama Radio Telescope. W40 is a blistered Hii region, and Serpens South is an infrared dark cloud accompanied by a young cluster. The relationship between these two regions which are separated by ~20' on the sky has not been clear so far. We found that the C18O emission is distributed smoothly throughout the W40 and Serpens South regions, and it seems that the two regions are physically connected. We divided the C18O emission into four groups in terms of the spatial distributions around the HII region which we call 5, 6, 7, and 8 km s^{-1} components according to their typical LSR velocities, and propose a three-dimensional model of the W40 and Serpens South complex. We found two elliptical structures in position-velocity diagrams, which can be explained as a part of two expanding shells. One of the shells is the small inner shell just around the HII region, and the other is the large outer shell corresponding to the boundary of the HII region. Dense gas associated with the young cluster of Serpens South is likely to be located at the surface of the outer shell, indicating that the natal clump of the young cluster is interacting with the outer shell being compressed by the expansion of the shell. We suggest that the expansion of the shell induced the formation of the young cluster.

We present Nobeyama 45m telescope C18O, 13CO, and 12CO(1-0) mapping observations towards an interstellar filament in the Taurus molecular cloud. We investigate the gas velocity structure along the filament and in its surrounding parent cloud. The filament is detected in the optically thin C18O emission as a single velocity component, ~1pc long, ~0.06pc wide structure. The C18O emission traces dust column densities larger than ~5x10^21 cm-2. The line-of-sight (LOS) velocity fluctuates along the filament crest with an average amplitude of ~0.2 km/s. The 13CO and 12CO integrated intensity maps show spatially extended emission around the elongated filament. We identify three extended structures with LOS velocities redshifted and blueshifted with respect to the average velocity of the filament identified in C18O. Based on combined analyses of velocity integrated channel maps and intensity variations of the optically thick 12CO spectra on and off the filament, we propose a 3-dimensional structure of the cloud surrounding the filament. We further suggest a multi-interaction scenario where sheet-like extended structures interact, in space and time, with the filament and are responsible for its compression and/or disruption, playing an important role in the star formation history of the filament. We also identify, towards the same field, a very faint filament showing a velocity field compatible with the filament formation process proposed by Inoue et al. (2017), where a filament is formed due to convergence of a flow of matter generated by the bending of the ambient magnetic field structure induced by an interstellar shock compression.

We present the discovery of expanding spherical shells around low to intermediate-mass young stars in the Orion A giant molecular cloud using observations of $^{12}$CO (1-0) and $^{13}$CO (1-0) from the Nobeyama Radio Observatory 45-meter telescope. The shells have radii from 0.05 to 0.85 pc and expand outward at 0.8 to 5 km/s. The total energy in the expanding shells is comparable to protostellar outflows in the region. Together, shells and outflows inject enough energy and momentum to maintain the cloud turbulence. The mass-loss rates required to power the observed shells are two to three orders of magnitude higher than predicted for line-driven stellar winds from intermediate-mass stars. This discrepancy may be resolved by invoking accretion-driven wind variability. We describe in detail several shells in this paper and present the full sample in the online journal.

The Orion Nebula Cluster toward the HII region M42 is the most outstanding young cluster at the smallest distance 410pc among the rich high-mass stellar clusters. By newly analyzing the archival molecular data of the 12CO(J=1-0) emission at 21" resolution, we identified at least three pairs of complementary distributions between two velocity components at 8km/s and 13km/s. We present a hypothesis that the two clouds collided with each other and triggered formation of the high-mass stars, mainly toward two regions including the nearly ten O stars, theta1 Ori and theta2 Ori, in M42 and the B star, NU Ori, in M43. The timescale of the collision is estimated to be ~0.1Myr by a ratio of the cloud size and velocity corrected for projection, which is consistent with the age of the youngest cluster members less than 0.1Myr. The majority of the low-mass cluster members were formed prior to the collision in the last one Myr. We discuss implications of the present hypothesis and the scenario of high-mass star formation by comparing with the other eight cases of triggered O star formation via cloud-cloud collision.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

We present an unbiased large-scale (9 deg$^2$) CN ($N$=1-0) and C$^{18}$O ($J$=1-0) survey of Cygnus-X conducted with the Nobeyama 45m Cygnus-X CO survey. CN and C$^{18}$O are detected in various objects towards the Cygnus-X North and South (e.g., DR17, DR18, DR21, DR22, DR23, and W75N). We find that CN/C$^{18}$O integrated intensity ratios are systematically different from region to region, and are especially enhanced in DR17 and DR18 which are irradiated by the nearby OB stars. This result suggests that CN/C$^{18}$O ratios are enhanced via photodissociation reactions. We investigate the relation between the CN/C$^{18}$O ratio and strength of the UV radiation field. As a result, we find that CN/C$^{18}$O ratios correlate with the far-UV intensities, $G_0$. We also find that CN/C$^{18}$O ratios decrease inside molecular clouds, where the interstellar UV radiation is reduced due to the interstellar dust extinction. We conclude that the CN/C$^{18}$O ratio is controlled by the UV radiation, and is a good probe of photon-dominated regions.

Recent studies have suggested that star formation (SF) in dense gas may be governed by essentially the same law in Galactic clouds and external galaxies. This conclusion remains controversial, however, because different tracers have been used to probe the dense gas mass in Galactic and extragalactic studies. We conducted observations of Aquila, Oph., and Orion B in HCN(1-0), HCO$^+$(1-0), and their isotopomers to calibrate the HCN and HCO$^+$ lines used as dense gas tracers in extragalactic studies and to test the possible universality of the star formation efficiency in dense gas, SFE$_{dense}$. H$^{13}$CO$^+$ and H$^{13}$CN were observed to be good tracers of the filaments detected with Herschel. Comparing the luminosities $L_{HCN}$ and $L_{HCO^+}$ with the reference masses $M_{Herschel}^{Av>8}$, the empirical conversion factors $\alpha_{Herschel-HCN}$ and $\alpha_{Herschel-HCO^+}$ were found to be anti-correlated with the local FUV strength. In agreement with Pety et al. (2017), HCN and HCO$^+$ were also found to trace gas down to Av>2. As a result, published extragalactic HCN studies must be tracing all of the moderate density gas. Estimating the contribution of this moderate density gas from the typical column density PDFs in nearby clouds, we obtained the following $G_0$-dependent HCN conversion factor for external galaxies:$\alpha_{Herschel-HCN}^{fit'}=64\times G_0^{-0.34}$. Re-estimating the dense gas masses in external galaxies with $\alpha_{Herschel-HCN}^{fit'}$, we found that SFE$_{dense}$ is remarkably constant over 8 orders of magnitude in dense gas mass. Our results confirm that SFE$_{dense}$ of galaxies is quasi-universal on a wide range of scales from ~1-10pc to >10kpc. Based on the tight link between SF and filamentary structure found in Herschel studies of nearby clouds, we argue that SFE$_{dense}$ is primarily set by the microphysics of SF in filaments.

Stars form out of the densest parts of molecular clouds. Far-IR emission can be used to estimate the Star Formation Rate (SFR) and high dipole moment molecules, typically HCN, trace the dense gas. A strong correlation exists between HCN and Far-IR emission, with the ratio being nearly constant, over a large range of physical scales. A few recent observations have found HCN to be weak with respect to the Far-IR and CO in subsolar metallicity (low-Z) objects. We present observations of the Local Group galaxies M33, IC10, and NGC6822 with the IRAM 30meter and NRO 45m telescopes, greatly improving the sample of low-Z galaxies observed. HCN, HCO$^+$, CS, C$_2$H, and HNC have been detected. Compared to solar metallicity galaxies, the Nitrogen-bearing species are weak (HCN, HNC) or not detected (CN, HNCO, N$_2$H$^+$) relative to Far-IR or CO emission. HCO$^+$ and C$_2$H emission is normal with respect to CO and Far-IR. While $^{13}$CO is the usual factor 10 weaker than $^{12}$CO, C$^{18}$O emission was not detected down to very low levels. Including earlier data, we find that the HCN/HCO$^+$ ratio varies with metallicity (O/H) and attribute this to the sharply decreasing Nitrogen abundance. The dense gas fraction, traced by the HCN/CO and HCO$^+$/CO ratios, follows the SFR but in the low-Z objects the HCO$^+$ is much easier to measure. Combined with larger and smaller scale measurements, the HCO$^+$ line appears to be an excellent tracer of dense gas and varies linearly with the SFR for both low and high metallicities.

We report a new 1-pc (30") resolution CS($J=2-1$) line map of the central 30 pc of the Galactic Center (GC), made with the Nobeyama 45m telescope. We revisit our previous study of the extraplanar feature called polar arc (PA), which is a molecular cloud located above SgrA* with a velocity gradient perpendicular to the Galactic plane. We find that the PA can be traced back to the Galactic disk. This provides clues of the launching point of the PA , roughly $6\times10^{6}$ years ago. Implications of the dynamical time scale of the PA might be related to the Galactic Center Lobe (GCL) at parsec scale. Our results suggest that in the central 30 pc of the GC, the feedback from past explosions could alter the orbital path of the molecular gas down to the central tenth of parsec. In the follow-up work of our new CS($J=2-1$) map, we also find that near the systemic velocity, the molecular gas shows an extraplanar hourglass-shaped feature (HG-feature) with a size of $\sim$13 pc. The latitude-velocity diagrams show that the eastern edge of the HG-feature is associated with an expanding bubble B1, $\sim$7 pc away from SgrA*. The dynamical time scale of this bubble is $\sim3\times10^{5}$ years. This bubble is interacting with the 50 km s$^{-1}$ cloud. Part of the molecular gas from the 50 km s$^{-1}$ cloud was swept away by the bubble to $b=-0.2deg$. The western edge of the HG-feature seems to be the molecular gas entrained from the 20 km s$^{-1}$ cloud towards the north of the Galactic disk. Our results suggest a fossil explosion in the central 30 pc of the GC a few 10$^{5}$ years ago.

Herschel observations of nearby molecular clouds suggest that interstellar filaments and prestellar cores represent two fundamental steps in the star formation process. The observations support a picture of low-mass star formation according to which ~ 0.1 pc-wide filaments form first in the cold interstellar medium, probably as a result of large-scale compression of interstellar matter by supersonic turbulent flows, and then prestellar cores arise from gravitational fragmentation of the densest filaments. Whether this scenario also applies to regions of high-mass star formation is an open question, in part because Herschel data cannot resolve the inner width of filaments in the nearest regions of massive star formation. We used the bolometer camera ArTeMiS on the APEX telescope to map the central part of the NGC6334 complex at a factor of > 3 higher resolution than Herschel at 350 microns. Combining ArTeMiS data with Herschel data allowed us to study the structure of the main filament of the complex with a resolution of 8" or < 0.07 pc at d ~ 1.7 kpc. Our study confirms that this filament is a very dense, massive linear structure with a line mass ranging from ~ 500 Msun/pc to ~ 2000 Msun/pc over nearly 10 pc. It also demonstrates that its inner width remains as narrow as W ~ 0.15 +- 0.05 pc all along the filament length, within a factor of < 2 of the characteristic 0.1 pc value found with Herschel for lower-mass filaments in the Gould Belt. While it is not completely clear whether the NGC 6334 filament will form massive stars or not in the future, it is two to three orders of magnitude denser than the majority of filaments observed in Gould Belt clouds, and yet has a very similar inner width. This points to a common physical mechanism for setting the filament width and suggests that some important structural properties of nearby clouds also hold in high-mass star forming regions.

L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.

We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $\alpha_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.

We performed the first spectral-line survey at 82--106 GHz and 335--355 GHz toward the outflow-shocked region, OMC 2-FIR 4, the outflow driving source, FIR 3, and the northern outflow lobe, FIR 3N. We detected 120 lines of 20 molecular species. The line profiles are found to be classifiable into two types: one is a single Gaussian component with a narrow ($<$ 3 km s$^{-1}$) width and another is two Gaussian components with narrow and wide ($>$ 3km s$^{-1}$) widths. The narrow components for the most of the lines are detected at all positions, suggesting that they trace the ambient dense gas. For CO, CS, HCN, and HCO$^{+}$, the wide components are detected at all positions, suggesting the outflow origin. The wide components of C$^{34}$S, SO, SiO, H$^{13}$CN, HC$^{15}$N, H$_2^{13}$CO, H$_2$CS, HC$_3$N, and CH$_3$OH are detected only at FIR 4, suggesting the outflow-shocked gas origin. The rotation diagram analysis revealed that the narrow components of C$_2$H and H$^{13}$CO$^+$ show low temperatures of 12.5$\pm$1.4 K, while the wide components show high temperatures of 20--70 K. This supports our interpretation that the wide components trace the outflow and/or outflow-shocked gas. We compared observed molecular abundances relative to H$^{13}$CO$^+$ with those of the outflow-shocked region, L1157 B1, and the hot corino, IRAS 16293-2422. Although we cannot exclude a possibility that the chemical enrichment in FIR 4 is caused by the hot core chemistry, the chemical compositions in FIR 4 are more similar to those in L1157 B1 than those in IRAS 16293-2422.

Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^{\alpha})$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $\sigma_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $\sigma_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

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

We report results of the CO(J=3-2) and HCO+(J=4-3) observations of the W40 HII region with the ASTE 10 m telescope (HPBW~22 arcsec) to search for molecular outflows and dense clumps.We found that the velocity field in the region is highly complex, consisting of at least four distinct velocity components at V LSR ~ 3, 5, 7, and 10 km/s. The ~7 km/s component represents the systemic velocity of cold gas surroundingthe entire region, and causes heavy absorption in the CO spectra over the velocity range 6 <V LSR< 9 km/s. The ~5 and ~10 km/s components exhibit high CO temperature (>40 K) and are found mostly around the HII region, suggesting that these components are likely to be tracing dense gas interacting with the expanding shell around the HII region. Based on the CO data, we identified 13 regions of high velocity gas which we interpret as candidate outflow lobes. Using the HCO+ data, we also identified six clumps and estimated their physical parameters. On the basis of the ASTE data and near-infrared images from 2MASS, we present an updated three-dimensional model of this region. In order to investigate molecular outflows in W40, the SiO (J=1-0, v=0) emission line and some other emission lines at 40 GHz were also observed with the 45 m telescope at the Nobeyama Radio Observatory, but they were not detected at the present sensitivity.

We performed single point [C I] $^3$P$_1$-$^3$P$_0$ and CO J=4-3 observations toward three T Tauri stars, DM Tau, LkCa 15, and TW Hya, using the Atacama Large Millimeter/submillimeter Array (ALMA) Band 8 qualification model receiver installed on the Atacama Submillimeter Telescope Experiment (ASTE). Two protostars in the Taurus L1551 region, L1551 IRS 5 and HL Tau, were also observed. We successfully detected [C I] emission from the protoplanetary disk around DM Tau as well as the protostellar targets. The spectral profile of the [C I] emission from the protoplanetary disk is marginally single-peaked, suggesting that atomic carbon (C) extends toward the outermost disk. The detected [C I] emission is optically thin and the column densities of C are estimated to be <~10$^{16}$ cm$^{-2}$ and ~10$^{17}$ cm$^{-2}$ for the T Tauri star targets and the protostars, respectively. We found a clear difference in the total mass ratio of C to dust, $M$(C)/$M$(dust), between the T Tauri stars and protostellar targets; the $M$(C)/$M$(dust) ratio of the T Tauri stars is one order of magnitude smaller than that of the protostars. The decrease of the estimated $M$(C)/$M$(dust) ratios for the disk sources is consistent with a theoretical prediction that the atomic C can survive only in the near surface layer of the disk and C$^+$/C/CO transition occurs deeper into the disk midplane.

We present Orion A giant molecular cloud core catalogs, which are based on 1.1 mm map with an angular resolution of 36 arcsec (sim 0.07 pc) and C18O (1-0) data with an angular resolution of 26.4 arcsec (sim 0.05 pc). We have cataloged 619 dust cores in the 1.1 mm map using the Clumpfind method. The ranges of the radius, mass, and density of these cores are estimated to be 0.01 - 0.20 pc, 0.6 - 1.2 times 10^2 Msun, and 0.3 times 10^4 - 9.2 times 10^6 cm^{-3}, respectively. We have identified 235 cores from the C18O data. The ranges of the radius, velocity width, LTE mass, and density are 0.13 -- 0.34 pc, 0.31 - 1.31 km s^{-1}, 1.0 - 61.8 Msun, and (0.8 - 17.5) times 10^3 cm^{-3}, respectively. From the comparison of the spatial distributions between the dust and C18O cores, four types of spatial relations were revealed: (1) the peak positions of the dust and C18O cores agree with each other (32.4% of the C18O cores), (2) two or more C18O cores are distributed around the peak position of one dust core (10.8% of the C18O cores), (3) 56.8% of the C18O cores are not associated with any dust cores, and (4) 69.3% of the dust cores are not associated with any C18O cores. The data sets and analysis are public.

The Serpens South infrared dark cloud consists of several filamentary ridges, some of which fragment into dense clumps. On the basis of CCS ($J_N=4_3-3_2$), HC$_3$N ($J=5-4$), N$_2$H$^+$ ($J=1-0$), and SiO ($J=2-1, v=0$) observations, we investigated the kinematics and chemical evolution of these filamentary ridges. We find that CCS is extremely abundant along the main filament in the protocluster clump. We emphasize that Serpens South is the first cluster-forming region where extremely-strong CCS emission is detected. The CCS-to-N$_2$H$^+$ abundance ratio is estimated to be about 0.5 toward the protocluster clump, whereas it is about 3 in the other parts of the main filament. We identify six dense ridges with different $V_{\rm LSR}$. These ridges appear to converge toward the protocluster clump, suggesting that the collisions of these ridges may have triggered cluster formation. The collisions presumably happened within a few $\times \ 10^5$ yr because CCS is abundant only in such a short time. The short lifetime agrees with the fact that the number fraction of Class I objects, whose typical lifetime is $0.4 \times \ 10^5$ yr, is extremely high as about 70 percent in the protocluster clump. In the northern part, two ridges appear to have partially collided, forming a V-shape clump. In addition, we detected strong bipolar SiO emission that is due to the molecular outflow blowing out of the protostellar clump, as well as extended weak SiO emission that may originate from the filament collisions.

L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC/ASTE telescope, which we combined with archive data from near- and mid-infrared (2MASS, WISE) to mm observations (Planck). The Herschel observations, combined with other data, show a sequence of objects from a cold clump to young stellar objects at different evolutionary stages. Source B-3 (2MASS J04351455-1414468) appears to be a YSO forming inside the L1642 cloud, instead of a foreground brown dwarf, as previously classified. Herschel data reveal striation in the diffuse dust emission around L1642. The western region shows striation towards NE and has a steeper column density gradient on its southern side. The densest central region has a bow-shock like structure showing compression from the west and a filamentary tail extending towards east. The differences suggest that these may be spatially distinct structures, aligned only in projection. We derive values of the dust emission cross-section per H nucleon for different regions of the cloud. Modified black-body fits to the spectral energy distribution of Herschel and Planck data give emissivity spectral index beta values 1.8-2.0 for the different regions. The compact sources have lower beta values and show an anticorrelation between T and beta. Markov chain Monte Carlo calculations demonstrate the strong anticorrelation between beta and T errors and the importance of mm Planck data in constraining the estimates. L1642 reveals a more complex structure and sequence of star formation than previously known.

Aims. We derive physical properties such as the optical depths and the column densities of $^{13}$CO and C$^{18}$O to investigate the relationship between the far ultraviolet (FUV) radiation and the abundance ratios between $^{13}$CO and C$^{18}$O. Method. We have carried out wide-field (0.4 deg$^2$) observations with an angular resolution of 25.8 arcsec ($\sim$ 0.05 pc) in $^{13}$CO ($J$=1--0) and C$^{18}$O ($J$=1--0) toward the Orion-A giant molecular cloud using the Nobeyama 45 m telescope in the on-the-fly mode. Results. Overall distributions and velocity structures of the $^{13}$CO and C$^{18}$O emissions are similar to those of the $^{12}$CO ($J$=1--0) emission. The optical depths of the $^{13}$CO and C18O emission lines are estimated to be 0.05 $<$ $\tau_{\rm ^{13}CO}$ $<$ 1.54 and 0.01 $<$ $\tau_{\rm C^{18}O}$ $<$ 0.18, respectively. The column densities of the $^{13}$CO and C$^{18}$O emission lines are estimated to be 0.2 $\times$ 10$^{16}$ $<$ $N_{\rm ^{13}CO}$ $<$ 3.7 $\times$ 10$^{17}$ cm$^{-2}$ and 0.4 $\times$ 10$^{15}$ $<$ $N_{\rm C^{18}O}$ $<$ 3.5 $\times$ 10$^{16}$ cm$^{-2}$, respectively. The abundance ratios between $^{13}$CO and C$^{18}$O, $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$, are found to be 5.7 - 33.0. The mean value of $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$ in the nearly edge-on photon-dominated regions is found to be 16.47 $\pm$ 0.10, which is a third larger than that the solar system value of 5.5. The mean value of $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$ in the other regions is found to be 12.29 $\pm$ 0.02. The difference of the abundance ratio is most likely due to the selective FUV photodissociation of C$^{18}$O.

We present the results of N$_2$H$^+$ ($J=1-0$) observations toward Serpens South, the nearest cluster-forming, infrared dark cloud. The physical quantities are derived by fitting the hyperfine structure of N$_2$H$^+$. The Herschel and 1.1-mm continuum maps show that a pc-scale filament fragments into three clumps with radii of $0.1-0.2$ pc and masses of $40-230M_\odot$. We find that the clumps contain smaller-scale ($\sim 0.04$ pc) structures, i.e., dense cores. We identify 70 cores by applying CLUMPFIND to the N$_2$H$^+$ data cube. In the central cluster-forming clump, the excitation temperature and line-width tend to be large, presumably due to protostellar outflow feedback and stellar radiation. However, for all the clumps, the virial ratios are evaluated to be $0.1-0.3$, indicating that the internal motions play only a minor role in the clump support. The clumps exhibit no free-fall, but low-velocity infall, and thus the clumps should be supported by additional forces. The most promising force is the globally-ordered magnetic field observed toward this region. We propose that the Serpens South filament was close to magnetically-critical and ambipolar diffusion triggered the cluster formation. We find that the northern clump, which shows no active star formation, has a mass and radius comparable to the central cluster-forming clump, and therefore, it is a likely candidate of a {\it pre-protocluster clump}. The initial condition for cluster formation is likely to be a magnetically-supported clump of cold, quiescent gas. This appears to contradict the accretion-driven turbulence scenario, for which the turbulence in the clumps is maintained by the accretion flow.

We have carried out wide-field (0.17 degree^2) and high-angular resolution (21.3 arcsec ~ 0.04 pc) observations in [CI] line toward the Orion-A giant molecular cloud with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in the On-The-Fly (OTF) mode. Overall features of the [CI] emission are similar to those of the CO (1--0) emission in Shimajiri et al. (2011); the total intensity ratio of the [CI] to CO emission ranges from 0.05 to 0.2. The optical depth of the [CI] emission is found to be 0.1 -- 0.75, suggesting optically thin emission. The column density of the [CI] emission is estimated to be (1.0 -- 19) x 10^17 cm^-2. These results are consistent with the results of the previous [CI] observations with a low-angular resolution of 2.2 arcmin (e.g. Ikeda et al. 1999). In the nearly edge-on PDRs and their candidates of the Orion Bar, DLSF, M 43 Shell, and Region D, the distributions of the [CI] emission coincide with those of the CO emission, inconsistent with the prediction by the plane-parallel PDR model (Hollenbach & Tielens 1999). In addition, the [CI] distribution in the Orion A cloud is found to be more similar to those of the ^{13}CO (1--0), C^{18}O (1--0), and H^{13}CO^+ (1--0) lines than that of the CO (1--0) line, suggesting that the [CI] emission is not limited to the cloud surface, but is tracing the dense, inner parts of the cloud.

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

OMC-2 FIR 4 is one of the closest known young intermediate-mass protoclusters, located at a distance of 420 pc in Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and IRAM 30-m telescopes. In order to investigate the morphology of this region, we used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as in DCO+(2-1), DCN(2-1), C34S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us to resolve structures of 5", equivalent to 2000 AU. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified displays a DCO+(2-1) emission peak, and weak dust continuum emission. Its higher abundance of DCO+ relative to the other two regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. It may alternatively be part of the colder envelope of OMC-2 FIR 4. Our interferometric observations evidence the complexity of this region, where multiple cores, chemical differentiation and an ionised region all coexist within an area of only 10000 AU.

We present high spatial resolution $^{12}$CO ($J=1-0$) images taken by the Nobeyama 45m telescope toward a $48' \times 48'$ area including the L1641-N cluster. The effective spatial resolution of the maps is $21"$, corresponding to 0.04 pc at a distance of 400 pc. A recent 1.1 mm dust continuum map reveals that the dense gas is concentrated in several thin filaments. We find that a few dust filaments are located at the parts where $^{12}$CO ($J=1-0$) emission drops sharply. Furthermore, the filaments have two-components with different velocities. The velocity difference between the two-components is about 3 km s$^{-1}$, corresponding to a Mach number of 10, significantly larger than the local turbulent velocity in the cloud. These facts imply that the collision of the two components (hereafter, the cloud-cloud collision) possibly contributed to the formation of these filaments. Since the two components appear to overlap toward the filaments on the plane of the sky, the collision may have occurred almost along the line of sight. Star formation in the L1641-N cluster was probably triggered by such a collision. We also find several parsec-scale CO shells whose centers are close to either the L1641-N cluster or V 380 Ori cluster. We propose that these shells were created by multiple winds and/or outflows from cluster YSOs, i.e., "protocluster winds." One exceptional dust filament located at the western cloud edge lies along a shell; it is presumably a part of the expanding shell. Both the cloud-cloud collision and protocluster winds are likely to influence the cloud structure and kinematics in this region.

We present the results of CO ($J=3-2$) and HCO$^+$ ($J=4-3$) mapping observations toward a nearby embedded cluster, Serpens South, using the ASTE 10 m telescope. Our CO ($J=3-2$) map reveals that many outflows are crowded in the dense cluster-forming clump that can be recognized as a HCO$^+$ clump with a size of $\sim$ 0.2 pc and mass of $\sim$ 80 M$_\odot$. The clump contains several subfragments with sizes of $\sim$ 0.05 pc. By comparing the CO ($J=3-2$) map with the 1.1 mm dust continuum image taken by AzTEC on ASTE, we find that the spatial extents of the outflow lobes are sometimes anti-correlated with the distribution of the dense gas and some of the outflow lobes apparently collide with the dense gas. The total outflow mass, momentum, and energy are estimated at 0.6 $M_\odot$, 8 $M_\odot$ km s$^{-1}$, and 64 $M_\odot$ km$^2$ s$^{-2}$, respectively. The energy injection rate due to the outflows is comparable to the turbulence dissipation rate in the clump, implying that the protostellar outflows can maintain the supersonic turbulence in this region. The total outflow energy seems only about 10 percent the clump gravitational energy. We conclude that the current outflow activity is not enough to destroy the whole cluster-forming clump, and therefore star formation is likely to continue for several or many local dynamical times.

We present the discovery of a cold massive dust disk around the T Tauri star V1094 Sco in the Lupus molecular cloud from the 1.1 millimeter continuum observations with AzTEC on ASTE. A compact ($r\lesssim$320 AU) continuum emission coincides with the stellar position having a flux density of 272 mJy which is largest among T Tauri stars in Lupus. We also present the detection of molecular gas associated with the star in the five-point observations in $^{12}$CO J=3--2 and $^{13}$CO J=3--2. Since our $^{12}$CO and $^{13}$CO observations did not show any signature of a large-scale outflow or a massive envelope, the compact dust emission is likely to come from a disk around the star. The observed SED of V1094 Sco shows no distinct turnover from near infrared to millimeter wavelengths, which can be well described by a flattened disk for the dust component, and no clear dip feature around 10 $\micron$ suggestive of absence of an inner hole in the disk. We fit a simple power-law disk model to the observed SED. The estimated disk mass ranges from 0.03 to $\gtrsim$0.12 $M_\sun$, which is one or two orders of magnitude larger than the median disk mass of T Tauri stars in Taurus.

We present new, wide and deep images in the 1.1 mm continuum and the $^{12}$CO ($J$=1-0) emission toward the northern part of the Orion A Giant Molecular Cloud (Orion-A GMC). The 1.1 mm data were taken with the AzTEC camera mounted on the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in Chile, and the $^{12}$CO ($J$=1-0) data were with the 25 beam receiver (BEARS) on the NRO 45 m telescope in the On-The-Fly (OTF) mode. The present AzTEC observations are the widest $(\timeform{1.D7}$ $\times$ $\timeform{2.D3}$, corresponding to 12 pc $\times$ 17 pc) and the highest-sensitivity ($\sim$9 mJy beam$^{-1}$) 1.1 mm dust-continuum imaging of the Orion-A GMC with an effective spatial resolution of $\sim$ 40$\arcsec$. The $^{12}$CO ($J$=1-0) image was taken over the northern $\timeform{1D.2} \times\timeform{1D.2}$ (corresponding 9 pc $\times$ 9 pc) area with a sensitivity of 0.93 K in $T_{\rm MB}$, a velocity resolution of 1.0 km s$^{-1}$, and an effective spatial resolution of 21$\arcsec$. With these data, together with the MSX 8 $\mu$m, Spitzer 24 $\mu$m and the 2MASS data, we have investigated the detailed structure and kinematics of molecular gas associated with the Orion-A GMC and have found evidence for interactions between molecular clouds and the external forces that may trigger star formation. Two types of possible triggers were revealed; 1) Collision of the diffuse gas on the cloud surface, particularly at the eastern side of the OMC-2/3 region, and 2) Irradiation of UV on the pre-existing filaments and dense molecular cloud cores. Our wide-field and high-sensitivity imaging have provided the first comprehensive view of the potential sites of triggered star formation in the Orion-A GMC.

We present the results of CO (J=3-2) and CO (J=1-0) mapping observations toward the active cluster forming clump, L1688, in the rho Ophiuchi molecular cloud. From the CO (J=3-2) and CO (J=1-0) data cubes, we identify five outflows, whose driving sources are VLA 1623, EL 32, LFAM 26, EL 29, and IRS 44. Among the identified outflows, the most luminous outflow is the one from the prototypical Class 0 source, VLA 1623. We also discover that the EL 32 outflow located in the Oph B2 region has very extended blueshifted and redshifted lobes with wide opening angles. This outflow is most massive and have the largest momentum among the identified outflows in the CO (J=1-0) map. We estimate the total energy injection rate due to the molecular outflows identified by the present and previous studies to be about 0.2 L_solar, larger than or at least comparable to the turbulence dissipation rate [~(0.03 - 0.1) L_solar]. Therefore, we conclude that the protostellar outflows are likely to play a significant role in replenishing the supersonic turbulence in this clump.

We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Compact Array pointed observations of N2H+ 1-0 emission towards the clustered, low mass star forming Oph B Core within the Ophiuchus molecular cloud. We compare these data with previously published results of high resolution NH3 (1,1) and (2,2) observations in Oph B. We use 3D Clumpfind to identify emission features in the single-dish N2H+ map, and find that the N2H+ `clumps' match well similar features previously identified in NH3 (1,1) emission, but are frequently offset to clumps identified at similar resolution in 850 micron continuum emission. Wide line widths in the Oph B2 sub-Core indicate non-thermal motions dominate the Core kinematics, and remain transonic at densities n ~ 3 x 10^5 cm^-3 with large scatter and no trend with N(H2). Non-thermal motions in Oph B1 and B3 are subsonic with little variation, but also show no trend with H2 column density. Over all Oph B, non-thermal N2H+ line widths are substantially narrower than those traced by NH3, making it unlikely NH3 and N2H+ trace the same material, but the v_LSR of both species agree well. We find evidence for accretion in Oph B1 from the surrounding ambient gas. The NH3/N2H+ abundance ratio is larger towards starless Oph B1 than towards protostellar Oph B2, similar to recent observational results in other star-forming regions. Small-scale structure is found in the ATCA N2H+ 1-0 emission, where emission peaks are again offset from continuum emission. In particular, the ~1 M_Sun B2-MM8 clump is associated with a N2H+ emission minimum and surrounded by a broken ring-like N2H+ emission structure, suggestive of N2H+ depletion. We find a strong general trend of decreasing N2H+ abundance with increasing N(H2) in Oph B which matches that found for NH3.

We performed simultaneous multiwavelength observations of OJ 287 with the Nobeyama Millimeter Array for radio, the KANATA telescope and the KVA telescope for optical, the Suzaku satellite for X-ray and the MAGIC telescope for very high energy (VHE) gamma-ray in 2007. The observations were conducted for a quiescent state in April and in a flaring state in November-December. We clearly observed increase of fluxes from radio to X-ray bands during the flaring state while MAGIC could not detect significant VHE gamma-ray emission from the source. We could derive an upper limit (95% confidence level) of 1.7% of the Crab Nebula flux above 150 GeV from about 41.2 hours of the MAGIC observation. A simple SSC model suggests that the observed flaring activity could be caused by evolutions in the distribution of the electron population rather than changes of the magnetic field strength or Doppler beaming factor in the jet.

We have conducted millimeter interferometric observations of the Orion Molecular Cloud-2 (OMC-2) FIR 6 region at an angular resolution of $\sim$ 4$\arcsec$ - 7$\arcsec$ with the Nobeyama Millimeter Array (NMA). In the 3.3 mm continuum emission we detected dusty core counterparts of the previously identified FIR sources (FIR 6a, 6b, 6c, and 6d), and moreover, resolved FIR 6a into three dusty cores. The size and mass of these cores are estimated to be 1100-5900 AU and 0.19-5.5 M$_{\odot}$, respectively. We found that in the $^{12}$CO ($J$=1--0) emission FIR 6b, 6c, and 6d eject the molecular outflow and that the FIR 6c outflow also exhibits at least two collimated jet-like components in the SiO ($J$=2--1) emission. At the tip of one of the SiO components there appears abrupt increase of the SiO line width ($\sim$ 15 km s$^{-1}$), where the three resolved cores in FIR 6a seem to delineate the tip. These results imply the presence of the interaction and the bowshock front between the FIR 6c molecular outflow and FIR 6a.

Suzaku observations of the blazar OJ 287 were performed in 2007 April 10--13 and November 7--9. They correspond to a quiescent and a flaring state, respectively. The X-ray spectra can be well described with single power-law models in both exposures. The derived X-ray photon index and the flux density at 1 keV were found to be Gamma = 1.65 +- 0.02 and S_{1 keV} = 215 +- 5 nJy, in the quiescent state. In the flaring state, the source exhibited a harder X-ray spectrum (Gamma = 1.50 +- 0.01) with a nearly doubled X-ray flux density S_{1 keV} = 404^{+6}_{-5} nJy. Moreover, significant hard X-ray signals were detected up to ~ 27 keV. In cooperation with the Suzaku, simultaneous radio, optical, and very-high-energy gamma-ray observations were performed with the Nobeyama Millimeter Array, the KANATA telescope, and the MAGIC telescope, respectively. The radio and optical fluxes in the flaring state (3.04 +- 0.46 Jy and 8.93 +- 0.05 mJy at 86.75 Hz and in the V-band, respectively) were found to be higher by a factor of 2--3 than those in the quiescent state (1.73 +- 0.26 Jy and 3.03 +- 0.01 mJy at 86.75 Hz and in the V-band, respectively). No notable gamma-ray events were detected in either observation. The spectral energy distribution indicated that the X-ray spectrum was dominated by inverse Compton radiation in both observations, while synchrotron radiation exhibited a spectral cutoff around the optical frequency. Furthermore, no significant difference in the synchrotron cutoff frequency was found between the quiescent and flaring states. According to a simple synchrotron self-Compton model, the change of the spectral energy distribution is due to an increase in the energy density of electrons with small changes of both the magnetic field strength and the maximum Lorentz factor of electrons.

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

We have carried out the observations of the OMC-2 FIR 3/4 region with the NMA and ASTE in the H$^{13}$CO$^{+}$ (1--0), $^{12}$CO (3--2, 1--0), SiO ($v$=0, $J$=2--1), CS (2--1), and CH$_3$OH ($J_K$=7$_K$--6$_K$) lines and in the 3.3 mm continuum emission. Our NMA observations in the H$^{13}$CO$^{+}$ emission have revealed 0.07 pc-scale dense gas associated with FIR 4. The $^{12}$CO (3--2,1--0) emission shows high-velocity blue and red shifted components at the both north-east and south-west of FIR 3, suggesting a molecular outflow nearly along the plane of the sky driven by FIR 3. The SiO and the CH$_{3}$OH emission are detected around the interface between the outflow and the dense gas. Furthermore, the $^{12}$CO (1--0) emission shows an L-shaped structure in the P-V diagram. These results imply presence of the shock due to the interaction between the molecular outflow driven by FIR 3 and the dense gas associated with FIR 4. Moreover, our high angular-resolution observations of FIR 4 in the 3.3 mm continuum emission have first found that FIR 4 consists of eleven dusty cores. The separation among these cores is on the same order of the Jeans length, suggesting that the fragmentation into these cores has been caused by the gravitational instability. The time scale of the fragmentation is similar to the time scale of the interaction between the molecular outflow and the dense gas. We suggest that the interaction between the molecular outflow from FIR 3 and the dense gas associated with FIR 4 triggered the fragmentation into these dusty cores, and hence the next generation the cluster formation.