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 introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $\Sigma$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$\mu$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By comparing with the ATLASGAL catalog, we found that the majority of the $^{13}$CO clouds with 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$ are likely located at distances of 2 kpc (Sagittarius arm) and 3 kpc (Scutum arm), respectively. Analyzing the spatial configuration of the identified clouds and their velocity structures, we attempt to reveal the origin of the cloud structure in this region. Here we discuss three possibilities: (1) overlapping with different velocities, (2) cloud oscillation, and (3) cloud-cloud collision. From the position-velocity diagrams, we found spatially-extended faint emission between $\sim$ 20 km s$^{-1}$ and $\sim$ 35 km s$^{-1}$, which is mainly distributed in the spatially-overlapped areas of the clouds. We also found that in some areas where clouds with different velocities overlapped, the magnetic field orientation changes abruptly. The distribution of the diffuse emission in the position-position-velocity space and the bending magnetic fields appear to favor the cloud-cloud collision scenario compared to other scenarios. In the cloud-cloud collision scenario, we propose that two $\sim$35 km s$^{-1}$ foreground clouds are colliding with clouds at $\sim$20 km s$^{-1}$ with a relative velocity of 15 km s$^{-1}$. These clouds may be substructures of two larger clouds having velocities of $\sim$ 35 km s$^{-1}$ ($\gtrsim 10^3 $ M$_{\odot}$) and $\sim$ 20 km s$^{-1}$ ($\gtrsim 10^4 $ M$_{\odot}$), respectively.

We present maps in several molecular emission lines of a 1 square-degree region covering the W40 and Serpens South molecular clouds belonging to the Aquila Rift complex. The observations were made with the 45 m telescope at the Nobeyama Radio Observatory. We found that the $^{12}$CO and $^{13}$CO emission lines consist of several velocity components with different spatial distributions. The component that forms the main cloud of W40 and Serpens South, which we call the "main component", has a velocity of $V_{\rm LSR}\simeq 7$ km s$^{-1}$. There is another significant component at $V_{\rm LSR}\simeq 40$ km s$^{-1}$, which we call the "40 km s$^{-1}$ component". The latter component is mainly distributed around two young clusters: W40 and Serpens South. Moreover, the two components look spatially anti-correlated. Such spatial configuration suggests that the star formation in W40 and Serpens South was induced by the collision of the two components. We also discuss a possibility that the 40 km s$^{-1}$ component consists of gas swept up by superbubbles created by SNRs and stellar winds from the Scorpius-Centaurus Association.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

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.

We investigate the formation of stars within giant molecular clouds (GMCs) evolving in environments of different global magnetic field strength and large-scale dynamics. Building upon a series of magnetohydrodynamic (MHD) simulations of non-colliding and colliding GMCs, we employ density- and magnetically-regulated star formation sub-grid models in clouds which range from moderately magnetically supercritical to near critical. We examine gas and star cluster morphologies, magnetic field strengths and relative orientations, pre-stellar core densities, temperatures, mass-to-flux ratios and velocities, star formation rates and efficiencies over time, spatial clustering of stars, and kinematics of the stars and natal gas. The large scale magnetic criticality of the region greatly affects the overall gas evolution and star formation properties. GMC collisions enhance star formation rates and efficiencies in magnetically supercritical conditions, but may actually inhibit them in the magnetically critical case. This may have implications for star formation in different Galactic environments such as the Galactic Center and the main Galactic disk.

The non-uniform distribution of gas and protostars in molecular clouds is caused by combinations of various physical processes that are difficult to separate. We explore this non-uniform distribution in the M17 molecular cloud complex that hosts massive star formation activity using the $^{12}$CO ($J=1-0$) and $^{13}$CO ($J=1-0$) emission lines obtained with the Nobeyama 45m telescope. Differences in clump properties such as mass, size, and gravitational boundedness reflect the different evolutionary stages of the M17-H{\scriptsize II} and M17-IRDC clouds. Clumps in the M17-H{\scriptsize II} cloud are denser, more compact, and more gravitationally bound than those in M17-IRDC. While M17-H{\scriptsize II} hosts a large fraction of very dense gas (27\%) that has column density larger than the threshold of $\sim$ 1 g cm$^{-2}$ theoretically predicted for massive star formation, this very dense gas is deficient in M17-IRDC (0.46\%). Our HCO$^+$ ($J=1-0$) and HCN ($J=1-0$) observations with the TRAO 14m telescope, {\nlqb trace all gas with column density higher than $3\times 10^{22}$ cm$^{-2}$}, confirm the deficiency of high density ($\gtrsim 10^5$ cm$^{-3}$) gas in M17-IRDC. Although M17-IRDC is massive enough to potentially form massive stars, its deficiency of very dense gas and gravitationally bound clumps can explain the current lack of massive star formation.

We have carried out mapping observations of molecular emission lines of HC$_{3}$N and CH$_{3}$OH toward two massive cluster-forming clumps, NGC2264-C and NGC2264-D, using the Nobeyama 45-m radio telescope. We derive an $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) integrated intensity ratio map, showing a higher value at clumps including 2MASS point sources at the northern part of NGC2264-D. Possible interpretations of the $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) ratio are discussed. We have also observed molecular emission lines from CCS and N$_{2}$H$^{+}$ toward five positions in each clump. We investigate the $N$(N$_{2}$H$^{+}$)/$N$(CCS) and $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) column density ratios among the ten positions in order to test whether they can be used as chemical evolutionary indicators in these clumps. The $N$(N$_{2}$H$^{+}$)/$N$(CCS) ratio shows a very high value toward a bright embedded IR source (IRS1), whereas the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) ratio at IRS1 is comparable with those at the other positions. These results suggest that UV radiation affects the chemistry around IRS1. We find that there are positive correlations between these column density ratios and the excitation temperatures of N$_{2}$H$^{+}$, which implies the chemical evolution of clumps. These chemical evolutionary indicators likely reflect the combination of evolution along the filamentary structure and evolution of each clump.

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

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

We 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 report a first clear detection of the Zeeman splitting of a CCS emission line at 45 GHz toward a nearby prestellar dense filament, Taurus Molecular Cloud-1. We observed HC$_3$N non-Zeeman line simultaneously as the CCS line, and did not detect any significant splitting of HC$_3$N line. Thus, we conclude that our detection of the CCS Zeeman splitting is robust. The derived \textcolor{black}{line-of-sight} magnetic field strength is about 117 $\pm$ 21 $\mu$G, which corresponds to the normalized mass-to-magnetic flux ratio of 2.2 if we adopt the inclination angle of 45$^\circ$. Thus, we conclude that the TMC-1 filament is magnetically supercritical. Recent radiative transfer calculations of CCS and HC$_3$N lines along the line of sight suggest that the filament is collapsing with a speed of $\sim$ 0.6 km s$^{-1}$, which is comparable to three times the isothermal sound speed. This infall velocity appears to be consistent with the evolution of a gravitationally-infalling core.

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.

We conducted near-infrared (JHKs) imaging polarimetry toward the infrared dark cloud (IRDC) M17 SWex, including almost all of the IRDC filaments as well as its outskirts, with the polarimeter SIRPOL on the IRSF 1.4 m telescope. We revealed the magnetic fields of M17 SWex with our polarization-detected sources that were selected by some criteria based on their near-IR colors and the column densities toward them, which were derived from the Herschel data. The selected sources indicate not only that the ordered magnetic field is perpendicular to the cloud elongation as a whole, but also that at both ends of the elongated cloud the magnetic field appears to bent toward its central part, i.e., large-scale hourglass-shaped magnetic field perpendicular to the cloud elongation. In addition to this general trend, the elongations of the filamentary subregions within the dense parts of the cloud appear to be mostly perpendicular to their local magnetic fields, while the magnetic fields of the outskirts appear to follow the thin filaments that protrude from the dense parts. The magnetic strengths were estimated to be ~70-300 microG in the subregions, of which lengths and average number densities are ~3-9 pc and ~2-7x10^3 cm^{-3}, respectively, by the Davis-Chandrasekhar-Fermi method with the angular dispersion of our polarization data and the velocity dispersion derived from the C^{18}O (J=1-0) data obtained by the Nobeyama 45 m telescope. These field configurations and our magnetic stability analysis of the subregions imply that the magnetic field have controlled the formation/evolution of the M17 SWex cloud.

We present the results of mapping observations toward a nearby starless filamentary cloud, the Taurus Molecular Cloud 1 (TMC-1), in the CCS(JN=43-32, 45.379033 GHz) emission line, using the Nobeyama 45-m telescope. The map shows that the TMC-1 filament has a diameter of ~0.1 pc and a length of ~0.5 pc at a distance of 140 pc. The position-velocity diagrams of CCS clearly indicate the existence of velocity-coherent substructures in the filament. We identify 21 substructures that are coherent in the position-position-velocity space by eye. Most of the substructures are elongated along the major axis of the TMC-1 filament. The line densities of the subfilaments are close to the critical line density for the equilibrium (~17 Mo/pc for the excitation temperature of 10 K), suggesting that self-gravity should play an important role in the dynamics of the subfilaments.

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

