Source author record

Fumitaka Nakamura

Fumitaka Nakamura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Cluster Formation in GGD12-15: Infall Motion with Rotation of the Natal Clump

We report results of observations of the GGD12-15 region, where cluster formation is actively taking place, using various molecular emission lines. The C18O (J= 1-0) emission line reveals a massive clump in the region with a mass of ~2800 Mo distributed over ~2 pc. The distribution of the C18O(J= 3-2) emission is similar to that of a star cluster forming therein, with an elliptical shape of ~1 pc in size. A bipolar molecular outflow driven by IRS 9Mc, a constituent star of the cluster, is blowing in a direction perpendicular to the elongated C18O (J= 3-2) distribution, covering the entire clump. There is a massive core with a radius of 0.3 pc and a mass of 530 Mo in the center of the clump. There are two velocity components around the core, which are prominent in a position-velocity (PV) diagram along the major axis of the clump. In addition, a PV diagram along the minor axis of the clump, which is parallel to the outflow, shows a velocity gradient opposite to that of the outflow. The velocity structure strongly indicates the infalling motion of the clump. Comparison of the observational data with a simple model of infalling oblate clumps indicates that the clump is undergoing gravitational contraction with rotation.

preprint2022arXiv

Vibrationally-excited Lines of HC$_{3}$N Associated with the Molecular Disk around the G24.78+0.08 A1 Hyper-compact H$_{\rm {II}}$ Region

We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH$_{3}$CN vibrationally-excited line ($v_{8}=1$), which indicates that the HC$_{3}$N emission is tracing the disk around the G24 HC H$_{\rm {II}}$ region previously identified by the CH$_{3}$CN lines. We derive the $^{13}$CH$_{3}$CN/HC$^{13}$CCN abundance ratios around G24 and compare them to the CH$_{3}$CN/HC$_{3}$N abundance ratios in disks around Herbig Ae and T Tauri stars. The $^{13}$CH$_{3}$CN/HC$^{13}$CCN ratios around G24 ($\sim 3.0-3.5$) are higher than the CH$_{3}$CN/HC$_{3}$N ratios in the other disks ($\sim 0.03-0.11$) by more than one order of magnitude. The higher CH$_{3}$CN/HC$_{3}$N ratios around G24 suggest that the thermal desorption of CH$_{3}$CN in the hot dense gas and efficient destruction of HC$_{3}$N in the region irradiated by the strong UV radiation are occurring. Our results indicate that the vibrationally-excited HC$_{3}$N lines can be used as a disk tracer of massive protostars at the HC H$_{\rm {II}}$ region stage, and the combination of these nitrile species will provide information of not only chemistry but also physical conditions of the disk structures.

preprint2021arXiv

Magnetic Fields in Massive Star-Forming Regions (MagMaR) II. Tomography Through Dust and Molecular Line Polarization in NGC 6334I(N)

Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to the field. By considering previous findings, the field retains a pinched shape which can be traced to clump length-scales from the envelope scales traced by ALMA, suggesting that the field is dynamically important across multiple length-scales in this region. The CS total intensity emission is found to be optically thick ($τ_{\mathrm{CS}} = 32 \pm 12$) while the C$^{33}$S emission appears to be optically thin ($τ_{\mathrm{C^{33}S}} = 0.1 \pm 0.01$). This suggests that sources of anisotropy other than large velocity gradients, i.e. anisotropies in the radiation field are required to explain the polarized emission from CS seen by ALMA. By using four variants of the Davis-Chandrasekhar-Fermi technique and the angle dispersion function methods (ADF), we obtain an average of estimates for the magnetic field strength onto the plane of the sky of $\left< \mathrm{B}_{\mathrm{pos}} \right> = 16$ mG from the dust and $\left< \mathrm{B}_{\mathrm{pos}} \right> \sim 2$ mG from the CS emission, where each emission traces different molecular hydrogen number densities. This effectively enables a tomographic view of the magnetic field within a single ALMA observation.

preprint2021arXiv

The ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES). V. Deuterated Molecules in the 70 $μ$m dark IRDC G14.492-00.139