The RCW 106 molecular cloud complex is an active massive star-forming region where a ministarburst is taking place. We examined its magnetic structure by near-IR polarimetric observations with the imaging polarimeter SIRPOL on the IRSF 1.4 m telescope. The global magnetic field is nearly parallel to the direction of the Galactic plane and the cloud elongation. We derived the magnetic field strength of $\sim100$-$1600~\mu$G for 71 clumps with the Davis-Chandrasekhar-Fermi method. We also evaluated the magnetic stability of these clumps and found massive star-forming clumps tend to be magnetically unstable and gravitationally unstable. Therefore, we propose a new criterion to search for massive star-forming clumps. These details suggest that the process enhancing the clump density without an increase of the magnetic flux is essential for the formation of massive stars and the necessity for accreting mass along the magnetic field lines.

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

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.

Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to map the $\mathrm{N_2H^+} \; (1-0)$ emission in two hub-filament systems with a spatial resolution of $\sim 0.034 \; \mathrm{pc}$. Kinematics are derived from sophisticated spectral fitting algorithm that accounts for line blending, large optical depth, and multiple velocity components. We identify five velocity coherent filaments and derive their velocity gradients with principal component analysis. The mass accretion rates along the filaments are up to $10^{-4} \; \mathrm{M_\odot \, \mathrm{yr^{-1}}}$ and are significant enough to affect the hub dynamics within one free-fall time ($\sim 10^5 \; \mathrm{yr}$). The $\mathrm{N_2H^+}$ filaments are in equilibrium with virial parameter $\alpha_\mathrm{vir} \sim 1.2$. We compare $\alpha_\mathrm{vir}$ measured in the $\mathrm{N_2H^+}$ filaments, $\mathrm{NH_3}$ filaments, $870 \; \mu\mathrm{m}$ dense clumps, and $3 \; \mathrm{mm}$ dense cores. The decreasing trend in $\alpha_\mathrm{vir}$ with decreasing spatial scales persists, suggesting an increasingly important role of gravity at small scales. Meanwhile, $\alpha_\mathrm{vir}$ also decreases with decreasing non-thermal motions. In combination with the absence of high-mass protostars and massive cores, our results are consistent with the global hierarchical collapse scenario.

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.

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.

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.

Linear polarization maps of the Carina Nebula were obtained at 250, 350, and 500 $\mu$m during the 2012 flight of the BLASTPol balloon-borne telescope. These measurements are combined with Planck 850 $\mu$m data in order to produce a submillimeter spectrum of the polarization fraction of the dust emission, averaged over the cloud. This spectrum is flat to within $\pm$15% (relative to the 350 $\mu$m polarization fraction). In particular, there is no evidence for a pronounced minimum of the spectrum near 350 $\mu$m, as suggested by previous ground-based measurements of other molecular clouds. This result of a flat polarization spectrum in Carina is consistent with recently-published BLASTPol measurements of the Vela C molecular cloud, and also agrees with a published model for an externally-illuminated, dense molecular cloud by Bethell and collaborators. The shape of the spectrum in Carina does not show any dependence on the radiative environment of the dust, as quantified by the Planck-derived dust temperature or dust optical depth at 353 GHz.

We have carried out observations in the 42$-$46 and 82$-$103 GHz bands with the Nobeyama 45-m radio telescope, and in the 338.2$-$339.2 and 348.45$-$349.45 GHz bands with the ASTE 10-m telescope toward three high-mass star-forming regions containing massive young stellar objects (MYSOs), G12.89+0.49, G16.86$-$2.16, and G28.28$-$0.36. We have detected HC$_{3}$N including its $^{13}$C and D isotopologues, CH$_{3}$OH, CH$_{3}$CCH, and several complex organic molecules (COMs). Combining our previous results of HC$_{5}$N in these sources, we compare the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratios in the three observed sources. The ratio in G28.28$-$0.36 is derived to be $0.091^{+0.109}_{-0.039}$, which is higher than that in G12.89+0.49 by one order of magnitude, and that in G16.86$-$2.16 by a factor of $\sim 5$. We investigate the relationship between the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratio and the $N$(CH$_{3}$CCH)/$N$(CH$_{3}$OH) ratio. The relationships of the two column density ratios in G28.28$-$0.36 and G16.86$-$2.16 are similar to each other, while HC$_{5}$N is less abundant when compared to CH$_{3}$CCH in G12.89+0.49. These results imply a chemical diversity in the lukewarm ($T \sim 20-30$ K) envelope around MYSOs. Besides, several spectral lines from complex organic molecules, including very-high-excitation energy lines, have been detected toward G12.89+0.49, while the line density is significantly low in G28.28$-$0.36. These results suggest that organic-poor MYSOs are surrounded by a carbon-chain-rich lukewarm envelope (G28.28$-$0.36), while organic-rich MYSOs, namely hot cores, are surrounded by a CH$_{3}$OH-rich lukewarm envelope (G12.89+0.49 and G16.86$-$2.16).

We have carried out interferometric observations of cyanopolyynes, HC$_{3}$N, HC$_{5}$N, and HC$_{7}$N, in the 36 GHz band toward the G28.28$-$0.36 high-mass star-forming region using the Karl G. Jansky Very Large Array (VLA) Ka-band receiver. The spatial distributions of HC$_{3}$N and HC$_{5}$N are obtained. HC$_{5}$N emission is coincident with a 450 $\mu$m dust continuum emission and this clump with a diameter of $\sim 0.2$ pc is located at the east position from the 6.7 GHz methanol maser by $\sim 0.15$ pc. HC$_{7}$N is tentatively detected toward the clump. The HC$_{3}$N : HC$_{5}$N : HC$_{7}$N column density ratios are estimated at 1.0 : $\sim 0.3$ : $\sim 0.2$ at an HC$_{7}$N peak position. We discuss possible natures of the 450 $\mu$m continuum clump associated with the cyanopolyynes. The 450 $\mu$m continuum clump seems to contain deeply embedded low- or intermediate-mass protostellar cores, and the most possible formation mechanism of the cyanopolyynes is the warm carbon chain chemistry (WCCC) mechanism. In addition, HC$_{3}$N and compact HC$_{5}$N emission is detected at the edge of the 4.5 $\mu$m emission, which possibly implies that such emission is the shock origin.

We clarify the line-of-sight structure of the Taurus Molecular Cloud 1 (TMC-1) on the basis of the CCS($J_N=4_3-3_2$) and HC$_3$N($J=5-4$) spectral data observed at a very high velocity resolution and sensitivity of $\Delta V \simeq 0.0004$ km s$^{-1}$ ($=61$ Hz) and $\Delta T_{\rm mb} \simeq 40$ mK. The data were obtained toward the cyanopolyyne peak with $\sim$30 hours integration using the Z45 receiver and the PolariS spectrometer installed in the Nobeyama 45m telescope. Analyses of the optically thin $F=4-4$ and $5-5$ hyperfine lines of the HC$_3$N emission show that the spectra consist of four distinct velocity components with a small line width ($\lesssim 0.1$ km s$^{-1}$) at $V_{\rm LSR}=$5.727, 5.901, 6.064, and 6.160 km s$^{-1}$, which we call A, B, C, and D, respectively, in the order of increasing LSR velocities. Utilizing the velocity information of the four velocity components, we further analyzed the optically thicker CCS spectrum and the other hyperfine lines of the HC$_3$N emission by solving the radiative transfer to investigate how the four velocity components overlap along the line of sight. Results indicate that they are located in the order of A, B, C, and D from far side to near side to the observer, indicating that TMC-1 is shrinking, moving inward as a whole.

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

Using Galactic Plane surveys, we have selected a massive (1200 M$_\odot$), cold (14 K) 3.6-70 $\mu$m dark IRDC G331.372-00.116. This IRDC has the potential to form high-mass stars and, given the absence of current star formation signatures, it seems to represent the earliest stages of high-mass star formation. We have mapped the whole IRDC with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.1 and 1.3 mm in dust continuum and line emission. The dust continuum reveals 22 cores distributed across the IRDC. In this work, we analyze the physical properties of the most massive core, ALMA1, which has no molecular outflows detected in the CO (2-1), SiO (5-4), and H$_2$CO (3-2) lines. This core is relatively massive ($M$ = 17.6 M$_\odot$), subvirialized (virial parameter $\alpha_{vir}=M_{vir}/M=0.14$), and is barely affected by turbulence (transonic Mach number of 1.2). Using the HCO$^+$ (3-2) line, we find the first detection of infall signatures in a relatively massive, prestellar core (ALMA1) with the potential to form a high-mass star. We estimate an infall speed of 1.54 km s$^{-1}$ and a high accretion rate of 1.96 $\times$ 10$^{-3}$ M$_\odot$ yr$^{-1}$. ALMA1 is rapidly collapsing, out of virial equilibrium, more consistent with competitive accretion scenarios rather than the turbulent core accretion model. On the other hand, ALMA1 has a mass $\sim$6 times larger than the clumps Jeans mass, being in an intermediate mass regime ($M_{J}=2.7<M\lesssim$ 30 M$_\odot$), contrary to what both the competitive accretion and turbulent core accretion theories predict.

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 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 report results of the observations of 15 regions in several molecular lines for a statistical study of massive cluster-forming clumps. We identified 24 clumps based on the C18O (J=1-0) data obtained by the NRO 45 m telescope, and found that 16 of them are associated with young clusters. The clumps associated with clusters have a typical mass, radius, and molecular density of ~1 X 10^3 Mo, ~0.5 pc, ~1 X 10^5 cm^-3, respectively. We categorized the clumps and clusters into four types according to the spatial coincidence of gas and star density, and discussed their evolutions: Clumps without clusters (Type 1), clumps showing good correlations with clusters (Type 2), clumps showing poor correlations with clusters (Type 3), and clusters with no associated clumps (Type 4). Analyses of the velocity structures and the chemical compositions imply that the clump + cluster systems should evolve from Type 1 to Type 4. We found that some of the Type 2 clumps are infalling on the clump-scale to form clusters at the clump center, which should commonly occur in the beginning of cluster formation. Interestingly, all of the identified Type 1 clumps are likely to be older than the Type 2 clumps in terms of chemical compositions, suggesting that they have been gravitationally stable for a long time possibly being supported by the strong magnetic field of > 1 mG.Type 1 clumps younger than the observed Type 2 clumps should be very rare to find because of their short lifetime.

We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for injection of turbulent momentum into high density gas at high rates.