We have observed the 70 $μ$m dark infrared dark cloud (IRDC) G14.492-00.139 in the N$_2$D$^+$ $J$=3--2, DCO$^+$ $J$=3--2, DCN $J$=3--2, and C$^{18}$O $J$=2--1 lines, using the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES). We find that the spatial distribution is different among the observed emission from the deuterated molecular lines. The N$_2$D$^+$ emission traces relatively quiescent regions, while both the DCO$^+$ and DCN emission emanates mainly from regions with signs of active star-formation. In addition, the DCO$^+$/N$_2$D$^+$ ratio is found to be lower in several dense cores than in starless cores embedded in low-mass star-forming regions. By comparing the observational results with chemical model calculations, we discuss the origin of the low DCO$^+$/N$_2$D$^+$ ratio in this IRDC clump. The low DCO$^+$/N$_2$D$^+$ ratio can be explained if the temperature of the dense cores is in the range between the sublimation temperature of N$_2$ ($\sim$20 K) and CO ($\sim$25 K). The results suggest that the dense cores in G14.492-00.139 are warmer and denser than the dense cores in low-mass star-forming regions.

preprint2020arXiv

A detailed analysis on the cloud structure and dynamics in Aquila Rift

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.

preprint2020arXiv

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios, in the MYSOs are derived to be $\sim 15$. This value is lower than that in a low-mass warm carbon chain chemistry (WCCC) source by more than one order of magnitude. We compare the observed CCH/HC$_{5}$N ratios with hot-core model calculations (Taniguchi et al. 2019). The observed ratios in the MYSOs can be best reproduced by models when the gas temperature is $\sim 85$ K, which is higher than in L1527, a low-mass WCCC source ($\sim 35$ K). These results suggest that carbon-chain molecules detected around the MYSOs exist at least partially in higher temperature regions than those in low-mass WCCC sources. There is no significant difference in column density among the three $^{13}$C isotopologues of HC$_{3}$N in G12.89+0.49 and G16.86-2.16, while HCC$^{13}$CN is more abundant than the others in G28.28--0.36. We discuss carbon-chain chemistry around the three MYSOs based on the CCH/HC$_{5}$N ratio and the $^{13}$C isotopic fractionation of HC$_{3}$N.

preprint2020arXiv

Cloud structures in M17 SWex : Possible cloud-cloud collision

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.

preprint2020arXiv

From diffuse gas to dense molecular cloud cores

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.

preprint2020arXiv

GMC Collisions as Triggers of Star Formation. VII. The Effect of Magnetic Field Strength on Star Formation

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.

preprint2020arXiv

Large-scale molecular gas distribution in the M17 cloud complex: dense gas conditions of massive star formation?

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.

preprint2020arXiv

Star Cluster Formation in Orion A

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 $Σ$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$μ$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.

preprint2020arXiv

The CARMA-NRO Orion Survey: Filament Formation via Collision-Induced Magnetic Reconnection -- The Stick in Orion A

A unique filament is identified in the {\it Herschel} maps of the Orion A giant molecular cloud. The filament, which, we name the Stick, is ruler-straight and at an early evolutionary stage. Transverse position-velocity diagrams show two velocity components closing in on the Stick. The filament shows consecutive rings/forks in C$^{18}$O(1-0) channel maps, which is reminiscent of structures generated by magnetic reconnection. We propose that the Stick formed via collision-induced magnetic reconnection (CMR). We use the magnetohydrodynamics (MHD) code Athena++ to simulate the collision between two diffuse molecular clumps, each carrying an anti-parallel magnetic field. The clump collision produces a narrow, straight, dense filament with a factor of $>$200 increase in density. The production of the dense gas is seven times faster than free-fall collapse. The dense filament shows ring/fork-like structures in radiative transfer maps. Cores in the filament are confined by surface magnetic pressure. CMR can be an important dense-gas-producing mechanism in the Galaxy and beyond.

preprint2020arXiv

The CARMA-NRO Orion Survey: Protostellar Outflows, Energetics, and Filamentary Alignment

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.

preprint2020arXiv

The molecular cloud lifecycle

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.

preprint2019arXiv

First clear detection of the CCS Zeeman splitting toward the pre-stellar core, Taurus Molecular Cloud-1

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

preprint2019arXiv

Interaction between Northern Coal Sack in the Cyg OB 7 cloud complex and the multiple super nova remnants including HB 21

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.

preprint2019arXiv

Investigation of Chemical Differentiation among the NGC2264 Cluster-Forming Clumps

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.

preprint2019arXiv

Nobeyama 45-m mapping observations toward nearby molecular clouds, Orion A, Aquila Rift, and M17: Project overview

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.