We present calculations of molecular, atomic and ionic line emission from simulations of giant molecular cloud (GMC) collisions. We post-process snapshots of the magneto-hydrodynamical simulations presented in an earlier paper in this series by Wu et al. (2017) of colliding and non-colliding GMCs. Using photodissociation region (PDR) chemistry and radiative transfer we calculate the level populations and emission properties of $^{12}$CO $J=1-0$, [CI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ at $609\,\mu$m, [CII] $158\,\mu$m and [OI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ transition at $63\,\mu$m. From integrated intensity emission maps and position-velocity diagrams, we find that fine-structure lines, particularly the [CII] $158\,\mu$m, can be used as a diagnostic tracer for cloud-cloud collision activity. These results hold even in more evolved systems in which the collision signature in molecular lines has been diminished.

Polarized emission from aligned dust is a crucial tool for studies of magnetism in the ISM and a troublesome contaminant for studies of CMB polarization. In each case, an understanding of the significance of the polarization signal requires well-calibrated physical models of dust grains. Despite decades of progress in theory and observation, polarized dust models remain largely underconstrained. During its 2012 flight, the balloon-borne telescope BLASTPol obtained simultaneous broad-band polarimetric maps of a translucent molecular cloud at 250, 350, and 500 microns. Combining these data with polarimetry from the Planck 850 micron band, we have produced a submillimeter polarization spectrum for a cloud of this type for the first time. We find the polarization degree to be largely constant across the four bands. This result introduces a new observable with the potential to place strong empirical constraints on ISM dust polarization models in a previously inaccessible density regime. Comparing with models by Draine and Fraisse (2009), our result disfavors two of their models for which all polarization arises due only to aligned silicate grains. By creating simple models for polarized emission in a translucent cloud, we verify that extinction within the cloud should have only a small effect on the polarization spectrum shape compared to the diffuse ISM. Thus we expect the measured polarization spectrum to be a valid check on diffuse ISM dust models. The general flatness of the observed polarization spectrum suggests a challenge to models where temperature and alignment degree are strongly correlated across major dust components.

We carried out line survey observations at the 26-30 GHz band toward the four high-mass star-forming regions containing hot cores, G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Robert C. Byrd Green Bank Telescope. We have detected HC5N from all of the sources, and HC7N from the three sources, except for G10.30-0.15. We further conducted observations of HC5N at the 42-46 GHz and 82-103 GHz bands toward the three sources, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Nobeyama 45 m radio telescope. The rotational lines of HC5N with the high excitation energies (Eu/k=63-100 K), which are hardly excited in the cold dark clouds, have been detected from the three sources. The rotational temperatures of HC5N are found to be 13-20 K in the three sources. The detection of the lines with the high excitation energies and the derived rotational temperatures indicate that HC5N exists in the warm gas within 0.07-0.1 pc radii around massive young stellar objects. The column densities of HC5N in the three sources are derived to be (2.0+-2.8)*10^(13) cm-2 . We compare the ratios between N(HC5N) the column density of HC5N and W(CH3OH) the integrated intensity of the thermal CH3OH emission line among the three high-mass star-forming regions. We found a possibility of the chemical differentiation in the three high-mass star-forming regions; G28.28-0.36 shows the largest N (HC5N)/W (CH3OH) ratio of > 8.0*10^(14) in units of (K km s-1 )-1 cm-2 , while G12.89+0.49 and G16.86-2.16 show the smaller values (~2*10^(13) ).

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

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics.

We present a large-scale combination of near-infrared (near-IR) interstellar polarization data from background starlight with polarized emission data at submillimeter (sub-mm) wavelengths for the Vela C molecular cloud. The near-IR data consist of more than 6700 detections probing a range of visual extinctions between $2$ and $20\,$mag in and around the cloud. The sub-mm data was collected in Antartica by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This is the first direct combination of near-IR and sub-mm polarization data for a molecular cloud aimed at measuring the "polarization efficiency ratio" ($R_{\mathrm{eff}}$), a quantity that is expected to depend only on grain intrinsic physical properties. It is defined as $p_{500}/(p_{I}/\tau_{V})$, where $p_{500}$ and $p_{I}$ are polarization fractions at $500\,\mu$m and $I$-band, respectively, and $\tau_{V}$ is the optical depth. To ensure that the same column density of material is producing both polarization from emission and from extinction, we conducted a careful selection of near-background stars using 2MASS, $Herschel$ and $Planck$ data. This selection excludes objects contaminated by the Galactic diffuse background material as well as objects located in the foreground. Accounting for statistical and systematic uncertainties, we estimate an average $R_{\mathrm{eff}}$ value of $2.4\pm0.8$, which can be used to test the predictions of dust grain models designed for molecular clouds when such predictions become available. $R_{\mathrm{eff}}$ appears to be relatively flat as a function of the cloud depth for the range of visual extinctions probed.

We present results of wide-field $^{12}$CO ($J = 2 - 1$) and $^{13}$CO ($J = 2 - 1$) observations toward the Aquila Rift and Serpens molecular cloud complexes (25$^\circ < l < 33^\circ$ and $1^\circ < b < 6^\circ$) at an angular resolution of 3$'$.4 ($\approx$ 0.25 pc) and at a velocity resolution of 0.079 km s$^{-1}$ with the velocity coverage of $-5$ km s$^{-1} < V_{\rm LSR} <$ 35 km s$^{-1}$. We found that the $^{13}$CO emission better traces the structures seen in the extinction map and derived the $X_{\rm ^{13}CO}$-factor of this region. Applying \texttt{SCIMES} to the $^{13}$CO data cube, we identified 61 clouds and derived their masses, radii, and line widths. The line-width-radius relation of the identified clouds basically follows those of nearby molecular clouds. Majority of the identified clouds are close to virial equilibrium although the dispersion is large. By inspecting the $^{12}$CO channel maps by eye, we found several arcs which are spatially extended to 0.2 $-$ 3 degree in length. In the longitude-velocity diagrams of $^{12}$CO, we also found the two spatially-extended components which appear to converge toward Serpens South and W40 region. The existence of two components with different velocities and arcs suggests that large-scale expanding bubbles and/or flows play a role in the formation and evolution of the Serpens South and W40 cloud.

We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.

We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $\mu$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $\alpha=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse.

We report the discovery of infalling motion with rotation of S235AB the massive cluster-forming clump (~10^3 Mo) in the S235 region. Our C18O observations with the 45m telescope at the Nobeyama Radio Observatory have revealed the elliptical shape of the clump. Position-velocity (PV) diagram taken along its major axis exhibits two well-defined peaks symmetrically located with respect to the clump center, which is similar to that found for a dynamically infalling envelope with rotation around a single protostar modeled by N. Ohashi and his collaborators, indicating that the cluster-forming clump is also collapsing by the self-gravity toward the clump center. With analogue to Ohashi's model, we made a simple model of an infalling, rotating clump to fit the observed data. Based on the inferred model parameters as well as results of earlier observations and simulations in the literature, we discuss structures of the clump such as the relation among the global mass infall rate (~10^-3 Mo/yr), formation of a compact core (with a mass and size of ~4 Mo and < 0.1 pc) at the center, and a massive star (~11 Mo) forming in the core.

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at \lambda ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

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 present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

It has been shown that a realistic level of magnetization of dense molecular cloud cores can suppress the formation of a rotationally supported disk (RSD) through catastrophic magnetic braking in the axisymmetric ideal MHD limit. In this study, we present conditions for the formation of RSDs through non-ideal MHD effects computed self-consistently from an equilibrium chemical network. We find that removing from the standard MRN distribution the large population of very small grains (VSGs) of ~10 $\AA$ to few 100 $\AA$ that dominate the coupling of the bulk neutral matter to the magnetic field increases the ambipolar diffusivity by ~1--2 orders of magnitude at densities below 10$^{10}$ cm$^{-3}$. The enhanced ambipolar diffusion (AD) in the envelope reduces the amount of magnetic flux dragged by the collapse into the circumstellar disk-forming region. Therefore, magnetic braking is weakened and more angular momentum can be retained. With continuous high angular momentum inflow, RSDs of tens of AU are able to form, survive, and even grow in size, depending on other parameters including cosmic-ray ionization rate, magnetic field strength, and rotation speed. Some disks become self-gravitating and evolve into rings in our 2D (axisymmetric) simulations, which have the potential to fragment into (close) multiple systems in 3D. We conclude that disk formation in magnetized cores is highly sensitive to chemistry, especially to grain sizes. A moderate grain coagulation/growth to remove the large population of VSGs, either in the prestellar phase or during free-fall collapse, can greatly promote AD and help formation of tens of AU RSDs.

We observed the J = 9-8 and 16-15 rotational transitions of the normal species and five 13C isotopologues of HC5N to study its formation mechanisms toward the cyanopolyyne peak in Taurus Molecular Cloud-1, with the 45-m radio telescope of Nobeyama Radio Observatory. We detected the five 13C isotopologues with high signal-to-noise ratios between 12 and 20, as well as the normal species. The abundance ratios of the five 13C isotopologues of HC5N are found to be 1.00:0.97:1.03:1.05:1.16 (0.19) for [H13CCCCCN]:[HC13CCCCN]:[HCC13CCCN]:[HCCC13CCN]:[HCCCC13CN]. We do not find any significant differences among the five 13C isotopologues. The averaged [HC5N]/[13C isotopologues] abundance ratio is determined to be 94 (6), which is slightly higher than the local interstellar elemental 12C/13C ratio of 60-70. Possible formation pathways are discussed on the basis of these results.

We discuss the role of protostellar outflow feedback in clustered star formation using the observational data of recent molecular outflow surveys toward nearby cluster-forming clumps. We found that for almost all clumps, the outflow momentum injection rate is significantly larger than the turbulence dissipation rate. Therefore, the outflow feedback is likely to maintain supersonic turbulence in the clumps. For less massive clumps such as B59, L1551, and L1641N, the outflow kinetic energy is comparable to the clump gravitational energy. In such clumps, the outflow feedback probably affects significantly the clump dynamics. On the other hand, for clumps with masses larger than about 200 M$_\odot$, the outflow kinetic energy is significantly smaller than the clump gravitational energy. Since the majority of stars form in such clumps, we conclude that outflow feedback cannot destroy the whole parent clump. These characteristics of the outflow feedback support the scenario of slow star formation.

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.

We developed a dual-linear-polarization HEMT (High Electron Mobility Transistor) amplifier receiver system of the 45-GHz band (hereafter Z45), and installed it in the Nobeyama 45-m radio telescope. The receiver system is designed to conduct polarization observations by taking the cross correlation of two linearly-polarized components, from which we process full-Stokes spectroscopy. We aim to measure the magnetic field strength through the Zeeman effect of the emission line of CCS ($J_N=4_3-3_2$) toward pre-protostellar cores. A linear-polarization receiver system has a smaller contribution of instrumental polarization components to the Stokes $V$ spectra than that of the circular polarization system, so that it is easier to obtain the Stokes $V$ spectra. The receiver has an RF frequency of 42 $-$ 46 GHz and an intermediate frequency (IF) band of 4$-$8 GHz. The typical noise temperature is about 50 K, and the system noise temperature ranges from 100 K to 150K over the frequency of 42 $-$ 46 GHz. The receiver system is connected to two spectrometers, SAM45 and PolariS. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. PolariS is a newly-developed digital spectrometer with a finest frequency resolution of 60 Hz, having a capability to process the full-Stokes spectroscopy. The Half Power Beam Width (HPBW) of the beam was measured to be 37$"$ at 43 GHz. The main beam efficiency of the Gaussian main beam was derived to be 0.72 at 43 GHz. The SiO maser observations show that the beam pattern is reasonably round at about 10 \% of the peak intensity and the side-lobe level was less than 3 \% of the peak intensity. Finally, we present some examples of astronomical observations using Z45.

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

We have developed a software-based polarization spectrometer, PolariS, to acquire full-Stokes spectra with a very high spectral resolution of 61 Hz. The primary aim of PolariS is to measure the magnetic fields in dense star-forming cores by detecting the Zeeman splitting of molecular emission lines. The spectrometer consists of a commercially available digital sampler and a Linux computer. The computer is equipped with a graphics processing unit (GPU) to process FFT and cross-correlation using the CUDA (Compute Unified Device Architecture) library developed by NVIDIA. Thanks to a high degree of precision in quantization of the analog-to-digital converter and arithmetic in the GPU, PolariS offers excellent performances in linearity, dynamic range, sensitivity, bandpass flatness and stability. The software has been released under the MIT License and is available to the public. In this paper, we report the design of PolariS and its performance verified through engineering tests and commissioning observations.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was a suborbital experiment designed to map magnetic fields in order to study their role in star formation processes. BLASTPol made detailed polarization maps of a number of molecular clouds during its successful flights from Antarctica in 2010 and 2012. We present the next-generation BLASTPol instrument (BLAST-TNG) that will build off the success of the previous experiment and continue its role as a unique instrument and a test bed for new technologies. With a 16-fold increase in mapping speed, BLAST-TNG will make larger and deeper maps. Major improvements include a 2.5 m carbon fiber mirror that is 40% wider than the BLASTPol mirror and ~3000 polarization sensitive detectors. BLAST-TNG will observe in three bands at 250, 350, and 500 microns. The telescope will serve as a pathfinder project for microwave kinetic inductance detector (MKID) technology, as applied to feedhorn coupled submillimeter detector arrays. The liquid helium cooled cryostat will have a 28-day hold time and will utilize a closed-cycle $^3$He refrigerator to cool the detector arrays to 270 mK. This will enable a detailed mapping of more targets with higher polarization resolution than any other submillimeter experiment to date. BLAST-TNG will also be the first balloon-borne telescope to offer shared risk observing time to the community. This paper outlines the motivation for the project and the instrumental design.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment designed to study the role magnetic fields play in star formation. BLASTPol has had two science flights from McMurdo Station, Antarctica in 2010 and 2012. These flights have produced thousands of polarization vectors at 250, 350 and 500 microns in several molecular cloud targets. We present the design, specifications, and progress towards the next-generation BLASTPol experiment (BLAST-TNG). BLAST-TNG will fly a 40% larger diameter primary mirror, with almost 8 times the number of polarization-sensitive detectors resulting in a factor of 16 increase in mapping speed. With a spatial resolution of 22 arcseconds and four times the field of view of BLASTPol, BLAST-TNG will bridge the angular scales between Planck's low resolution all-sky maps and ALMA's ultra-high resolution narrow fields. The new receiver has a larger cryogenics volume, allowing for a 28 day hold time. BLAST-TNG employs three arrays of Microwave Kinetic Inductance Detectors (MKIDs) with 30% fractional bandwidth at 250, 350 and 500 microns. In this paper, we will present the new BLAST-TNG instrument and science objectives.

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.

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.

Protostellar outflows have been shown theoretically to be capable of maintaining supersonic turbulence in cluster-forming clumps and keeping the star formation rate per free-fall time as low as a few percent. We aim to test two basic predictions of this outflow-regulated cluster formation model, namely (1) the clump should be close to virial equilibrium and (2) the turbulence dissipation rate should be balanced by the outflow momentum injection rate, using recent outflow surveys toward 8 nearby cluster-forming clumps (B59, L1551, L1641N, Serpens Main Cloud, Serpens South, $\rho$ Oph, IC 348, and NGC 1333). We find, for almost all sources, that the clumps are close to virial equilibrium and the outflow momentum injection rate exceeds the turbulence momentum dissipation rate. In addition, the outflow kinetic energy is significantly smaller than the clump gravitational energy for intermediate and massive clumps with $M_{\rm cl} \gtrsim {\rm a \ few} \times 10^2 M_\odot$, suggesting that the outflow feedback is not enough to disperse the clump as a whole. The number of observed protostars also indicates that the star formation rate per free-fall time is as small as a few percent for all clumps. These observationally-based results strengthen the case for outflow-regulated cluster formation.

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

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

We present combined Submillimeter-Array (SMA) + single-dish images of the (sub)millimeter dust continuum emission toward two prestellar cores SM1 and B2-N5 in the nearest star cluster forming region, $\rho$ Ophiuchus. Our combined images indicate that SM1 and B2-N5 consist of three and four condensations, respectively, with masses of $10^{-2}-10^{-1}M_\odot$ and sizes of a few hundred AU. The individual condensations have mean densities of $10^8-10^9$ cm$^{-3}$ and the masses are comparable to or larger than the critical Bonner-Ebert mass, indicating that the self-gravity plays an important role in the dynamical evolution of the condensations. The coalescence timescale of these condensations is estimated to be about $10^4$ yr, which is comparable to the local gravitational collapse timescale, suggesting that merging of the condensations, instead of accretion, plays an essential role in the star formation process. These results challenge the standard theory of star formation, where a single, rather featureless prestellar core collapses to form at most a couple of condensations, each of which potentially evolves into a protostar that is surrounded by a rotating disk where planets are created.

Stars form predominantly in clusters inside dense clumps of turbulent, magnetized molecular clouds. The typical size and mass of the cluster-forming clumps are \sim 1 pc and \sim 10^2 - 10^3 M_\odot, respectively. Here, we discuss some recent progress on theoretical and observational studies of clustered star formation in such parsec-scale clumps with emphasis on the role of protostellar outflow feedback. Recent simulations indicate that protostellar outflow feedback can maintain supersonic turbulence in a cluster-forming clump, and the clump can keep a virial equilibrium long after the initial turbulence has decayed away. In the clumps, star formation proceeds relatively slowly; it continues for at least several global free-fall times of the parent dense clump (t_{ff}\sim a few x 10^5 yr). The most massive star in the clump is formed at the bottom of the clump gravitational potential well at later times through the filamentary mass accretion streams that are broken up by the outflows from low-mass cluster members. Observations of molecular outflows in nearby cluster-forming clumps appear to support the outflow-regulated cluster formation model.

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.

Near-infrared (IR) imaging polarimetry in the J, H, and Ks bands has been carried out for the protostellar cluster region around NGC 2264 IRS 2 in the Monoceros OB1 molecular cloud. Various infrared reflection nebulae clusters (IRNCs) associated with NGC 2264 IRS 2 and IRAS 12 S1 core were detected as well as local infrared reflection nebulae (IRNe). The illuminating sources of the IRNe were identified with known or new near- and mid-IR sources. In addition, 314 point-like sources were detected in all three bands and their aperture polarimetry was studied. Using a color-color diagram, reddened field stars and diskless pre-main sequence stars were selected to trace the magnetic field (MF) structure of the molecular cloud. The mean polarization position angle of the point-like sources is 81 \pm 29 degree in the cluster core, and 58 \pm 24 degree in the perimeter of the cluster core, which is interpreted as the projected direction on the sky of the MF in the observed region of the cloud. The Chandrasekhar-Fermi method gives a rough estimate of the MF strength to be about 100 {\mu}G. A comparison with recent numerical simulations of the cluster formation implies that the cloud dynamics is controlled by the relatively strong MF. The local MF direction is well associated with that of CO outflow for IRAS 12 S1 and consistent with that inferred from submillimeter polarimetry. In contrast, the local MF direction runs roughly perpendicular to the Galactic MF direction.

We investigate the physical properties of dense cores formed in turbulent, magnetized, parsec-scale clumps of molecular clouds, using three-dimensional numerical simulations that include protostellar outflow feedback. The dense cores are identified in the simulated density data cube through a clumpfind algorithm. We find that the core velocity dispersion does not show any clear dependence on the core size, in contrast to Larson's linewidth-size relation, but consistent with recent observations. In the absence of a magnetic field, the majority of the cores have supersonic velocity dispersions. A moderately-strong magnetic field reduces the dispersion to a subsonic or at most transonic value typically. Most of the cores are out of virial equilibrium, with the external pressure dominating the self-gravity. The implication is that the core evolution is largely controlled by the outflow-driven turbulence. Even an initially-weak magnetic field can retard star formation significantly, because the field is amplified by the outflow-driven turbulence to an equipartition strength, with the distorted field component dominating the uniform one. In contrast, for a moderately-strong field, the uniform component remains dominant. Such a difference in the magnetic structure is evident in our simulated polarization maps of dust thermal emission; it provides a handle on the field strength. Recent polarization measurements show that the field lines in cluster-forming clumps are spatially well-ordered. It is indicative of a moderately-strong, dynamically important, field which, in combination with outflow feedback, can keep the rate of star formation in embedded clusters at the observationally-inferred, relatively-slow rate of several percent per free-fall time.

Stars form in dense cores of magnetized molecular clouds. If the magnetic flux threading the cores is dragged into the stars, the stellar field would be orders of magnitude stronger than observed. This well-known "magnetic flux problem" demands that most of the core magnetic flux be decoupled from the matter that enters the star. We carry out the first exploration of what happens to the decoupled magnetic flux in 3D, using an MHD version of the ENZO adaptive mesh refinement code. The field-matter decoupling is achieved through a sink particle treatment, which is needed to follow the protostellar accretion phase of star formation. We find that the accumulation of the decoupled flux near the accreting protostar leads to a magnetic pressure buildup. The high pressure is released anisotropically, along the path of least resistance. It drives a low-density expanding region in which the decoupled magnetic flux is expelled. This decoupling-enabled magnetic structure has never been seen before in 3D MHD simulations of star formation. It generates a strong asymmetry in the protostellar accretion flow, potentially giving a kick to the star. In the presence of an initial core rotation, the structure presents an obstacle to the formation of a rotationally supported disk, in addition to magnetic braking, by acting as a rigid magnetic wall that prevents the rotating gas from completing a full orbit around the central object. We conclude that the decoupled magnetic flux from the stellar matter can strongly affect the protostellar collapse dynamics.

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.

The Serpens South embedded cluster, which is located at the constricted part in a long filamentary infrared dark cloud, is believed to be in very early stage of cluster formation. We present results of near-infrared (JHKs) polarization observations toward the filamentary cloud. Our polarization measurements of near-infrared point sources indicate a well-ordered global magnetic field that is perpendicular to the main filament, implying that the magnetic field is likely to have controlled the formation of the main filament. On the other hand, the sub-filaments, which converge on the central part of the cluster, tend to run along the magnetic field. The global magnetic field appears to be curved in the southern part of the main filament. Such morphology is consistent with the idea that the global magnetic field is distorted by gravitational contraction along the main filament toward the northern part that contains larger mass. Applying the Chandrasekhar-Fermi method, the magnetic field strength is roughly estimated to be a few x 100 microgauss, suggesting that the filamentary cloud is close to magnetically critical as a whole.

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.

Stars form predominantly in clusters inside dense clumps of molecular clouds that are both turbulent and magnetized. The typical size and mass of the cluster-forming clumps are $\sim 1$ pc and $\sim 10^2 - $ 10$^3$ M$_\odot$, respectively. Here, we discuss some recent progress on numerical simulations of clustered star formation in such parsec-scale dense clumps with emphasis on the role of magnetic fields. The simulations have shown that magnetic fields tend to slow down global gravitational collapse and thus star formation, especially in the presence of protostellar outflow feedback. Even a relatively weak can retard star formation significantly, because the field is amplified by supersonic turbulence to an equipartition strength. However, in such a case, the distorted field component dominates the uniform one. In contrast, if the field is moderately strong, the uniform component remains dominant. Such a difference in the magnetic structure is observed in simulated polarization maps of dust thermal emission. Recent polarization measurements show that the field lines in nearby cluster-forming clumps are spatially well-ordered, indicative of a rather strong field. In such strongly-magnetized clumps, star formation should proceed relatively slowly; it continues for at least several global free-fall times of the parent dense clump ($t_{\rm ff}\sim $ a few $\times 10^5$ yr).

Magnetic fields are generally expected to increase the characteristic mass of stars formed in stellar clusters, because they tend to increase the effective Jeans mass. We test this expectation using adaptive mesh refinement (AMR) magnetohydrodynamic simulations of cluster formation in turbulent magnetized clumps of molecular clouds, treating stars as accreting sink particles. We find that, contrary to the common expectation, a magnetic field of strength in the observed range decreases, rather than increases, the characteristic stellar mass. It (1) reduces the number of intermediate-mass stars that are formed through direct turbulent compression, because sub-regions of the clump with masses comparable to those of stars are typically magnetically subcritical and cannot be compressed directly into collapse, and (2) increases the number of low-mass stars that are produced from the fragmentation of dense filaments. The filaments result from mass accumulation along the field lines. In order to become magnetically supercritical and fragment, the filament must accumulate a large enough column density (proportional to the field strength), which yields a high volume density (and thus a small thermal Jeans mass) that is conducive to forming low-mass stars. We find, in addition, that the characteristic stellar mass is reduced further by outflow feedback. The conclusion is that both magnetic fields and outflow feedback are important in shaping the stellar initial mass function (IMF).

We made deep NIR imaging polarimetry toward the Serpens cloud core. The polarization vector maps enable us to newly detect 24 small IR reflection nebulae with YSOs. Polarization measurements of NIR point sources indicate an hourglass-shaped magnetic field, of which symmetry axis is nearly perpendicular to the elongation of the C18O (J=1-0) or submillimeter continuum emission. The bright part of C18O (J=1-0), submillimeter continuum cores as well as many class 0/I objects are located just toward the constriction region of the hourglass-shaped magnetic field. Applying the CF method, the magnetic field strength was estimated to be ~100 muG, suggesting that the ambient region of the Serpens cloud core is moderately magnetically supercritical. These suggest that the Serpens cloud core first contracted along the magnetic field to be an elongated cloud, which is perpendicular to the magnetic field, and that then the central part contracted cross the magnetic field due to the high density in the central region of the cloud core, where star formation is actively continuing. Comparison of this magnetic field with the previous observations of molecular gas and large-scale outflows suggests a possibility that the cloud dynamics is controlled by the magnetic field, protostellar outflows and gravitational inflows. This appears to be in good agreement with the outflow-driven turbulence model and implies the importance of the magnetic field to continuous star formation in the center region of the cluster forming region.

Using the archive data of the H13CO+ (J=1-0) line emission taken with the Nobeyama 45 m radio telescope with a spatial resolution of about 0.01pc, we have identified 68 dense cores in the central dense region of the rho Ophiuchi main cloud. The H13CO+ data also indicates that the fractional abundance of H13CO+ relative to H2 is roughly inversely proportional to the square root of the H2 column density with a mean of 1.72 x 10^{-11}. The mean radius, FWHM line width, and LTE mass of the identified cores are estimated to be 0.045 +- 0.011 pc, 0.49 +- 0.14 km/s, and 3.4 +- 3.6 Msolar, respectively. The majority of the identified cores have subsonic internal motions. The virial ratio, the ratio of the virial mass to the LTE mass, tends to decrease with increasing the LTE mass and about 60 percent of the cores have virial ratios smaller than 2, indicating that these cores are not transient structures but self-gravitating. The detailed virial analysis suggests that the surface pressure often dominates over the self-gravity and thus plays a crucial role in regulating core formation and evolution. By comparing the rho Oph cores with those in the Orion A molecular cloud observed with the same telescope, we found that the statistical properties of the core physical quantities are similar between the two clouds if the effect of the different spatial resolutions is corrected. The line widths of the rho Oph cores appear to be nearly independent of the core radii over the range of 0.01 - 0.1 pc and deviate upwards from the Heyer & Brunt relation. This may be evidence that turbulent motions are driven by protostellar outflows in the cluster environment.

(Abridged) We investigate massive star formation in turbulent, magnetized, parsec-scale clumps of molecular clouds including protostellar outflow feedback using Enzo-based MHD simulations with accreting sink particles and effective resolution $2048^3$. We find that, in the absence of regulation by magnetic fields and outflow feedback, massive stars form readily in a turbulent, moderately condensed clump of $\sim 1,600$ solar masses, along with a cluster of hundreds of lower mass stars. The massive stars are fed at high rates by (1) transient dense filaments produced by large-scale turbulent compression at early times, and (2) by the clump-wide global collapse resulting from turbulence decay at late times. In both cases, the bulk of the massive star's mass is supplied from outside a 0.1 pc-sized "core" that surrounds the star. In our simulation, the massive star is clump-fed rather than core-fed. The need for large-scale feeding makes the massive star formation prone to regulation by outflow feedback, which directly opposes the feeding processes. The outflows reduce the mass accretion rates onto the massive stars by breaking up the dense filaments that feed the massive star formation at early times, and by collectively slowing down the global collapse that fuel the massive star formation at late times. The latter is aided by a moderate magnetic field of strength in the observed range. We conclude that the massive star formation in our simulated turbulent, magnetized, parsec-scale clump is outflow-regulated and clump-fed (ORCF for short). An important implication is that the formation of low-mass stars in a dense clump can affect the formation of massive stars in the same clump, through their outflow feedback on the clump dynamics.

We present the results of C18O observations by the Nobeyama Millimeter Array toward dense clumps with radii of ~ 0.3 pc in six cluster-forming regions including massive (proto)stars. We identified 171 cores, whose radius, line width, and molecular mass range from 0.01 to 0.09 pc, 0.43 to 3.33 km/s, and 0.5 to 54.1 Mo, respectively. Many cores with various line widths exist in one clump, and the index of the line width-radius relationship of the cores and the parental clump differs from core to core in the clump. This indicates that the degree of dissipation of the turbulent motion varies for each core in one clump. Although the mass of the cores increases with the line width, most cores are gravitationally bound by the external pressure. In addition, the line width and the external pressure of the cores tend to decrease with the distance from the center of the clump, and these dependencies may be caused by the inner H2 density structure of the clump that affects the physical properties of the cores. Moreover, the number density of the cores and the number density of young (proto)stars have a similar relationship to the average H2 density of the clumps. Thus, our findings suggest that the cluster is formed in the clump through the formation of such multiple cores, whose physical properties would have been strongly related to the H2 density structure and the turbulent motion of the clump.

We carry out three-dimensional MHD simulations of star formation in turbulent, magnetized clouds, including ambipolar diffusion and feedback from protostellar outflows. The calculations focus on relatively diffuse clouds threaded by a strong magnetic field capable of resisting severe tangling by turbulent motions and retarding global gravitational contraction in the cross-field direction. They are motivated by observations of the Taurus molecular cloud complex (and, to a lesser extent, Pipe Nebula), which shows an ordered large-scale magnetic field, as well as elongated condensations that are generally perpendicular to the large-scale field. We find that stars form in earnest in such clouds when enough material has settled gravitationally along the field lines that the mass-to-flux ratios of the condensations approach the critical value. Only a small fraction (of order 1% or less) of the nearly magnetically-critical, condensed material is turned into stars per local free-fall time, however. The slow star formation takes place in condensations that are moderately supersonic; it is regulated primarily by magnetic fields, rather than turbulence. The quiescent condensations are surrounded by diffuse halos that are much more turbulent, as observed in the Taurus complex. Strong support for magnetic regulation of star formation in this complex comes from the extremely slow conversion of the already condensed, relatively quiescent C$^{18}$O gas into stars, at a rate two orders of magnitude below the maximum, free-fall value. We analyze the properties of dense cores, including their mass spectrum, which resembles the stellar initial mass function.

The majority of stars are thought to form in clusters. Cluster formation in dense clumps of molecular clouds is strongly influenced, perhaps controlled, by supersonic turbulence. We have previously shown that the turbulence in regions of active cluster formation is quickly transformed by the forming stars through protostellar outflows, and that the outflow-driven protostellar turbulence is the environment in which most of the cluster members form. Here, we take initial steps in quantifying the global properties of the protostellar turbulence through 3D MHD simulations. We find that collimated outflows are more efficient in driving turbulence than spherical outflows that carry the same amounts of momentum. Gravity plays an important role in shaping the turbulence, generating infall motions in the cluster forming region that balance the outward motions driven by outflows. The resulting quasi-equilibrium state is maintained through a slow rate of star formation, with a fraction of the total mass converted into stars per free fall time as low as a few percent. Magnetic fields are dynamically important even in magnetically supercritical clumps, provided that their initial strengths are not far below the critical value for static cloud support. We find that the mass weighted PDF of the volume density is often, although not always, approximately lognormal. The PDFs of the column density deviate more strongly from lognormal distributions. There is a prominent break in the velocity power spectrum of the protostellar turbulence, which may provide a way to distinguish it from other types of turbulence.

Most, perhaps all, stars go through a phase of vigorous outflow during formation. We examine, through 3D MHD simulation, the effects of protostellar outflows on cluster formation. We find that the initial turbulence in the cluster-forming region is quickly replaced by motions generated by outflows. The protostellar outflow-driven turbulence (``protostellar turbulence'' for short) can keep the region close to a virial equilibrium long after the initial turbulence has decayed away. We argue that there exist two types of turbulence in star-forming clouds: a primordial (or ``interstellar'') turbulence and a protostellar turbulence, with the former transformed into the latter mostly in embedded clusters such as NGC 1333. Since the majority of stars are thought to form in clusters, an implication is that the stellar initial mass function is determined to a large extent by the stars themselves, through outflows which individually limit the mass accretion onto forming stars and collectively shape the environments (density structure and velocity field) in which most cluster members form. We speculate that massive cluster-forming clumps supported by protostellar turbulence gradually evolve towards a highly centrally condensed ``pivotal'' state, culminating in rapid formation of massive stars in the densest part through accretion.

The effect of smooth cloud boundaries on the interaction of steady planar shock waves with interstellar clouds is studied using a high-resolution local AMR technique with a second-order accurate axisymmetric Godunov hydrodynamic scheme. A 3D calculation is also done to confirm the results of the 2D ones. We consider an initially spherical cloud whose density distribution is flat near the cloud center and has a power-law profile in the cloud envelope. When an incident shock is transmitted into a smooth cloud, velocity gradients in the cloud envelope steepen the smooth density profile at the upstream side, resulting in a sharp density jump having an arc-like shape. Such a ``slip surface'' forms immediately when a shock strikes a cloud with a sharp boundary. For smoother boundaries, the formation of slip surface and therefore the onset of hydrodynamic instabilities are delayed. Since the slip surface is subject to the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the shocked cloud is eventually destroyed in $\sim 3-10$ cloud crushing times. After complete cloud destruction, small blobs formed by fragmentation due to hydrodynamic instabilities have significant velocity dispersions of the order of 0.1 $v_b$, where $v_b$ is the shock velocity in the ambient medium. This suggests that turbulent motions generated by shock-cloud interaction are directly associated with cloud destruction. The interaction of a shock with a cold HI cloud should lead to the production of a spray of small HI shreds, which could be related to the small cold clouds recently observed by Stanimirovic & Heiles (2005). The linewidth-size relation obtained from our 3D simulation is found to be time-dependent. A possibility for gravitational instability triggered by shock compression is also discussed.

The efficiency of star formation, defined as the ratio of the stellar to total (gas and stellar) mass, is observed to vary from a few percent in regions of dispersed star formation to about a third in cluster-forming cores. This difference may reflect the relative importance of magnetic fields and turbulence in controlling star formation. We investigate the interplay between supersonic turbulence and magnetic fields using numerical simulations, in a sheet-like geometry. We demonstrate that star formation with an efficiency of a few percent can occur over several gravitational collapse times in moderately magnetically subcritical clouds that are supersonically turbulent. The turbulence accelerates star formation by reducing the time for dense core formation. The dense cores produced are predominantly quiescent, with subsonic internal motions. These cores tend to be moderately supercritical. They have lifetimes long compared with their local gravitational collapse time. Some of the cores collapse to form stars, while others disperse away without star formation. In turbulent clouds that are marginally magnetically supercritical, the star formation efficiency is higher, but can still be consistent with the values inferred for nearby embedded clusters. If not regulated by magnetic fields at all, star formation in a multi-Jeans mass cloud endowed with a strong initial turbulence proceeds rapidly, with the majority of cloud mass converted into stars in a gravitational collapse time. The efficiency is formally higher than the values inferred for nearby cluster-forming cores, indicating that magnetic fields are dynamically important even for cluster formation.

The evolution of a gas shell, swept by the supernova remnant of a massive first generation star, is studied with H_2 and HD chemistry taken into account. When a first-generation star explodes as a supernova, H_2 and HD molecules are formed in the swept gas shell and effectively cool the gas shell to temperatures of 32 K - 154 K. If the supernova remnant can sweep to gather the ambient gas, the gas shell comes to be dominated by its self-gravity, and hence, is expected to fragment. Our result shows that for a reasonable range of temperatures (200 K - 1000 K) of interstellar gas, the formation of second-generation stars can be triggered by a single supernova or hypernova.

We investigate numerically the combined effects of supersonic turbulence, strong magnetic fields and ambipolar diffusion on cloud evolution leading to star formation. We find that, in clouds that are initially magnetically subcritical, supersonic turbulence can speed up star formation, through enhanced ambipolar diffusion in shocks. The speedup overcomes a major objection to the standard scenario of low-mass star formation involving ambipolar diffusion, since the diffusion time scale at the average density of a molecular cloud is typically longer than the cloud life time. At the same time, the strong magnetic field can prevent the large-scale supersonic turbulence from converting most of the cloud mass into stars in one (short) turbulence crossing time, and thus alleviate the high efficiency problem associated with the turbulence-controlled picture for low-mass star formation. We propose that relatively rapid but inefficient star formation results from supersonic collisions of somewhat subcritical gas in strongly magnetized, turbulent clouds. The salient features of this shock-accelerated, ambipolar diffusion-regulated scenario are demonstrated with numerical experiments.

We study the non-axisymmetric evolution of magnetized clouds, using a 2D MHD code based on the physically motivated thin-disk approximation. We found that such clouds become unstable to non-axisymmetric perturbations after the supercritical cores are formed due to ambipolar diffusion. We show that for a wide range of initial cloud parameters, the $m=2$ mode grows nonlinearly into a bar during the isothermal collapse after the supercritical core formation. The supercritical core can break up into fragments during or after the isothermal phase of cloud evolution. The outcome of fragmentation depends on the initial cloud conditions, such as the magnetic field strength, rotation rate, amount of cloud mass and mass distribution. It is classified into three different types: (1) separate core formation, (2) bar fragmentation, and (3) disk fragmentation. These three types of fragmentation loosely correspond to the empirical classification of embedded binary and multiple systems of Looney, Mundy, & Welch, based on millimeter dust continuum observations. The well-studied starless core, L1544, appears to belong to the bar fragmentation type. We expect it to produce a highly elongated, opaque bar at the center in the future, which should break up into fragments of initial masses in the substellar regime.

Isolated low-mass stars are formed in dense cores of molecular clouds. In the standard picture, the cores are envisioned to condense out of strongly magnetized clouds through ambipolar diffusion. Most previous calculations based on this scenario are limited to axisymmetric cloud evolution leading to a single core, which collapses to form an isolated star or stellar system at the center. These calculations are here extended to the nonaxisymmetric case under thin-disk approximation, which allows for a detailed investigation into the process of fragmentation, fundamental to binary, multiple system, and cluster formation. We have shown previously that initially axisymmetric, magnetically subcritical clouds with an $m=2$ density perturbation of modest fractional amplitude ($\sim 5%$) can develop highly elongated bars, which facilitate binary and multiple system formation. In this paper, we show that in the presence of higher order ($m\ge 3$) perturbations of similar amplitude such clouds are capable of breaking up into a set of discrete dense cores. These multiple cores are magnetically supercritical. They are expected to collapse into single stars or stellar systems individually and, collectively, to form a small stellar group. Our calculations demonstrate that the standard scenario for single star formation involving magnetically subcritical clouds and ambipolar diffusion can readily produce more than one star, provided that the cloud mass is well above the Jeans limit and relatively uniformly distributed. The fragments develop in the central part of the cloud, after the region has become magnetically supercritical but before rapid collapse sets in. It is enhanced by the flattening of mass distribution along the field lines and by the magnetic tension force.

In the standard scenario of isolated low-mass star formation, strongly magnetized molecular clouds are envisioned to condense gradually into cores, driven by ambipolar diffusion. Once the cores become magnetically supercritical, they collapse to form stars. Most previous studies based on this scenario are limited to axisymmetric calculations leading to single supercritical core formation. The assumption of axisymmetry has precluded a detailed investigation of cloud fragmentation, generally thought to be a necessary step in the formation of binary and multiple stars. In this contribution, we describe the non-axisymmetric evolution of initially magnetically subcritical clouds using a newly-developed MHD code. It is shown that non-axisymmetric perturbations of modest fractional amplitude ($\sim 5%$) can grow nonlinearly in such clouds during the supercritical phase of cloud evolution, leading to the production of either a highly elongated bar or a set of multiple dense cores.

In the context of star formation through fragmentation of an extremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary gas clouds is explored with 1D numerical hydrodynamics coupled with non-equilibrium chemistry of H$_2$ and HD. It is found that the cloud evolution is governed mainly by the initial central density ($n_{\rm c,0}$) and H$_2$ abundance ($x_{\rm H_2,0}$). In particular, the evolution of low-density filaments bifurcates at a threshold H$_2$ abundance of $x_{\rm H_2,cr}\simeq 3 \times 10^{-3}$, beyond which HD cooling overwhelms H$_2$ cooling. The contraction of a low density filament with $x_{\rm H_2, 0}\gtrsim x_{\rm H_2,cr}$ is strongly decelerated when $n_{\rm c}$ reaches a critical density of HD, and the filament is expected to fragment at $\sim 10^{7}$ cm$^{-3}$. The fragment mass is lowered to be $\approx 10M_\odot$. In contrast, the contraction of a low density filament with $x_{\rm H_2, 0}\lesssim x_{\rm H_2,cr}$ is regulated by H$_2$ cooling. In this case, the fragment mass is as high as $\approx 10^2M_\odot$. For a high-density filament, the cloud evolution is governed by H$_2$ cooling. The fragmentation is not expected to take place until the cloud becomes opaque in H$_2$ lines at $n_{\rm c,0}\sim 10^{12-13}$ cm$^{-3}$, so that the fragment mass is reduced to 1-2 M$_\odot$. As a result, the stellar IMF could be bimodal and deficient in sub-solar mass stars, where the high mass peak is around $10M_\odot$ or $10^2M_\odot$, dependently on $n_{\rm c,0}$ and $x_{\rm H_2,0}$.

We have begun a systematic numerical study of the nonlinear growth of nonaxisymmetric perturbations during the ambipolar diffusion-driven evolution of initially magnetically subcritical molecular clouds, with an eye on the formation of binaries, multiple stellar systems and small clusters. In this initial study, we focus on the $m=2$ (or bar) mode, which is shown to be unstable during the dynamic collapse phase of cloud evolution after the central region has become magnetically supercritical. We find that, despite the presence of a strong magnetic field, the bar can grow fast enough that for a modest initial perturbation (at 5% level) a large aspect ratio is obtained during the isothermal phase of cloud collapse. The highly elongated bar is expected to fragment into small pieces during the subsequent adiabatic phase. Our calculations suggest that the strong magnetic fields observed in some star-forming clouds and envisioned in the standard picture of single star formation do not necessarily suppress bar growth and fragmentation; on the contrary, they may actually promote these processes, by allowing the clouds to have more than one (thermal) Jeans mass to begin with without collapsing promptly. Nonlinear growth of the bar mode in a direction perpendicular to the magnetic field, coupled with flattening along field lines, leads to the formation of supercritical cores that are triaxial in general. It removes a longstanding objection to the standard scenario of isolated star formation involving subcritical magnetic field and ambipolar diffusion based on the likely prolate shape inferred for dense cores. Continuted growth of the bar mode in already elongated starless cores, such as L1544, may lead to future binary and multiple star formation.

In the context of the star formation through the fragmentation of an extremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary gas clouds is explored with H$_2$ and HD chemistry. It is found by 1D hydrodynamical simulations that the cloud evolution is prescribed mainly by the initial density ($n_0$) and H$_2$ abundance ($x_{\rm H_2,0}$). In particular, it turns out that the evolution of low-density filaments ($n_0 \lesssim 10^5$ cm$^{-3}$) bifurcates at a critical H$_2$ abundance of $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$, beyond which HD cooling overwhelms H$_2$ cooling. The numerical results indicate that the stellar IMF is likely to be double-peaked and deficient in sub-solar mass stars, where the high mass peak of the IMF is around $10M_\odot$ or $10^2M_\odot$, dependently on the initial density and H$_2$ abundance. If the gas in protogalactic clouds is photoionized by UV radiation or shock-heated, the H$_2$ abundance could exceed $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$ by H$^-$ reactions. Then, the high mass peak would be $O(10) M_\odot$.

The collapse and fragmentation of filamentary primordial gas clouds are explored using 1D and 2D hydrodynamical simulations coupled with the nonequilibrium processes of H2 formation. The simulations show that depending upon the initial density,there are two occasions for the fragmentation of primordial filaments. If a filament has relatively low initial density, the radial contraction is slow due to less effective H2 cooling. This filament tends to fragment into dense clumps before the central density reaches $10^{8-9}$ cm$^{-3}$, where H2 cooling by three-body reactions is effective and the fragment mass is more massive than some tens $M_\odot$. In contrast, if a filament is initially dense, the more effective H2 cooling with the help of three-body reactions allows the filament to contract up to $n\sim 10^{12}$ cm$^{-3}$. After the density reaches $n\sim 10^{12}$ cm$^{-3}$, the filament becomes optically thick to H2 lines and the radial contraction subsequently almost stops. At this final hydrostatic stage, the fragment mass is lowered down to $\approx 1M_\odot$ because of the high density of the filament. The dependence of the fragment mass upon the initial density could be translated into the dependence on the local amplitude of random Gaussian density fields or the epoch of the collapse of a parent cloud. Hence, it is predicted that the initial mass function of Population III stars is likely to be bimodal with peaks of $\approx 10^2 M_\odot$ and $\approx 1M_\odot$, where the relative heights could be a function of the collapse epoch.

We examined the gravitational contraction of isothermal molecular cloud cores with slow rotation by means of two-dimensional numerical simulations. Applying a sink-cell method, we followed the evolution of the cloud cores up to the stages at which most of the matter accretes onto the central region (i.e., a protostar and a rotationally-supported circumstellar disk). We show that both an infalling disklike envelope and a rotationally-supported disk around the central star are natural outcome of the gravitational contraction of a {\it prolate} cloud core with slow rotation. The early evolution of the infalling envelopes resembles sheet models recently proposed by Hartmann and coworkers. In the infalling disklike envelope, the radial profiles of the density, radial velocity, and azimuthal velocity can be approximated by $\rho\propto r^{-1.5}$, $v_r\propto r^{-0.5}$, and $v_\phi\propto r^{-1}$, respectively. The fate of the infalling envelopes is also discussed.

Performing 1D hydrodynamical calculations coupled with non-equilibrium processes for H2 formation, we pursue the thermal and dynamical evolution of filamentary primordial clouds and attempt to make an estimate on the mass of population III stars. It is found that, almost independent of initial conditions, a filamentary cloud continues to collapse nearly isothermally due to H_2 cooling until the cloud becomes optically thick against the H_2 lines. During the collapse the cloud structure separates into two parts, i.e., a denser spindle and a diffuse envelope. The spindle contracts quasi-statically, and thus the line mass of the spindle keeps a characteristic value determined solely by the temperature ($\sim 800$ K). Applying a linear theory, we find that the spindle is unstable against fragmentation during the collapse. The wavelength of the fastest growing perturbation lessens as the collapse proceeds. Consequently, successive fragmentation could occur. When the central density exceeds $n_c \sim 10^{10-11} cm^{-3}$, the successive fragmentation may cease since the cloud becomes opaque against the H_2 lines and the collapse decelerates appreciably. The mass of the first star is then expected to be typically $\sim 3 M_\odot$, which may grow up to $\sim 16 M_\odot$ by accreting the diffuse envelope. Thus, the first-generation stars are anticipated to be massive but not supermassive.

The fraction of star-forming galaxies in rich clusters of galaxies increases rapidly with the redshift. This is interpreted as the result of a rapid evolution of cluster galaxies, though its mechanism is not yet clear. One hypothesis is that if galaxies run into the dense intracluster medium (ICM) regions, starbursts can be induced due to compression of molecular clouds in the galaxies by a raised "external" ICM pressure. In a merging cluster, there would be dense ICM regions and some galaxies may experience a rapid increase of external pressure. Thus, if the mechanism works, starbursts would occur in such galaxies. Such a situation is considered to be realized in Coma cluster, and galaxies showing a recent starburst have a characteristic spatial distribution; they populate an elongated region between two sub-clusters. We examined the above hypothesis by measuring the spatial distribution of blue galaxies, regarded as star-forming galaxies in Abell~168 to check whether a case like Coma is realized in this recently merging cluster. However, we could not find an enhanced blue-galaxy distribution between two sub-clusters. We discuss the causes of our result, including possibilities of surviving the hypothesis.

Recent X-ray observations have established that collisions between subclusters of galaxies are rather common phenomena. Prompted by such observations, we have performed N-body simulations of two equal-mass subclusters of galaxies, which are going to merge. We first have confirmed that only a part of kinetic energy associated with the relative motion of two subclusters is converted to internal energy of each subcluster at the first encounter. We performed simulations for a variety of dark-matter (DM) distributions, and find that the time scale for washing out the double peak structures after the first encounter strongly depends on the distribution of the density and velocity distributions of the DM. It takes longer when the DM is spatially extended. According to our calculation it takes more than $4 imes 10^9$yr after the first encounter until the density contour shows only a single peak. To explain the high fraction of the clusters with substructures among nearby clusters, Richstone,Loeb Turner(1992; hereafter RLT) required recent($$ $10^9 h^{-1}$yr) cluster formation. However, our results show that the timescale for subcluster merging is still uncertain and possibly much longer than the time scale assumed by RLT. Caution should be exercised when concluding that the density of the universe is high by using RLT's method.

The dynamical instability of a self-gravitating magnetized filamentary cloud was investigated while taking account of rotation around its axis. The filamentary cloud of our model is supported against self-gravity in part by both a magnetic field and rotation. The density distribution in equilibrium was assumed to be a function of the radial distance from the axis, $ \rho _0 (r) \, = \, \rho _{\rm c} \, ( 1 \, + \, r ^2 / 8 H ^2 ) ^{-2} $, where $ \rho _{\rm c} $ and $ H $ are model parameters specifying the density on the axis and the length scale, respectively; the magnetic filed was assumed to have both longitudinal ($z$-) and azimuthal ($\varphi$-) components with a strength of $ B _0 (r) \, \propto \, \sqrt{ \rho _0 (r) } $. The rotation velocity was assumed to be $ v _{0\varphi} \, = \, \Omega _{\rm c} \, r \, (1 \, + \, r ^2 / 8 H ^2 ) ^{-1/2} $. We obtained the growth rate and eigenfunction numerically for (1) axisymmetric $ ( m \, = \, 0 ) $ perturbations imposed on a rotating cloud with a longitudinal magnetic field, (2) non-axisymmetric $ ( m \, = \, 1 ) $ perturbations imposed on a rotating cloud with a longitudinal magnetic filed, and (3) axisymmetric perturbations imposed on a rotating cloud with a helical magnetic field. The fastest growing perturbation is an axisymmetric one for all of the model clouds studied. Its wavelength is $ \lambda _{\rm max} \, \le \, 11.14 \, H $ for a non-rotating cloud without a magnetic field, and is shorter