preprint2016arXiv

Balloon-Borne Submillimeter Polarimetry of the Vela C Molecular Cloud: Systematic Dependence of Polarization Fraction on Column Density and Local Polarization-Angle Dispersion

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 μm. In this initial paper, we show our 500 μ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.

preprint2016arXiv

Dense core properties in the Infrared Dark cloud G14.225-0.506 revealed by ALMA

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 $μ$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 $α=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.

preprint2016arXiv

Discovery of Infalling Motion with Rotation of the Cluster-Forming Clump S235AB and Its Implication to the Clump Structures

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.

preprint2016arXiv

Protostellar Disk Formation Enabled by Removal of Small Dust Grains

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 $Å$ to few 100 $Å$ 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.

preprint2016arXiv

Submillimeter Polarization Spectrum in the Vela C Molecular Cloud

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

preprint2016arXiv

The Deuterium Fraction in Massive Starless Cores and Dynamical Implications

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, $α_{\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.

preprint2016arXiv

The intrinsic abundance ratio and X-factor of CO isotopologues in L1551 shielded from FUV photodissociation

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.

preprint2015arXiv

Dense Clumps and Candidates for Molecular Outflows in W40

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.

preprint2015arXiv

Implication of Formation Mechanisms of HC5N in TMC-1 as Studied by 13C Isotopic Fractionation

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.

preprint2015arXiv

Observations of Protostellar Outflow Feedback in Clustered Star Formation

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.

preprint2015arXiv

Z45: A New 45-GHz Band Dual-Polarization HEMT Receiver for the NRO 45-m Radio Telescope

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.

preprint2014arXiv

Cluster Formation Triggered by Filament Collisions in Serpens South

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.

preprint2014arXiv

Confronting Outflow-Regulated Cluster Formation Model with Observations

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, $ρ$ 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.

preprint2014arXiv

High abundance ratio of $^{13}$CO to C$^{18}$O toward photon-dominated regions in the Orion-A giant molecular cloud

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 $<$ $τ_{\rm ^{13}CO}$ $<$ 1.54 and 0.01 $<$ $τ_{\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.

preprint2014arXiv

Software Polarization Spectrometer "PolariS"

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.

preprint2014arXiv

Star Cluster Formation and Feedback

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.

preprint2014arXiv

The Next Generation BLAST Experiment

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.

preprint2014arXiv

The next-generation BLASTPol experiment

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.

preprint2013arXiv

The Dynamical State of The Serpens South Filamentary Infrared Dark Cloud

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.

preprint2012arXiv

Present-Day Star Formation: Protostellar Outflows and Clustered Star Formation

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.

preprint2012arXiv

Substellar-Mass Condensations in Prestellar Cores

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, $ρ$ 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.

preprint2011arXiv

Clustered Star Formation in Magnetic Clouds: Properties of Dense Cores Formed in Outflow-Driven Turbulence

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.

preprint2011arXiv

Complex Scattered Radiation Fields and Multiple Magnetic Fields in the Protostellar Cluster in NGC 2264

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

preprint2011arXiv

Evidence For Cloud-Cloud Collision and Parsec-Scale Stellar Feedback Within the L1641-N Region

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.

preprint2011arXiv

Magnetic Flux Expulsion in Star Formation

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.

preprint2011arXiv

Molecular Outflows From the Protocluster, Serpens South

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.

preprint2010arXiv

Lowering the Characteristic Mass of Cluster Stars by Magnetic Fields and Outflow Feedback

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

preprint2010arXiv

Near-Infrared Imaging Polarimetry of the Serpens Cloud Core: Magnetic Field Structure, Outflows, and Inflows in A Cluster Forming Clump

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.

preprint2010arXiv

Physical Properties of Dense Cores in the Rho Ophiuchi Main Cloud and A Significant Role of External Pressures in Clustered Star Formation

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.

preprint2010arXiv

The Molecular Outflows in the rho Ophiuchi Main Cloud: Implications For Turbulence Generation

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.

preprint2010arXiv

Theory of Cluster Formation: Effects of Magnetic Fields

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

preprint2009arXiv

Outflow Feedback Regulated Massive Star Formation in Parsec-Scale Cluster Forming Clumps

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

preprint2007arXiv

Protostellar Turbulence Driven by Collimated Outflows

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.

preprint2001arXiv

The Stellar IMF in Very Metal-Deficient Gas

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