Source author record

Yasuo Fukui

Yasuo Fukui 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

86works
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

86 published item(s)

preprint2022arXiv

A Multiwavelength Study of the Sgr B Region: Contiguous Cloud-Cloud Collisions Triggering Widespread Star Formation Events?

The Sgr\,B region, including Sgr\,B1 and Sgr\,B2, is one of the most active star-forming regions in the Galaxy. Hasegawa et al. (1994) originally proposed that Sgr\,B2 was formed by a cloud-cloud collision (CCC) between two clouds with velocities of $\sim$45 km~s$^{-1}$ and $\sim$75 km~s$^{-1}$. However, some recent observational studies conflict with this scenario. We have re-analyzed this region, by using recent, fully sampled, dense-gas data and by employing a recently developed CCC identification methodology, with which we have successfully identified more than 50 CCCs and compared them at various wavelengths. We found two velocity components that are widely spread across this region and that show clear signatures of a CCC, each with a mass of $\sim$10$^6$ $M_\odot$. Based on these observational results, we suggest an alternative scenario, in which contiguous collisions between two velocity features with a relative velocity of $\sim$20 km~s$^{-1}$ created both Sgr\,B1 and Sgr\,B2. The physical parameters, such as the column density and the relative velocity of the colliding clouds, satisfy a relation that has been found to apply to the most massive Galactic CCCs, meaning that the triggering of high-mass star formation in the Galaxy and starbursts in external galaxies can be understood as being due to the same physical CCC process.

preprint2022arXiv

An ALMA study of the massive molecular clump N159W-North in the Large Magellanic Cloud: A possible gas flow penetrating one of the most massive protocluster systems in the Local Group

Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation of populous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angular resolution of $\sim$0."25($\sim$0.07 pc) revealed more than five protostellar sources with CO outflows within the main ridge clump. One of the thermal continuum sources, MMS-2, shows especially massive/dense nature whose total H$_2$ mass and peak column density are $\sim$10$^{4}$ $M_{\odot}$ and $\sim$10$^{24}$ cm$^{-2}$, respectively, and harbors massive ($\sim$100 $M_{\odot}$) starless core candidates identified as its internal substructures. The main ridge containing this source can be categorized as one of the most massive protocluster systems in the Local Group. The CO high-resolution observations found several distinct filamentary clouds extending southward from the star-forming spots. The CO (1-0) data set with a larger field of view reveals a conical-shaped, $\sim$30 pc long complex extending toward the northern direction. These features indicate that a large-scale gas compression event may have produced the massive star-forming complex. Based on the striking similarity between the N159W-North complex and the previously reported other two high-mass star-forming clouds in the nearby regions, we propose a $"$teardrops inflow model$"$ that explains the synchronized, extreme star formation across $>$50 pc, including one of the most massive protocluster clumps in the Local Group.

preprint2022arXiv

Charge Exchange X-ray Emission Detected in Multiple Shells of Supernova Remnant G296.1-0.5

Recent high-resolution X-ray spectroscopy revealed possible presence of charge exchange (CX) X-ray emission in supernova remnants (SNRs). Although CX is expected to take place at outermost edges of SNR shells, no significant measurement has been reported so far due to the lack of nearby SNR samples. Here we present an X-ray study of SNR G296.1$-$0.5, which has a complicated multiple-shell structure, with the Reflection Grating Spectrometer (RGS) onboard XMM-Newton. We select two shells in different regions and find that in both regions O VII line shows a high forbidden-to-resonance ($f/r$) ratio that cannot be reproduced by a simple thermal model. Our spectral analysis suggests a presence of CX and the result is also supported by our new radio observation, where we discover evidence of molecular clouds associated with these shells. Assuming G296.1$-$0.5 has a spherical shock, we estimate that CX is dominant in a thin layer with a thickness of 0.2--0.3\% of the shock radius. The result is consistent with a previous theoretical expectation and we therefore conclude that CX occurs in G296.1$-$0.5.

preprint2022arXiv

The First Detection of a Protostellar CO Outflow in the Small Magellanic Cloud with ALMA

Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$.

preprint2022arXiv

Three-dimensional geometry and dust/gas ratios in massive star forming regions over the entire LMC as revealed by IRSF/SIRIUS survey

We derive the entire dust extinction (Av) map for the Large Magellanic Cloud (LMC) estimated from the color excess at near-infrared wavelengths. Using the percentile method we recently adopted to evaluate Av distribution along the line of sight, we derive the three-dimensional(3D)Av maps of the three massive star forming regions of N44, N79 and N11 based on the IRSF/SIRIUS point source catalog. The 3D AV maps are compared with the hydrogen column densities N(H) of three different velocity components where one is of the LMC disk velocity and the other two are of velocities lower than the disk velocity. As a result, we obtain 3D dust geometry suggesting that gas collision is on-going between the different velocity components. We also find difference in the timing of the gas collision between the massive star forming regions, which indicates that the gas collision in N44, N79 and N11 occurred later than that in 30 Doradus. In addition, difference by a factor of two in Av/N(H) is found between the velocity components for N44, while significant difference is not found for N79 and N11. From the 3D geometry and Av/N(H) in each star forming region, we suggest that the massive star formation in N44 was induced by an external trigger of tidal interaction between the LMC and the SMC, while that in N79 and N11 is likely to have been induced by internal triggers such as gas converging from the galactic spiral arm and expansion of a supershell although the possibility of tidal interaction cannot be ruled out.

preprint2020arXiv

A low-velocity bipolar outflow from a deeply embedded object in Taurus revealed by the Atacama Compact Array

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the Atacama Large Millimeter/submillimeter Array$-$Atacama Compact Array (ALMA$-$ACA) stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitters in our sample. All of the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

preprint2020arXiv

ALMA CO Observations of a Giant Molecular Cloud in M33: Evidence for High-Mass Star Formation Triggered by Cloud-Cloud Collisions

We report the first evidence for high-mass star formation triggered by collisions of molecular clouds in M33. Using the Atacama Large Millimeter/submillimeter Array, we spatially resolved filamentary structures of giant molecular cloud 37 in M33 using $^{12}$CO($J$ = 2-1), $^{13}$CO($J$ = 2-1), and C$^{18}$O($J$ = 2-1) line emission at a spatial resolution of $\sim$2 pc. There are two individual molecular clouds with a systematic velocity difference of $\sim$6 km s$^{-1}$. Three continuum sources representing up to $\sim$10 high-mass stars with the spectral types of B0V-O7.5V are embedded within the densest parts of molecular clouds bright in the C$^{18}$O($J$ = 2-1) line emission. The two molecular clouds show a complementary spatial distribution with a spatial displacement of $\sim$6.2 pc, and show a V-shaped structure in the position-velocity diagram. These observational features traced by CO and its isotopes are consistent with those in high-mass star-forming regions created by cloud-cloud collisions in the Galactic and Magellanic Cloud HII regions. Our new finding in M33 indicates that the cloud-cloud collision is a promising process to trigger high-mass star formation in the Local Group.

preprint2020arXiv

ALMA Observations of Giant Molecular Clouds in M33 I: Resolving Star Formation Activities in the Giant Molecular Filaments Possibly Formed by a Spiral Shock

We report molecular line and continuum observations toward one of the most massive giant molecular clouds (GMCs), GMC-16, in M33 using ALMA with an angular resolution of 0$''$44 $\times$ 0$''$27 ($\sim$2 pc $\times$ 1 pc). We have found that the GMC is composed of several filamentary structures in $^{12}$CO and $^{13}$CO ($J$ = 2-1). The typical length, width, and total mass are $\sim$50-70 pc, $\sim$5-6 pc, and $\sim$10$^{5}$ $M_{\odot}$, respectively, which are consistent with those of giant molecular filaments (GMFs) as seen in the Galactic GMCs. The elongations of the GMFs are roughly perpendicular to the direction of the galaxy's rotation, and several H$\;${\sc ii} regions are located at the downstream side relative to the filaments with an offset of $\sim$10-20 pc. These observational results indicate that the GMFs are considered to be produced by a galactic spiral shock. The 1.3 mm continuum and C$^{18}$O ($J$ = 2-1) observations detected a dense clump with the size of $\sim$2 pc at the intersection of several filamentary clouds, which is referred to as the $"$hub filament,$"$ possibly formed by a cloud-cloud collision. A strong candidate for protostellar outflow in M33 has also been identified at the center of the clump. We have successfully resolved the parsec-scale local star formation activity in which the galactic scale kinematics may induce the formation of the parental filamentary clouds.

preprint2020arXiv

An ALMA view of the Galactic super star cluster RCW38 at 270-AU resolution

We report millimeter/submillimeter continuum and molecular line observations of the Galactic super star cluster RCW 38, obtained from the Atacama Large Millimeter/Submillimeter Array with a minimum angular resolution of $0''.17\times0''.15$ ($\simeq289\,{\rm AU}\times255\,{\rm AU}$). The C$^{18}$O image reveal many massive condensations embedded within filamentary structures extending along the northwest-southeast direction in the center of cluster. The condensations have sizes of 0.01-0.02 pc, H$_2$ column densities of $10^{23}$-$10^{24}$ cm$^{-2}$, and H$_2$ masses of 10-130 $M_\odot$. In addition, the 233-GHz continuum image reveals two dense, small millimeter-sources with radii of 460 and 200 AU (Source A and Source B). Source A is embedded within the most massive C$^{18}$O condensation, whereas no counterpart is seen for Source B. The masses of Source A and Source B are estimated as 13 and 3 $M_\odot$ at the optically-thin limit, respectively. The C$^{18}$O emission shows a velocity gradient of 2 km s$^{-1}$ at the central 2000 AU of Source A, which could be interpreted as a Keplerian rotation with a central mass of a few $M_\odot$ or infall motion of gas. Further, the ALMA $^{12}$CO data reveal that Source A and Source B are associated with molecular outflows exhibiting maximum velocities of $\sim$30-70 km s$^{-1}$. The outflows have short dynamical timescales of $<$1000 yr and high mass outflow rates of $\sim10^{-4}$-$10^{-3}$ $M_\odot$ yr$^{-1}$. These observational signatures suggest an early evolutionary phase of the massive star formation in Source A and Source B.

preprint2020arXiv

FOREST unbiased Galactic plane imaging survey with the Nobeyama 45 m telescope (FUGIN): Possible evidence of cloud-cloud collisions triggering high-mass star formation in the giant molecular cloud M16 (Eagle Nebula)

M16, the Eagle Nebula, is an outstanding \HII \ region which exhibits extensive high-mass star formation and hosts remarkable "pillars". We herein obtained new $^{12}$CO $J=$1-0 data for the region observed with NANTEN2, which were combined with the $^{12}$CO $J=$1-0 data obtained using FUGIN survey. These observations revealed that a giant molecular cloud (GMC) of $\sim 1.3 \times 10^5$ \Msun \ is associated with M16, which is elongated by over 30 pc and is perpendicular to the galactic plane, at a distance of 1.8 kpc. This GMC can be divided into the northern (N) cloud, the eastern (E) filament, the southeast (SE) cloud, the southeast (SE) filament, and the southern (S) cloud. We also found two velocity components (blue and red shifted component) in the N cloud. The blue-shifted component shows a ring-like structure, as well as the red-shifted component coincides with the intensity depression of the ring-like structure. The position-velocity diagram of the components showed a V-shaped velocity feature. The spatial and velocity structures of the cloud indicated that two different velocity components collided with each other at a relative velocity of 11.6 \kms. The timescale of the collision was estimated to be $\sim 4 \times 10^5$ yr. The collision event reasonably explains the formation of the O9V star ALS15348, as well as the shape of the Spitzer bubble N19. A similar velocity structure was found in the SE cloud, which is associated with the O7.5V star HD168504. In addition, the complementary distributions of the two velocity components found in the entire GMC suggested that the collision event occurred globally. On the basis of the above results, we herein propose a hypothesis that the collision between the two components occurred sequentially over the last several $10^{6}$ yr and triggered the formation of O-type stars in the NGC6611 cluster.

preprint2020arXiv

FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45 m telescope (FUGIN). VI. Dense gas and mini-starbursts in the W43 giant molecular cloud complex

We performed new large-scale $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=$1--0 observations of the W43 giant molecular cloud complex in the tangential direction of the Scutum arm ($l\sim {30^\circ}$) as a part of the FUGIN project. The low-density gas traced by $^{12}$CO is distributed over 150 pc $\times$ 100 pc ($l \times b$), and has a large velocity dispersion (20-30 km s$^{-1}$). However, the dense gas traced by C$^{18}$O is localized in the W43 Main, G30.5, and W43 South (G29.96-0.02) high-mass star-forming regions in the W43 GMC complex, which have clumpy structures. We found at least two clouds with a velocity difference of $\sim$ 10-20 km s$^{-1}$, both of which are likely to be physically associated with these high-mass star-forming regions based on the results of high $^{13}$CO $J=$ 3-2 to $J =$ 1-0 intensity ratio and morphological correspondence with the infrared dust emission. The velocity separation of these clouds in W43 Main, G30.5, and W43 South is too large for each cloud to be gravitationally bound. We also revealed that the dense gas in the W43 GMC has a high local column density, while "the current SFE" of entire the GMC is low ($\sim 4\%$) compared with the W51 and M17 GMC. We argue that the supersonic cloud-cloud collision hypothesis can explain the origin of the local mini-starbursts and dense gas formation in the W43 GMC complex.

preprint2020arXiv

Formation of Massive Star Clusters by Fast HI Gas Collision

Young massive clusters (YMCs) are dense aggregates of young stars, which are essential to galaxy evolution, owing to their ultraviolet radiation, stellar winds, and supernovae. The typical mass and radius of YMCs are M~10^4 M_sun and R~1 pc, respectively, indicating that many stars are located in a small region. The formation of YMC precursor clouds may be difficult because a very compact massive cloud should be formed before stellar feedback blows off the cloud. Recent observational studies suggest that YMCs can be formed as a consequence of the fast HI gas collision with a velocity of ~100 km s^-1, which is the typical velocity of the galaxy-galaxy interaction. In this study, we examine whether the fast HI gas collision triggers YMC formation using three-dimensional magnetohydrodynamics simulations, which include the effects of self-gravity, radiative cooling/heating, and chemistry. We demonstrate that massive gravitationally bound gas clumps with M >10^4 M_sun and L~4 pc are formed in the shock compressed region induced by the fast HI gas collision, which massive gas clumps can evolve into YMCs. Our results show that the YMC precursors are formed by the global gravitational collapse of molecular clouds, and YMCs can be formed even in low-metal environments, such as the Magellanic Clouds. Additionally, the very massive YMC precursor cloud, with M >10^5 M_sun, can be created when we consider the fast collision of HI clouds, which may explain the origin of the very massive stellar cluster R136 system in the Large Magellanic Cloud.

preprint2020arXiv

FRagmentation and Evolution of dense cores Judged by ALMA (FREJA). I (Overview). Inner $\sim$1000 au structures of prestellar/protostellar cores in Taurus

We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the Atacama Large Millimeter/submillimeter Array-Atacama Compact Array (ALMA-ACA) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small-scale fragmentation/coalescence processes occur in a region smaller than 0.1 pc, which may determine the final core mass associated with individual protostar formation before starting the dynamical collapse of the core with central density of $\sim$(0.3-1) $\times$ 10$^6$ cm$^{-3}$.

preprint2020arXiv

High-mass star formation in Orion B triggered by cloud-cloud collision: Merging molecular clouds in NGC 2024

We performed new comprehensive $^{13}$CO($J$=2--1) observations toward NGC 2024, the most active star forming region in Orion B, with an angular resolution of $\sim$100'' obtained with NANTEN2. We found that the associated cloud consists of two independent velocity components. The components are physically connected to the H{\sc ii} region as evidenced by their close correlation with the dark lanes and the emission nebulosity. The two components show complementary distribution with a displacement of $\sim$0.6 pc. Such complementary distribution is typical to colliding clouds discovered in regions of high-mass star formation. We hypothesize that a cloud-cloud collision between the two components triggered the formation of the late O-type stars and early B stars localized within 0.3 pc of the cloud peak. The duration time of the collision is estimated to be 0.3 million years from a ratio of the displacement and the relative velocity $\sim$3 km s$^{-1}$ corrected for probable projection. The high column density of the colliding cloud $\sim$10$^{23}$ cm$^{-2}$ is similar to those in the other high-mass star clusters in RCW 38, Westerlund 2, NGC 3603, and M42, which are likely formed under trigger by cloud-cloud collision. The present results provide an additional piece of evidence favorable to high-mass star formation by a major cloud-cloud collision in Orion.

preprint2020arXiv

Massive star formation in the Carina nebula complex and Gum 31 -- I. The Carina nebula complex

Herein, we present results from observations of the 12CO (J=1-0), 13CO (J=1-0), and 12CO (J=2-1) emission lines toward the Carina nebula complex (CNC) obtained with the Mopra and NANTEN2 telescopes. We focused on massive-star-forming regions associated with the CNC including the three star clusters Tr14, Tr15, and Tr16, and the isolated WR-star HD92740. We found that the molecular clouds in the CNC are separated into mainly four clouds at velocities -27, -20, -14, and -8 km/s. Their masses are 0.7x10^4Msun, 5.0x10^4 Msun, 1.6x10^4 Msun, and 0.7x10^4 Msun, respectively. Most are likely associated with the star clusters, because of their high 12CO (J=2-1)/12CO (J=1-0) intensity ratios and their correspondence to the Spitzer 8 micron distributions. In addition, these clouds show the observational signatures of cloud--cloud collisions. In particular, there is a V-shaped structure in the position--velocity diagram and a complementary spatial distribution between the -20 km/s cloud and the -14 km/s cloud. Based on these observational signatures, we propose a scenario wherein the formation of massive stars in the clusters was triggered by a collision between the two clouds. By using the path length of the collision and the assumed velocity separation, we estimate the timescale of the collision to be ~1 Myr. This is comparable to the ages of the clusters estimated in previous studies.

preprint2020arXiv

Molecular and Atomic Clouds toward the Wolf-Rayet Nebula NGC 2359: Possible Evidence for Isolated High-Mass Star Formation Triggered by a Cloud-Cloud Collision

NGC 2359 is an HII region located in the outer Galaxy that contains the isolated Wolf-Rayet (WR) star HD 56925. We present millimeter/submillimeter observations of $^{12}$CO($J$ = 1-0, 3-2) line emission toward the entire nebula. We identified that there are three molecular clouds at VLSR $\sim$37, $\sim$54, and $\sim$67 km s$^{-1}$, and three HI clouds: two of them are at VLSR $\sim$54 km s$^{-1}$ and the other is at $\sim$63 km s$^{-1}$. These clouds except for the CO cloud at 67 km s$^{-1}$ are limb-brightened in the radio continuum, suggesting part of each cloud has been ionized. We newly found an expanding gas motion of CO/HI, whose center and expansion velocities are $\sim$51 and $\sim$4.5 km s$^{-1}$, respectively. This is consistent with large line widths of the CO and HI clouds at 54 km s$^{-1}$. The kinematic temperature of CO clouds at 37 and 54 km s$^{-1}$ are derived to be 17 and 61 K, respectively, whereas that of the CO cloud at 67 km s$^{-1}$ is only 6 K, indicating that the former two clouds have been heated by strong UV radiation. We concluded that the 37 and 54 km s$^{-1}$ CO clouds and three HI clouds are associated with NGC 2359, even if these clouds have different velocities. Although the velocity difference including the expanding motion are typical signatures of the stellar feedback from the exciting star, our analysis revealed that the observed large momentum for the 37 km s$^{-1}$ CO cloud cannot be explained only by the total wind momentum of the WR star and its progenitor. We therefore propose an alternative scenario that the isolated high-mass progenitor of HD 56925 was formed by a collision between the CO clouds at 37 and 54 km s$^{-1}$. If we apply the collision scenario, NGC 2359 corresponds to the final phase of the cloud-cloud collision.

preprint2020arXiv

Observations of the [CI] ($^3P_1$-$^3P_0$) emission toward the massive star-forming region RCW38: further evidence for highly-clumped density distribution of the molecular gas

We present observations of the $^3P_1$-$^3P_0$ fine-structure line of atomic carbon using the ASTE 10 m sub-mm telescope towards RCW38, the youngest super star cluster in the Milky Way. The detected [CI] emission is compared with the CO $J$ = 1-0 image cube presented in Fukui et al. (2016) which has an angular resolution of 40$^{\prime \prime}$ ($\sim$ 0.33 pc). The overall distribution of the [CI] emission in this cluster is similar to that of the $^{13}$CO emission. The optical depth of the [CI] emission was found to be $τ$ = 0.1-0.6, suggesting mostly optically thin emission. An empirical conversion factor from the [CI] integrated intensity to the H$_2$ column density was estimated as $X_{\rm [CI]}$ = 6.3 $\times$ 10$^{20}$ cm$^{-2}$ K$^{-1}$ km$^{-1}$ s (for visual extinction: $A_V$ $\le$ 10 mag) and 1.4 $\times$ 10$^{21}$ cm$^{-2}$ K$^{-1}$ km$^{-1}$ s (for $A_V$ of 10-100 mag). The column density ratio of the [CI] to CO ($N_{\rm [CI]}/N_{\rm CO}$) was derived as $\sim$ 0.1 for $A_V$ of 10-100 mag, which is consistent with that of the Orion cloud presented in Ikeda et al. (2002). However, our results cover an $A_V$ regime of up to 100 mag, which is wider than the coverage found in Orion, which reach up to $\sim$ 60 mag. Such a high [CI]/CO ratio in a high $A_V$ region is difficult to be explained by the plane-parallel photodissociation region (PDR) model, which predicts that this ratio is close to 0 due to the heavy shielding of the ultraviolet (UV) radiation. Our results suggest that the molecular gas in this cluster is highly clumpy, allowing deep penetration of UV radiation even at averaged $A_V$ values of 100 mag. Recent theoretical works have presented models consistent with such clumped gas distribution with a sub-pc clump size (e.g., Tachihara et al. 2018).

preprint2020arXiv

Rapid and efficient mass collection by a supersonic cloud-cloud collision as a major mechanism of high-mass star formation

A supersonic cloud-cloud collision produces a shock-compressed layer which leads to formation of high-mass stars via gravitational instability. We carried out a detailed analysis of the layer by using the numerical simulations of magneto-hydrodynamics which deal with colliding molecular flows at a relative velocity of 20 km s$^{-1}$ (Inoue & Fukui 2013). Maximum density in the layer increases from 1000 cm$^{-3}$ to more than $10^{5}$ cm$^{-3}$ within 0.3 Myrs by compression, and the turbulence and the magnetic field in the layer are amplified by a factor of $\sim 5$, increasing the mass accretion rate by two orders of magnitude to more than $10^{-4}$ $M_{\odot}$ yr$^{-1}$. The layer becomes highly filamentary due to gas flows along the magnetic field lines, and dense cores are formed in the filaments. The massive dense cores have size and mass of 0.03 -- 0.8 pc and 8 -- 50 $M_{\odot}$ and they are usually gravitationally unstable. The mass function of the dense cores is significantly top-heavy as compared with the universal IMF, indicating that the cloud-cloud collision triggers preferentially the formation of O and early B stars. We argue that the cloud-cloud collision is a versatile mechanism which creates a variety of stellar clusters from a single O star like RCW120 and M20 to tens of O stars of a super star cluster like RCW38 and a mini-starburst W43. The core mass function predicted by the present model is consistent with the massive dense cores obtained by recent ALMA observations in RCW38 (Torii et al. 2019) and W43 (Motte et al. 2018) considering the increasing evidence for collision-triggered high-mass star formation, we argue that cloud-cloud collision is a major mechanism of high mass star formation.

preprint2020arXiv

The formation of young massive clusters triggered by cloud-cloud collisions in the Antennae Galaxies NGC 4038/NGC 4039

The formation mechanism of super star clusters (SSCs), a present-day analog of the ancient globulars, still remains elusive. The major merger, the Antennae galaxies is forming SSCs and is one of the primary targets to test the cluster formation mechanism. We reanalyzed the archival ALMA CO data of the Antennae and found three typical observational signatures of a cloud-cloud collision toward SSC B1 and other SSCs in the overlap region; i. two velocity components with $\sim$100 km s$^{-1}$ velocity separation, ii. the bridge features connecting the two components, and iii. the complementary spatial distribution between them, lending support for collisions of the two components as a cluster formation mechanism. We present a scenario that the two clouds with 100 km s$^{-1}$ velocity separation collided, and SSCs having $\sim$10$^6$-10$^7$ $M_{\rm \odot}$ were formed rapidly during the time scale. {We compared the present results with the recent studies of star forming regions in the Milky Way and the LMC, where the SSCs having $\sim$10$^4$-10$^5$ $M_{\rm \odot}$ are located. As a result, we found that there is a positive correlation between the compressed gas pressure generated by collisions and the total stellar mass of SSC, suggesting that the pressure may be a key parameter in the SSC formation.

preprint2020arXiv

Triggered high-mass star formation in the HII region W28A2: A cloud-cloud collision scenario

We report on a study of the high-mass star formation in the the HII region W28A2 by investigating the molecular clouds extended over ~5-10 pc from the exciting stars using the 12CO and 13CO (J=1-0) and 12CO (J=2-1) data taken by the NANTEN2 and Mopra observations. These molecular clouds consist of three velocity components with the CO intensity peaks at V_LSR ~ -4 km s$^{-1}$, 9 km s$^{-1}$ and 16 km s$^{-1}$. The highest CO intensity is detected at V_LSR ~ 9 km s$^{-1}$, where the high-mass stars with the spectral types of O6.5-B0.5 are embedded. We found bridging features connecting these clouds toward the directions of the exciting sources. Comparisons of the gas distributions with the radio continuum emission and 8 um infrared emission show spatial coincidence/anti-coincidence, suggesting physical associations between the gas and the exciting sources. The 12CO J=2-1 to 1-0 intensity ratio shows a high value (> 0.8) toward the exciting sources for the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds, possibly due to heating by the high-mass stars, whereas the intensity ratio at the CO intensity peak (V_LSR ~ 9 km s$^{-1}$) lowers down to ~0.6, suggesting self absorption by the dense gas in the near side of the +9 km s$^{-1}$ cloud. We found partly complementary gas distributions between the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds, and the -4 km s$^{-1}$ and +16 km s$^{-1}$ clouds. The exciting sources are located toward the overlapping region in the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds. Similar gas properties are found in the Galactic massive star clusters, RCW 38 and NGC 6334, where an early stage of cloud collision to trigger the star formation is suggested. Based on these results, we discuss a possibility of the formation of high-mass stars in the W28A2 region triggered by the cloud-cloud collision.

preprint2019arXiv

Cloud-cloud collisions in the common foot point of molecular loops 1 and 2 in the Galactic Center

Recent large-area, deep CO surveys in the Galactic disk have revealed the formation of ~50 high-mass stars or clusters triggered by cloud-cloud collisions (CCCs). Although the Galactic Center (GC) -- which contains the highest volume density of molecular gas -- is the most favorable place for cloud collisions, systematic studies of CCCs in that region are still untouched. Here we report for the first time evidence of CCCs in the common foot point of molecular loops 1 and 2 in the GC. We have investigated the distribution of molecular gas toward the foot point by using a methodology for identifying CCCs, and we have discovered clear signatures of CCCs. Using the estimated displacements and relative velocities of the clouds, we find the elapsed time since the beginnings of the collisions to be 105-6 yr. We consider possible origins for previously reported peculiar velocity features in the foot point and discuss star formation triggered by CCCs in the GC.

preprint2019arXiv

Cluster formation induced by a cloud--cloud collision in [DBS2003]179

[DBS2003]179 is a super star cluster in the Galaxy discovered by deep near infrared observations. We carried out CO J=1-0 and J=3-2 observations of the region of [DBS2003]179 with NANTEN2, ASTE and the Mopra 22m telescope. We identified and mapped two molecular clouds which are likely associated with the cluster. The association is evidenced by the spatial correlation with the 8 micron Spitzer image, and a high ratio of the two transitions of 12CO (J=3-2 to J=1-0). The two clouds show complementary distribution in space and bridge features connecting them in velocity. We frame a hypothesis that the two clouds collided with each other 1-2 Myr ago, and the collision compressed the interface layer, triggering the formation of the cluster. This offers an additional piece of evidence for a super star cluster formed by cloud--cloud collision alongside of the four super star clusters including Wd2, NGC3603, RCW38 and R136. The findings indicate that the known super star clusters having closely associated dust emission are formed by cloud-cloud collision, lending support for the important role of cloud--cloud collision in high-mass star formation.

preprint2019arXiv

High-mass star formation in Orion possibly triggered by cloud-cloud collision III, NGC2068 and NGC2071

Using the NANTEN2 Observatory, we carried out a molecular line study of high-mass star forming regions with reflection nebulae, NGC 2068 and NGC 2071, in Orion in the 13CO(J=2-1) transition. The 13CO distribution shows that there are two velocity components at 9.0 and 10.5 km/s . The blue-shifted component is in the northeast associated with NGC 2071, whereas the red-shifted component is in the southwest associated with NGC 2068. The total intensity distribution of the two clouds shows a gap of ~1 pc, suggesting that they are detached at present. A detailed spatial comparison indicates that the two show complementary distributions. The blue-shifted component lies toward an intensity depression to the northwest of the red-shifted component, where we find that a displacement of 0.8 pc makes the two clouds fit well with each other. Furthermore, a new simulation of non-frontal collisions shows that observations from 60 degrees off the collisional axis agreed well with the velocity structure in this region. On the basis of these results, we hypothesize that the two components collided with each other at a projected relative velocity 3.0 km/s estimated to be 0.3 Myr for an assumed axis of the relative motion 60 degrees off the line of sight. We assume that the two most massive early B-type stars in the cloud, illuminating stars of the two reflection nebulae, were formed by collisional triggering at the interfaces between the two clouds. Given the other young high-mass star forming regions, namely, M42, M43, and NGC 2024 (Fukui et al. 2018b; Ohama et al. 2017a), it seems possible that collisional triggering has been independently working to form O-type and early B-type stars in Orion in the last Myr over a projected distance of ~80 pc.

preprint2019arXiv

Relations Between Molecular Cloud Structure Sizes and Line Widths in the Large Magellanic Cloud

We present a comparative study of the size-line width relation for substructures within six molecular clouds in the Large Magellanic Cloud (LMC) mapped with the Atacama Large Millimeter/submillimeter Array (ALMA). Our sample extends our previous study, which compared a Planck detected cold cloud in the outskirts of the LMC with the 30 Doradus molecular cloud and found the typical line width for 1 pc radius structures to be 5 times larger in 30 Doradus. By observing clouds with intermediate levels of star formation activity, we find evidence that line width at a given size increases with increasing local and cloud-scale 8$μ$m intensity. At the same time, line width at a given size appears to independently correlate with measures of mass surface density. Our results suggest that both virial-like motions due to gravity and local energy injection by star formation feedback play important roles in determining intracloud dynamics.

preprint2016arXiv

A quality check of the $AKARI$ mid-infrared all-sky diffuse map toward the massive star-forming regions NGC 6334 and NGC 6357

We present a comparative study of CO and polycyclic aromatic hydrocarbon (PAH) emission toward a region including the massive star-forming regions of NGC 6334 and NGC 6357. We use the NANTEN $^{12}$CO($J$ = 1--0) data and the $AKARI$ 9 $μ$m All-Sky diffuse map in order to evaluate the calibration accuracy of the $AKARI$ data. We confirm that the overall CO distribution shows a good spatial correspondence with the PAH emission, and their intensities exhibit a good power-law correlation with a spatial resolution down to 4$'$ over the region of 10$^\circ$$\times$10$^\circ$. We also reveal poorer correlation for small scale structures between the two quantities toward NGC 6357, due to strong UV radiation from local sources. Larger scatter in the correlation toward NGC 6357 indicates higher ionization degree and/or PAH excitation than that of NGC 6334.

preprint2016arXiv

Ammonia excitation imaging of shocked gas towards the W28 gamma-ray source HESS J1801-233

We present 12 mm Mopra observations of the dense (>10^3 cm^-3 ) molecular gas towards the north-east (NE) of the W28 supernova remnant (SNR). This cloud is spatially well-matched to the TeV gamma-ray source HESS J1801-233 and is known to be a SNR-molecular cloud interaction region. Shock-disruption is evident from broad NH3 (1,1) spectral line-widths in regions towards the W28 SNR, while strong detections of spatially-extended NH3(3,3), NH3(4,4) and NH3(6,6) inversion emission towards the cloud strengthen the case for the existence of high temperatures within the cloud. Velocity dispersion measurements and NH3(n,n)/(1,1) ratio maps, where n=2, 3, 4 and 6, indicate that the source of disruption is from the side of the cloud nearest to the W28 SNR, suggesting that it is the source of cloud-disruption. Towards part of the cloud, the ratio of ortho to para-NH3 is observed to exceed 2, suggesting gas-phase NH3 enrichment due to NH3 liberation from dust grain mantles. The measured NH3 abundance with respect to H2 is ~(1.2+/-0.5)*10^-9, which is not high, as might be expected for a hot, dense molecular cloud enriched by sublimated grain-surface molecules. The results are suggestive of NH3 sublimation and destruction in this molecular cloud, which is likely to be interacting with the W28 SNR shock.

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

Instrumental performance and results from testing of the BLAST-TNG receiver, submillimeter optics, and MKID arrays

Polarized thermal emission from interstellar dust grains can be used to map magnetic fields in star forming molecular clouds and the diffuse interstellar medium (ISM). The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) flew from Antarctica in 2010 and 2012 and produced degree-scale polarization maps of several nearby molecular clouds with arcminute resolution. The success of BLASTPol has motivated a next-generation instrument, BLAST-TNG, which will use more than 3000 linear polarization sensitive microwave kinetic inductance detectors (MKIDs) combined with a 2.5m diameter carbon fiber primary mirror to make diffraction-limited observations at 250, 350, and 500 $μ$m. With 16 times the mapping speed of BLASTPol, sub-arcminute resolution, and a longer flight time, BLAST-TNG will be able to examine nearby molecular clouds and the diffuse galactic dust polarization spectrum in unprecedented detail. The 250 $μ$m detector array has been integrated into the new cryogenic receiver, and is undergoing testing to establish the optical and polarization characteristics of the instrument. BLAST-TNG will demonstrate the effectiveness of kilo-pixel MKID arrays for applications in submillimeter astronomy. BLAST-TNG is scheduled to fly from Antarctica in December 2017 for 28 days and will be the first balloon-borne telescope to offer a quarter of the flight for "shared risk" observing by the community.

preprint2016arXiv

Mid- and far-infrared properties of Spitzer Galactic bubbles revealed by the AKARI all-sky surveys

We have carried out a statistical study on the mid- and far-infrared (IR) properties of Galactic IR bubbles observed by Spitzer. Using the Spitzer 8 $μ{\rm m}$ images, we estimated the radii and covering fractions of their shells, and categorized them into closed, broken and unclassified bubbles with our data analysis method. Then, using the AKARI all-sky images at wavelengths of 9, 18, 65, 90, 140 and 160 $μ{\rm m}$, we obtained the spatial distributions and the luminosities of polycyclic aromatic hydrocarbon (PAH), warm and cold dust components by decomposing 6-band spectral energy distributions with model fitting. As a result, 180 sample bubbles show a wide range of the total IR luminosities corresponding to the bolometric luminosities of a single B-type star to many O-type stars. For all the bubbles, we investigated relationships between the radius, luminosities and luminosity ratios, and found that there are overall similarities in the IR properties among the bubbles regardless of their morphological types. In particular, they follow a power-law relation with an index of $\sim$3 between the total IR luminosity and radius, as expected from the conventional picture of the Str$\rm{\ddot{o}}$mgren sphere. The exceptions are large broken bubbles; they indicate higher total IR luminosities, lower fractional luminosities of the PAH emission, and dust heating sources located nearer to the shells. We discuss the implications of those differences for a massive star-formation scenario.

preprint2016arXiv

Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. I. N159W

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

preprint2016arXiv

Revealing a detailed mass distribution of a high-density core MC27/L1521F in Taurus with ALMA

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

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.

preprint2015arXiv

A dense gas survey of the gamma-ray sources HESS J1731-347 and HESS J1729-345

The results of Mopra molecular spectral line observations towards the supernova remnant HESSJ1731-347 (G353.6-0.7) and the unidentified gamma-ray source HESSJ1729-345 are presented. Dense molecular gas in three different velocity-bands (corresponding to three Galactic arms) are investigated using the CS(1-0) line. The CS-traced component provides information about the dense target material in a hadronic scenario for gamma-ray production (cosmic rays interacting with gas) and an understanding of the dynamics. Furthermore, the effects of cosmic ray diffusion into dense gas may alter the gamma-ray spectrum to cause a flattening of spectra towards such regions. Dense molecular gas mass at a level of ~10^5 Mo was revealed in this survey, with mass of the order of ~10^3 Mo towards HESSJ1729-345 in each coincident Galactic arm, but no significant detection of dense molecular gas towards HESSJ1731-347 at the currently-preferred distance of ~5.2-6.2 kpc was discovered.

preprint2015arXiv

Evidence for Thermal X-Ray Line Emission from the Synchrotron-Dominated Supernova Remnant RX J1713.7-3946

We report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

preprint2015arXiv

High-mass star formation triggered by collision between CO filaments in N159 West in the Large Magellanic Cloud

We have carried out 13CO(J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0pc and a length of 5-10pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37 Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~6x10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.

preprint2015arXiv

Stochastic Noncircular Motion and Outflows Driven by Magnetic Activity in the Galactic Bulge Region

By performing a global magneto-hydrodynamical simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic center region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. In addition, the magnetic pressure-gradient force also drives radial flows in a similar manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram-shape of the central molecular zone. This is a natural extension into the central few 100 pc of the magnetic activity, which is observed as molecular loops at radii from a few 100 pc to 1 kpc. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time-dependent, which we discuss from a viewpoint of the outflow from the bulge.

preprint2014arXiv

ALMA observations of a high-density core in Taurus: dynamical gas interaction at the possible site of a multiple star formation

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

preprint2014arXiv

Dense Molecular Clumps associated with the LMC Supergiant Shells LMC 4 \& LMC 5

We investigate the effects of Supergiant Shells (SGSs) and their interaction on dense molecular clumps by observing the Large Magellanic Cloud (LMC) star forming regions N48 and N49, which are located between two SGSs, LMC 4 and LMC 5. $^{12}$CO ($J$=3-2, 1-0) and $^{13}$CO ($J$=1-0) observations with the ASTE and Mopra telescopes have been carried out towards these regions. A clumpy distribution of dense molecular clumps is revealed with 7 pc spatial resolution. Large velocity gradient analysis shows that the molecular hydrogen densities ($n({\rm H}_2)$) of the clumps are distributed from low to high density ($10^3$-$10^5$ cm$^{-3}$) and their kinetic temperatures ($T_{\rm kin}$) are typically high (greater than $50$ K). These clumps seem to be in the early stages of star formation, as also indicated from the distribution of H$α$, young stellar object candidates, and IR emission. We found that the N48 region is located in the high column density HI envelope at the interface of the two SGSs and the star formation is relatively evolved, whereas the N49 region is associated with LMC 5 alone and the star formation is quiet. The clumps in the N48 region typically show high $n({\rm H}_2)$ and $T_{\rm kin}$, which are as dense and warm as the clumps in LMC massive cluster-forming areas (30 Dor, N159). These results suggest that the large-scale structure of the SGSs, especially the interaction of two SGSs, works efficiently on the formation of dense molecular clumps and stars.

preprint2014arXiv

Detailed Investigation of the Gamma-Ray Emission in the Vicinity of SNR W28 with FERMI-LAT

We present a detailed investigation of the $γ$-ray emission in the vicinity of the supernova remnant (SNR) W28 (G6.4$-$0.1) observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. We detected significant $γ$-ray emission spatially coincident with TeV sources HESS J1800$-$240A, B, and C, located outside the radio boundary of the SNR. Their spectra in the 2-100 GeV band are consistent with the extrapolation of the power-law spectra of the TeV sources. We also identified a new source of GeV emission, dubbed Source W, which lies outside the boundary of TeV sources and coincides with radio emission from the western part of W28. All of the GeV $γ$-ray sources overlap with molecular clouds in the velocity range from 0 to 20 km s$^{-1}$. Under the assumption that the $γ$-ray emission towards HESS J1800-240A, B, and C comes from $π^0$ decay due to the interaction between the molecular clouds and cosmic rays (CRs) escaping from W28, they can be naturally explained by a single model in which the CR diffusion coefficient is smaller than the theoretical expectation in the interstellar space. The total energy of the CRs escaping from W28 is constrained through the same modeling to be larger than $\sim$ 2 $\times$ 10$^{49}$ erg. The emission from Source W can also be explained with the same CR escape scenario.

preprint2014arXiv

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

preprint2014arXiv

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. II. Gas-to-Dust Ratio Variations across ISM Phases

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

preprint2014arXiv

HI, CO, and Planck/IRAS dust properties in the high-latitude-cloud complex, MBM 53, 54, 55 and HLCG 92-35; Possible evidence for an optically thick HI envelope around the CO clouds

We present an analysis of the HI and CO gas in conjunction with the Planck/IRAS submillimeter/far-infrared dust properties toward the most outstanding high latitude clouds MBM 53, 54, 55 and HLCG 92-35 at b = -30 deg to -45 deg. The CO emission, dust opacity at 353 GHz (tau353), and dust temperature (Td) show generally good spatial correspondence. On the other hand, the correspondence between the HI emission and the dust properties is less clear than in CO. The integrated HI intensity WHI and tau353 show a large scatter with a correlation coefficient of ~0.6 for a Td range from 16 K to 22 K. We find, however, that WHI and tau353 show better correlation for smaller ranges of Td every 0.5 K, generally with a correlation coefficient of 0.7-0.9. We set up a hypothesis that the HI gas associated with the highest Td >= 21.5 K is optically thin, whereas the HI emission is generally optically thick for Td lower than 21.5 K. We have determined a relationship for the optically thin HI gas between atomic hydrogen column density and tau353, NHI (cm-2) = (1.5 x 10^26) x tau353, under the assumption that the dust properties are uniform and we have applied this to estimate NHI from tau353 for the whole cloud. NHI was then used to solve for Ts and tauHI over the region. The result shows that the HI is dominated by optically thick gas having a low spin temperature of 20-40 K and a density of 40-160 cm-3. The HI envelope has a total mass of ~1.2 x 10^4 Msol, an order of magnitude larger than that of the CO clouds. The HI envelope properties derived by this method do not rule out a mixture of HI and H2 in the dark gas, but we present indirect evidence that most of the gas mass is in the atomic state.

preprint2014arXiv

Revealing the physical properties of molecular gas in Orion with a large scale survey in J=2-1 lines of 12CO, 13CO and C18O

We present fully sampled ~3' resolution images of the 12CO(J=2-1), 13CO(J=2-1), and C18O(J=2-1) emission taken with the newly developed 1.85-m mm-submm telescope toward the entire area of the Orion A and B giant molecular clouds. The data were compared with the J=1-0 of the 12CO, 13CO, and C18O data taken with the Nagoya 4-m telescope and the NANTEN telescope at the same angular resolution to derive the spatial distributions of the physical properties of the molecular gas. We explore the large velocity gradient formalism to determine the gas density and temperature by using the line combinations of 12CO(J=2-1), 13CO(J=2-1), and 13CO(J=1-0) assuming uniform velocity gradient and abundance ratio of CO. The derived gas density is in the range of 500 to 5000 cm-3, and the derived gas temperature is mostly in the range of 20 to 50 K along the cloud ridge with a temperature gradient depending on the distance from the star forming region. We found the high-temperature region at the cloud edge facing to the HII region, indicating that the molecular gas is interacting with the stellar wind and radiation from the massive stars. In addition, we compared the derived gas properties with the Young Stellar Objects distribution obtained with the Spitzer telescope to investigate the relationship between the gas properties and the star formation activity therein. We found that the gas density and star formation efficiency are well positively correlated, indicating that stars form effectively in the dense gas region.

preprint2014arXiv

Spitzer View of Massive Star Formation in the Tidally Stripped Magellanic Bridge

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

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

CO Component Estimation Based on the Independent Component Analysis

Fast Independent Component Analysis (FastICA) is a component separation algorithm based on the levels of non-Gaussianity. Here we apply the FastICA to the component separation problem of the microwave background including carbon monoxide (CO) line emissions that are found to contaminate the PLANCK High Frequency Instrument (HFI) data. Specifically we prepare 100GHz, 143GHz, and 217GHz mock microwave sky maps including galactic thermal dust, NANTEN CO line, and the Cosmic Microwave Background (CMB) emissions, and then estimate the independent components based on the kurtosis. We find that the FastICA can successfully estimate the CO component as the first independent component in our deflection algorithm as its distribution has the largest degree of non-Gaussianity among the components. By subtracting the CO and the dust components from the original sky maps, we will be able to make an unbiased estimate of the cosmological CMB angular power spectrum.

preprint2013arXiv

Dense Gas Towards the RXJ1713.7-3946 Supernova Remnant

We present results from a Mopra 7mm-wavelength survey that targeted the dense gas-tracing CS(1-0) transition towards the young gamma-ray-bright supernova remnant, RXJ1713.7-3946 (SNR G347.3-0.5). In a hadronic gamma-ray emission scenario, where cosmic ray protons interact with gas to produce the observed gamma-ray emission, the mass of potential cosmic ray target material is an important factor. We summarise newly-discovered dense gas components, towards Cores G and L, and Clumps N1, N2, N3 and T1, which have masses of 1-10^4 solar masses. We argue that these components are not likely to contribute significantly to gamma-ray emission in a hadronic gamma-ray emission scenario. This would be the case if RXJ1713.7-3946 were at either the currently favoured distance of ~1kpc or an alternate distance (as suggested in some previous studies) of ~6kpc. This survey also targeted the shock-tracing SiO molecule. Although no SiO emission corresponding to the RXJ1713.7-3946 shock was observed, vibrationally-excited SiO(1-0) maser emission was discovered towards what may be an evolved star. Observations taken one year apart confirmed a transient nature, since the intensity, line-width and central velocity of SiO(J=1-0,v=1,2) emission varied significantly.

preprint2013arXiv

Detection of an ultra-bright submillimeter galaxy behind the Small Magellanic Cloud

We report the discovery of a new ultra-bright submillimeter galaxy (SMG) behind the Small Magellanic Cloud (SMC). This SMG is detected as a 43.3+-8.4 mJy point source (MM J01071-7302, hereafter MMJ0107) in the 1.1 mm continuum survey of the SMC by AzTEC on the ASTE telescope. MMJ0107 is also detected in the radio (843 MHz), Herschel/SPIRE, Spitzer MIPS 24 μm, all IRAC bands, Wide-field Infrared Survey Explorer, and near-infrared (J, H, KS). We find an optical (U, B, V) source, which might be the lensing object, at a distance of 1.4 arcsec from near-infrared and IRAC sources. Photometric redshift estimates for the SMG using representative spectral energy distribution templates show the redshifts of 1.4-3.9. We estimate total far-infrared luminosity of (0.3-2.2)x10^14 μ^-1 L_sun and a star formation rate of 5600-39, 000 μ^-1 M_sun yr^-1, where μ is the gravitational magnification factor. This apparent extreme star formation activity is likely explained by a highly magnified gravitational lens system.

preprint2013arXiv

Formation of Massive Molecular Cloud Cores by Cloud-cloud Collision

Recent observations of molecular clouds around rich massive star clusters including NGC3603, Westerlund 2, and M20 revealed that the formation of massive stars could be triggered by a cloud-cloud collision. By using three-dimensional, isothermal, magnetohydrodynamics simulations with the effect of self-gravity, we demonstrate that massive, gravitationally unstable, molecular cloud cores are formed behind the strong shock waves induced by the cloud-cloud collision. We find that the massive molecular cloud cores have large effective Jeans mass owing to the enhancement of the magnetic field strength by shock compression and turbulence in the compressed layer. Our results predict that massive molecular cloud cores formed by the cloud-cloud collision are filamentary and threaded by magnetic fields perpendicular to the filament.

preprint2013arXiv

Formation of Molecular Clouds and Global Conditions for Star Formation

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

preprint2013arXiv

Interstellar gas towards CTB 37A and the TeV gamma-ray source HESS J1714-385

Observations of dense molecular gas towards the supernova remnants CTB 37A (G348.5+0.1) and G348.5-0.0 were carried out using the Mopra and Nanten2 radio telescopes. We present CO(2-1) and CS(1-0) emission maps of a region encompassing the CTB 37A TeV gamma-ray emission, HESS J1714-385, revealing regions of dense gas within associated molecular clouds. Some gas displays good overlap with gamma-ray emission, consistent with hadronic gamma-ray emission scenarios. Masses of gas towards the HESS J1714-385 TeV gamma-ray emission region were estimated, and were of the order of 10^3-10^4 solar masses. In the case of a purely hadronic origin for the gamma-ray emission, the cosmic ray flux enhancement is ~80-1100 times the local solar value. This enhancement factor and other considerations allow a discussion of the age of CTB 37A, which is consistent with ~10^4 yr.

preprint2013arXiv

Lupus I Observations from the 2010 Flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 μm. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

preprint2013arXiv

Molecular and atomic gas in the young TeV γ-ray SNRs RX J1713.7-3946 and RX J0852.0-4622; evidence for the hadronic production of γ-rays

The interstellar molecular clouds are the site of star formation and also the target for the cosmic ray protons to produce γ-rays via the hadronic process. The interstellar atomic gas is enveloping the molecular clouds and may also be dense enough to affect the γ-ray production. In this Chapter, some of the basic properties of the interstellar gas both in molecular and atomic forms will be reviewed. Then, it is presented that two young TeV γ-ray SNRs, RX J1713.7-3946 and RX J0852.0-4622, show good spatial correspondence between the γ-rays and the interstellar protons. The good spatial correspondence provides a support for the hadronic origin of the γ-rays in these SNRs. It is emphasized that both molecular and atomic hydrogen plays a role as targets for cosmic ray (CR) protons. The clumpy distribution of the target interstellar medium (ISM) protons is crucial in the interaction of the supernova shocks with the ISM, whereas models with uniform ISM distribution are not viable. Finally, it is suggested that the dense atomic gas without molecules may occupy the dominant part of the dark gas in the local ISM.

preprint2013arXiv

Molecular Clouds Toward the Super Star Cluster NGC3603; Possible Evidence for a Cloud-Cloud Collision in Triggering the Cluster Formation

We present new large field observations of molecular clouds with NANTEN2 toward the super star cluster NGC3603 in the transitions 12CO(J=2-1, J=1-0) and 13CO(J=2-1, J=1-0). We suggest that two molecular clouds at 13 km s-1 and 28 km s-1 are associated with NGC3603 as evidenced by higher temperatures toward the H II region as well as morphological correspondence. The mass of the clouds is too small to gravitationally bind them, given their relative motion of ~20 km s-1. We suggest that the two clouds collided with each other a Myr ago to trigger the formation of the super star cluster. This scenario is able to explain the origin of the highest mass stellar population in the cluster which is as young as a Myr and is segregated within the central sub-pc of the cluster. This is the second super star cluster along side Westerlund2 where formation may have been triggered by a cloud-cloud collision.

preprint2012arXiv

3 to 12 millimetre studies of dense gas towards the western rim of supernova remnant RX J1713.7-3946

The young X-ray and gamma-ray-bright supernova remnant RXJ1713.7-3946 (SNR G347.3-0.5) is believed to be associated with molecular cores that lie within regions of the most intense TeV emission. Using the Mopra telescope, four of the densest cores were observed using high-critical density tracers such as CS(J=1-0,J=2-1) and its isotopologue counterparts, NH3(1,1) and (2,2) inversion transitions and N2H+(J=1-0) emission, confirming the presence of dense gas >10^4cm^-3 in the region. The mass estimates for Core C range from 40M_{\odot} (from CS(J=1-0)) to 80M_{\odot} (from NH3 and N2H+), an order of magnitude smaller than published mass estimates from CO(J=1-0) observations. We also modelled the energy-dependent diffusion of cosmic-ray protons accelerated by RXJ1713.7-3946 into Core C, approximating the core with average density and magnetic field values. We find that for considerably suppressed diffusion coefficients (factors χ=10^{-3} down to 10^{-5} the galactic average), low energy cosmic-rays can be prevented from entering the inner core region. Such an effect could lead to characteristic spectral behaviour in the GeV to TeV gamma-ray and multi-keV X-ray fluxes across the core. These features may be measurable with future gamma-ray and multi-keV telescopes offering arcminute or better angular resolution, and can be a novel way to understand the level of cosmic-ray acceleration in RXJ1713.7-3946 and the transport properties of cosmic-rays in the dense molecular cores.

preprint2012arXiv

Recombining Plasma and Hard X-ray Filament in the Mixed-Morphology Supernova Remnant W44

We report new features of the typical mixed-morphology (MM) supernova remnant (SNR) W44. In the X-ray spectra obtained with Suzaku, radiative recombination continua (RRCs) of highly ionized atoms are detected for the first time. The spectra are well reproduced by a thermal plasma in a recombining phase. The best-fit parameters suggest that the electron temperature of the shock-heated matters cooled down rapidly from $\sim1$,keV to $\sim 0.5$,keV, possibly due to adiabatic expansion (rarefaction) occurred $\sim20,000$ years ago. We also discover hard X-ray emission which shows an arc-like structure spatially-correlated with a radio continuum filament. The surface brightness distribution shows a clear anti-correlation with $^{12}$CO (J=2-1) emission from a molecular cloud observed with NANTEN2. While the hard X-ray is most likely due to a synchrotron enhancement in the vicinity of the cloud, no current model can quantitatively predict the observed flux.

preprint2012arXiv

The Carina Flare: What can fragments in the wall tell us?

$^{13}$CO(J=2--1) and C$^{18}$O(J=2--1) observations of the molecular cloud G285.90+4.53 (Cloud~16) in the Carina Flare supershell (GSH287+04-17) with the APEX telescope are presented. With an algorithm DENDROFIND we identify 51 fragments and compute their sizes and masses. We discuss their mass spectrum and interpret it as being the result of the shell fragmentation process described by the pressure assisted gravitational instability - PAGI. We conclude that the explanation of the clump mass function needs a combination of gravity with pressure external to the shell.

preprint2011arXiv

A Pan-Carina YSO Catalog: Intermediate-Mass Young Stellar Objects in the Carina Nebula Identified Via Mid-Infrared Excess Emission

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg^2 field surveyed by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data from the Spitzer Space Telescope Vela--Carina survey and the Two-Micron All Sky Survey. We model the 1--24 /mu m IR spectral energy distributions of the YSOs to constrain physical properties. Our Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 Msun < m < 10 Msun) objects with disks, including Herbig Ae/Be stars and their less evolved progenitors. The PCYC provides a valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is consistent with the presence of low-mass companions, but we also find that X-ray emission correlates with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo activity during the convective atmosphere phase, but this emission dies off as the stars approach the main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, we predict a total population of >2x10^4 YSOs and a present-day star formation rate (SFR) of >0.008 Msun/yr. The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been approximately constant.

preprint2011arXiv

Molecular and Atomic Gas toward HESS J1745-303 in the Galactic Center: Further Support for the Hadronic Scenario

We have compared the TeV gamma-rays with the new 12CO J=2-1 data toward HESS J1745-303 in the Galactic center and confirmed that the molecular gas MG358.9-0.5 toward (l, b)=(358.9, -0.5 at VLSR=-100-0 km s-1 shows a reasonable positional agreement with the primary peak (northern part) of the gamma-ray source. For the southern part of HESS J1745-303, we see no CO counterpart, whereas the HI gas in the Parkes 21 cm HI dataset shows a possible counterpart to the gamma-ray source. This HI gas may be optically thick as supported by the HI line shape similar to the optically thick 12CO. We estimate the total mass of interstellar protons including both the molecular and atomic gas to be 2x10^6 Mo and the cosmic-ray proton energy to be 6x10^{48} ergs in the hadronic scenario. We discuss possible origins of the cosmic-ray protons including the nearby SNR G359.1-0.5. The SNR may be able to explain the northern gamma-ray source but the southern source seems to be too far to be energized by the SNR. As an alternative, we argue that the second-order Fermi acceleration in the inter-clump space surrounded by randomly moving high-velocity clumps may offer a possible mechanism to accelerate protons. The large turbulent motion with velocity dispersion of ~15 km s-1 has energy density two orders of magnitude higher than in the solar vicinity and is viable as the energy source.

preprint2011arXiv

Molecular Clouds as Cosmic Ray Laboratories

We will here discuss how the gamma-ray emission from molecular clouds can be used to probe the cosmic ray flux in distant regions of the Galaxy and to constrain the highly unknown cosmic ray diffusion coefficient. In particular we will discuss the GeV to TeV emission from runaway cosmic rays penetrating molecular clouds close to young and old supernova remnants and in molecular clouds illuminated by the background cosmic ray flux.

preprint2011arXiv

Processing of polycyclic aromatic hydrocarbons in molecular-loop regions near the Galactic center revealed by AKARI

We present the AKARI mid-infrared diffuse map of an area of about 4 deg x 3 deg near the Galactic center in the 9 um band. The band intensity is mostly attributed to the aromatic hydrocarbon infrared emissions of carbonaceous grains at wavelengths of 6.2, 7,7, 8.6, and 11.3 um. We detect the 9 um emission structures extending from the Galactic plane up to the latitude of about 2.5 deg, which have spatial correspondence to the molecular loops revealed by the NANTEN 12CO (J=1-0) observations. We find that the surface brightness at 9 um is suppressed near the foot points of the CO loops. The ratios of the 9 um to the IRAS 100 um brightness show significant depression near such bright regions in the CO emission. With the AKARI near-IR 2.5--5 um spectroscopy, we find that the 3.3 um aromatic hydrocarbon emission is absent in the region associated with the loop. These suggest the processing and destruction of carbonaceous grains in the CO molecular loops.

preprint2011arXiv

Source Contamination in X-ray Studies of Star-Forming Regions: Application to the Chandra Carina Complex Project

We describe detailed simulations of X-ray-emitting populations to evaluate the levels of contamination by both Galactic and extragalactic X-ray sources unrelated to a star-forming region under study. For Galactic contaminations, we consider contribution from main-sequence stars and giants (not including cataclysmic variables and other classes of accretion-driven X-ray binary systems) as they make the dominant contribution at the position of the Carina Nebula. The simulations take into consideration a variety of technical factors involving a Galactic population synthesis model, stellar X-ray luminosity functions, Chandra telescope response, source detection methodology, and possible spatial variations in the X-ray background and absorption through molecular clouds. When applied to the 1.42 square-degree field of the Chandra Carina Complex Project (CCCP), the simulations predict ~5000 contaminating sources (1 source per square arcminute of the survey), evenly distributed across the field. The results of the simulations are further employed in a companion CCCP study to assign membership probabilities to individual sources.

preprint2011arXiv

Spitzer Characterization of Dust in the Ionized Medium of the Large Magellanic Cloud

A systematic investigation of dust emission associated with the ionized gas has so far been performed only in our Galaxy and for wavelengths longer than 60 μm. Newly available Spitzer data now offer the opportunity to carry out a similar analysis in the Large Magellanic Cloud (LMC). By cross-correlating Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) data with the ATCA/Parkes HI 21-cm data, the NANTEN 12CO (J=1-0) data, and both the SHASSA Hα and the Parkes 6-cm data, we investigate the physical properties of dust associated with the different phases of the gas (atomic, molecular and ionized). In particular, we study the presence and nature of dust from 3.6 to 160 μm and for various regimes of the ionized gas, spanning emission measures (EM) from \sim 1 pc cm-6 (diffuse component) to \sim 10^3 pc cm-6 (HII regions). Using a dust emission model, and testing our results with several radiation field spectra, we show that dust in the ionized gas is warmer than dust associated with other phases (atomic and molecular). We also find a decrease of the polycyclic aromatic hydrocarbons (PAH) relative abundance with respect to big grains (BGs), as well as an increase of the near infrared (NIR) continuum. These three results (e.g. warmer temperature, decrease of PAH abundance and increase of the NIR continuum) are found consistently for all regimes of the ionized gas. On the contrary, the molecular phase appears to provide favorable conditions for the survival of PAHs. Furthermore, the very small grain (VSG) relative abundance tends to increase in the ionized phase, especially in bright HII regions. Last but not least, our analysis shows that the emissivity of dust associated with the ionized gas is lower in the LMC than in our Galaxy, and that this difference is not accounted for by the lower metallicity of the LMC.

preprint2011arXiv

Stochastic Acceleration of Cosmic Rays in the Central Molecular Zone of the Galaxy

Particle acceleration in the inner ~ 200 pc of the Galaxy is discussed, where diffuse TeV gamma-rays have been detected by the High Energy Stereoscopic System (HESS) observation. The diffuse gamma-ray emission has a strong correlation with molecular clouds with large velocity dispersion, indicating the presence of turbulence. It is argued that the turbulence may contribute to the acceleration of cosmic rays via stochastic acceleration. The stochastic acceleration may energize cosmic-ray protons up to ~ 100 TeV and electrons to ~ 1 TeV in a relatively tenuous medium. The difference in the efficiency between protons and electrons supports the hadronic scenario of the diffuse TeV gamma-ray emission.

preprint2011arXiv

The CO-to-H2 Conversion Factor From Infrared Dust Emission Across the Local Group

We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are α_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and α_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.

preprint2011arXiv

The Galactic Census of High- and Medium-mass Protostars. I. Catalogues and First Results from Mopra HCO+ Maps

The Census of High- and Medium-mass Protostars (CHaMP) is the first large-scale, unbiased, uniform mapping survey at sub-parsec scale resolution of 90 GHz line emission from massive molecular clumps in the Milky Way. We present the first Mopra (ATNF) maps of the CHaMP survey region (300°>l>280°) in the HCO+ J=1-0 line, which is usually thought to trace gas at densities up to 10^11 m-3. In this paper we introduce the survey and its strategy, describe the observational and data reduction procedures, and give a complete catalogue of moment maps of the HCO+ J=1-0 emission from the ensemble of 301 massive molecular clumps. From these maps we also derive the physical parameters of the clumps, using standard molecular spectral-line analysis techniques. This analysis yields the following range of properties: integrated line intensity 1-30 K km s-1, peak line brightness 1-7 K, linewidth 1-10 km s-1, integrated line luminosity 0.5-200 K km s-1 pc^2, FWHM size 0.2-2.5 pc, mean projected axial ratio 2, optical depth 0.08-2, total surface density 30-3000 M{\sun} pc-2, number density 0.2-30 x 10^9 m-3, mass 15-8000 M{\sun}, virial parameter 1-55, and total gas pressure 0.3-700 pPa. We find that the CHaMP clumps do not obey a Larson-type size-linewidth relation. Among the clumps, there exists a large population of subthermally excited, weakly-emitting (but easily detectable) dense molecular clumps, confirming the prediction of Narayanan et al. (2008). These weakly-emitting clumps comprise 95% of all massive clumps by number, and 87% of the molecular mass, in this portion of the Galaxy; their properties are distinct from the brighter massive star-forming regions that are more typically studied. If the clumps evolve by slow contraction, the 95% of fainter clumps may represent a long-lived stage of pressure-confined, gravitationally stable massive clump evolution, while the CHaMP ... (abridged)

preprint2011arXiv

The Magellanic Mopra Assessment (MAGMA). I. The Molecular Cloud Population of the Large Magellanic Cloud

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

preprint2011arXiv

Toward Understanding the Cosmic-ray Acceleration at Young Supernova Remnants Interacting with Interstellar Clouds: Possible Applications to RX J1713.7-3946

Using three-dimensional magnetohydrodynamics simulations, we investigate general properties of a blast wave shock interacting with interstellar clouds. The shock-cloud interaction generates a turbulent shell through the vorticity generations. In the turbulent shell, the magnetic field is amplified as a result of turbulent dynamo action. In the case of a young supernova remnant, the corresponding strength of the magnetic field is approximately 1 mG. The propagation speed of the shock wave is significantly stalled in the clouds because of the high density, while the shock maintains a high velocity in the diffuse surrounding. In addition, when the shock wave hits the clouds, reflection shock waves are generated that propagate back into the shocked shell. From these simulation results, many observational characteristics of a young SNR RX J1713.7-3946 that is suggested to be interacting with molecular clouds, can be explained as follows: The reflection shocks can accelerate particles in the turbulent downstream region where the magnetic field strength reaches 1mG, which causes short-time variability of synchrotron X-rays. Since the shock velocity is stalled locally in the clouds, the temperature in the shocked cloud is suppressed below 1 keV. Thus, thermal X-ray line emissions would be faint even if the SNR is interacting with molecular clouds. We also find that the photon index of the $π^0$-decay gamma rays generated by cosmic-ray protons can be 1.5, because the penetration depth of high-energy particles into the clumpy clouds depends on their energy, which is consistent with the recent gamma-ray observation by Fermi space telescope.

preprint2010arXiv

A Detailed Observational Study of Molecular Loops 1 and 2 in the Galactic Center

Fukui et al. (2006) discovered two huge molecular loops in the Galactic center located in (l, b) ~ (355 deg-359 deg, 0 deg-2 deg) in a large velocity range of -180-40 km s^-1. Following the discovery, we present detailed observational properties of the two loops based on NANTEN 12CO(J=1-0) and 13CO(J=1-0) datasets at 10 pc resolution including a complete set of velocity channel distributions and comparisons with HI and dust emissions as well as with the other broad molecular features. We find new features on smaller scales in the loops including helical distributions in the loop tops and vertical spurs. The loops have counterparts of the HI gas indicating that the loops include atomic gas. The IRAS far infrared emission is also associated with the loops and was used to derive an X-factor of 0.7(+/-0.1){\times}10^20 cm^-2 (K km s^-1)^-1 to convert the 12CO intensity into the total molecular hydrogen column density. From the 12CO, 13CO, H I and dust datasets we estimated the total mass of loops 1 and 2 to be ~1.4 {\times} 106 Msun and ~1.9 {\times} 10^6 Msun, respectively, where the H I mass corresponds to ~10-20% of the total mass and the total kinetic energy of the two loops to be ~10^52 ergs. An analysis of the kinematics of the loops yields that the loops are rotating at ~47 km s-1 and expanding at ~141 km s^-1 at a radius of 670 pc from the center. Fukui et al. (2006) presented a model that the loops are created by the magnetic flotation due to the Parker instability with an estimated magnetic field strength of ~150 μG. We present comparisons with the recent numerical simulations of the magnetized nuclear disk by Machida et al. (2009) and Takahashi et al. (2009) and show that the theoretical results are in good agreements with the observations. The helical distributions also suggest that some magnetic instability plays a role similarly to the solar helical features.

preprint2010arXiv

Dense Clumps in Giant Molecular Clouds in the Large Magellanic Cloud: Density and Temperature Derived from $^{13}$CO($J=3-2$) Observations

In order to precisely determine temperature and density of molecular gas in the Large Magellanic Cloud, we made observations of optically thin $^{13}$CO($J=3-2$) transition by using the ASTE 10m telescope toward 9 peaks where $^{12}$CO($J=3-2$) clumps were previously detected with the same telescope. The molecular clumps include those in giant molecular cloud (GMC) Types I (with no signs of massive star formation), II (with HII regions only), and III (with HII regions and young star clusters). We detected $^{13}$CO($J=3-2$) emission toward all the peaks and found that their intensities are 3 -- 12 times lower than those of $^{12}$CO($J=3-2$). We determined the intensity ratios of $^{12}$CO($J=3-2$) to $^{13}$CO($J=3-2$), $R^{12/13}_{3-2}$, and $^{13}$CO($J=3-2$) to $^{13}$CO($J=1-0$), $R^{13}_{3-2/1-0}$, at 45$\arcsec$ resolution. These ratios were used for radiative transfer calculations in order to estimate temperature and density of the clumps. The parameters of these clumps range kinetic temperature $T\mathrm{_{kin}}$ = 15 -- 200 K, and molecular hydrogen gas density $n(\mathrm{H_2})$ = 8$\times 10^2$ -- 7$\times 10^3$ cm$^{-3}$. We confirmed that the higher density clumps show higher kinetic temperature and that the lower density clumps lower kinetic temperature at a better accuracy than in the previous work. The kinetic temperature and density increase generally from a Type I GMC to a Type III GMC. We interpret that this difference reflects an evolutionary trend of star formation in molecular clumps. The $R^{13}_{3-2/1-0}$ and kinetic temperature of the clumps are well correlated with H$α$ flux, suggesting that the heating of molecular gas $n(\mathrm{H_2})$ = $10^3$ -- $10^4$ cm$^{-3}$ can be explained by stellar FUV photons.

preprint2010arXiv

High Excitation Molecular Gas in the Galactic Center Loops; 12CO(J =2-1 and J =3-2) Observations

We have carried out 12CO(J =2-1) and 12CO(J =3-2) observations at spatial resolutions of 1.0-3.8 pc toward the entirety of loops 1 and 2 and part of loop 3 in the Galactic center with NANTEN2 and ASTE. These new results revealed detailed distributions of the molecular gas and the line intensity ratio of the two transitions, R3-2/2-1. In the three loops, R3-2/2-1 is in a range from 0.1 to 2.5 with a peak at ~ 0.7 while that in the disk molecular gas is in a range from 0.1 to 1.2 with a peak at 0.4. This supports that the loops are more highly excited than the disk molecular gas. An LVG analysis of three transitions, 12CO J =3-2 and 2-1 and 13CO J =2-1, toward six positions in loops 1 and 2 shows density and temperature are in a range 102.2 - 104.7 cm-3 and 15-100 K or higher, respectively. Three regions extended by 50-100 pc in the loops tend to have higher excitation conditions as characterized by R3-2/2-1 greater than 1.2. The highest ratio of 2.5 is found in the most developed foot points between loops 1 and 2. This is interpreted that the foot points indicate strongly shocked conditions as inferred from their large linewidths of 50-100 km s-1, confirming the suggestion by Torii et al. (2010b). The other two regions outside the foot points suggest that the molecular gas is heated up by some additional heating mechanisms possibly including magnetic reconnection. A detailed analysis of four foot points have shown a U shape, an L shape or a mirrored-L shape in the b-v distribution. It is shown that a simple kinematical model which incorporates global rotation and expansion of the loops is able to explain these characteristic shapes.

preprint2010arXiv

Temperature and Density in the Foot Points of the Molecular Loops in the Galactic Center; Analysis of Multi-J Transitions of 12CO(J=1-0, 3-2, 4-3, 7-6), 13CO(J=1-0) and C18O(J=1-0)

Fukui et al. (2006) discovered two molecular loops in the Galactic center and argued that the foot points of the molecular loops, two bright spots at both loops ends, represent the gas accumulated by the falling motion along the loops, subsequent to magnetic flotation by the Parker instability. We have carried out sensitive CO observations of the foot points toward l=356 deg at a few pc resolution in the six rotational transitions of CO; 12CO(J=1-0, 3-2, 4-3, 7-6), 13CO(J=1-0) and C18O(J=1-0). The high resolution image of 12CO (J=3-2) has revealed the detailed distribution of the high excitation gas including U shapes, the outer boundary of which shows sharp intensity jumps accompanying strong velocity gradients. An analysis of the multi-J CO transitions shows that the temperature is in a range from 30-100 K and density is around 10^3-10^4 cm^-3, confirming that the foot points have high temperature and density although there is no prominent radiative heating source such as high mass stars in or around the loops. We argue that the high temperature is likely due to the shock heating under C-shock condition caused by the magnetic flotation. We made a comparison of the gas distribution with theoretical numerical simulations and note that the U shape is consistent with numerical simulations. We also find that the region of highest temperature of ~100 K or higher inside the U shape corresponds to the spur having an upward flow, additionally heated up either by magnetic reconnection or bouncing in the interaction with the narrow neck at the bottom of the U shape. We note these new findings further reinforce the magnetic floatation interpretation.

preprint2009arXiv

A peculiar jet and arc of molecular gas toward the rich and young stellar cluster Westerlund 2 and a TeV gamma ray source

We have discovered remarkable jet- and arc-like molecular features toward the rich and young stellar cluster Westerlund2. The jet has a length of ~100 pc and a width of ~10 pc, while the arc shows a crescent shape with a radius of ~30 pc. These molecular features each have masses of ~10000 solar mass and show spatial correlations with the surrounding lower density HI gas. The jet also shows an intriguing positional alignment with the core of the TeV gamma ray source HESS J1023-575 and with the MeV/GeV gamma-ray source recently reported by the Fermi collaboration. We argue that the jet and arc are caused by an energetic event in Westerlund 2, presumably due to an anisotropic supernova explosion of one of the most massive member stars. While the origin of the TeV and GeV gamma-ray sources is uncertain, one may speculate that they are related to the same event via relativistic particle acceleration by strong shock waves produced at the explosion or by remnant objects such as a pulsar wind nebula or microquasar.

preprint2009arXiv

Discovery of Large-Scale Gravitational Infall in a Massive Protostellar Cluster

We report Mopra (ATNF), Anglo-Australian Telescope, and Atacama Submillimeter Telescope Experiment observations of a molecular clump in Carina, BYF73 = G286.21+0.17, which give evidence of large-scale gravitational infall in the dense gas. From the millimetre and far-infrared data, the clump has mass ~ 2 x 10^4 Msun, luminosity ~ 2-3 x 10^4 Lsun, and diameter ~ 0.9 pc. From radiative transfer modelling, we derive a mass infall rate ~ 3.4 x 10^-2 Msun yr-1. If confirmed, this rate for gravitational infall in a molecular core or clump may be the highest yet seen. The near-infrared K-band imaging shows an adjacent compact HII region and IR cluster surrounded by a shell-like photodissociation region showing H2 emission. At the molecular infall peak, the K imaging also reveals a deeply embedded group of stars with associated H2 emission. The combination of these features is very unusual and we suggest they indicate the ongoing formation of a massive star cluster. We discuss the implications of these data for competing theories of massive star formation.

preprint2009arXiv

Discovery of Molecular Loop 3 in the Galactic Center: Evidence for a Positive-Velocity Magnetically Floated Loop towards $L=355^\circ-359^\circ$

We have discovered a molecular dome-like feature towards $355^{\circ} \leq l \leq 359^{\circ}$ and $0^{\circ} \leq b \leq 2^{\circ}$. The large velocity dispersions of 50--100 km s$^{-1}$ of this feature are much larger than those in the Galactic disk and indicate that the feature is located in the Galactic center, probably within $\sim1$ kpc of Sgr A$^{*}$. The distribution has a projected length of $\sim600$ pc and height of $\sim300$ pc from the Galactic disk and shows a large-scale monotonic velocity gradient of $\sim130$ km s $^{-1}$ per $\sim600$ pc. The feature is also associated with HI gas having a more continuous spatial and velocity distribution than that of $^{12}$CO. We interpret the feature as a magnetically floated loop similar to loops 1 and 2 and name it "loop 3". Loop 3 is similar to loops 1 and 2 in its height and length but is different from loops 1 and 2 in that the inner part of loop 3 is filled with molecular emission. We have identified two foot points at the both ends of loop 3. HI, $^{12}$CO and $^{13}$CO datasets were used to estimate the total mass and kinetic energy of loop 3 to be $\sim3.0 \times 10^{6} \Mo$ and $\sim1.7 \times 10^{52}$ ergs. The huge size, velocity dispersions and energy are consistent with the magnetic origin the Parker instability as in case of loops 1 and 2 but is difficult to be explained by multiple stellar explosions. We argue that loop 3 is in an earlier evolutionary phase than loops 1 and 2 based on the inner-filled morphology and the relative weakness of the foot points. This discovery indicates that the western part of the nuclear gas disk of $\sim1$ kpc radius is dominated by the three well-developed magnetically floated loops and suggests that the dynamics of the nuclear gas disk is strongly affected by the magnetic instabilities.

preprint2009arXiv

Formation of Galactic Center Magnetic Loops

A survey for the molecular clouds in the Galaxy with NANTEN mm telescope has discovered molecular loops in the Galactic center region. The loops show monotonic gradients of the line of sight velocity along the loops and the large velocity dispersions towards their foot points. It is suggested that these loops are explained in terms of the buoyant rise of magnetic loops due to the Parker instability. We have carried out global three-dimensional magneto-hydrodynamic simulations of the gas disk in the Galactic center. The gravitational potential is approximated by the axisymmetric potential proposed by Miyamoto & Nagai (1975). At the initial state, we assume a warm (~ 10^4 K) gas torus threaded by azimuthal magnetic fields. Self-gravity and radiative cooling of the gas are ignored. We found that buoyantly rising magnetic loops are formed above the differentially rotating, magnetically turbulent disk. By analyzing the results of global MHD simulations, we have identified individual loops, about 180 in the upper half of the disk, and studied their statistical properties such as their length, width, height, and velocity distributions along the loops. Typical length and height of a loop are 1kpc and 200pc, respectively. The line of sight velocity changes linearly along a loop and shows large dispersions around the foot-points. Numerical results indicate that loops emerge preferentially from the region where magnetic pressure is large. We argue that these properties are consistent with those of the molecular loops discovered by NANTEN.

preprint2009arXiv

Search for Remnant Clouds Associated with the TW Hya Association

We report on a search for the parental molecular clouds of the TW Hya association (TWA), using CO emission and Na I absorption lines. TWA is the nearest young (~ 50 pc; ~ 10 Myr) stellar association, yet in spite of its youth, there are no detection of any associated natal molecular gas, as is the case for other typical young clusters. Using infrared maps as a guide, we conducted a CO cloud survey toward a region with a dust extinction of E(B-V) > 0.2 mag, or AV > 0.6 mag. CO emission is detected toward three IR dust clouds, and we reject one cloud from the TWA, as no interstellar Na absorption was detected from the nearby Hipparcos stars, implying that it is too distant to be related. The other two clouds exhibit only faint and small-scale CO emission. Interstellar Na I absorptions of Hipparcos targets, HIP 57809, HIP 64837, and HIP 64925 (at distances of 133, 81, and 101 pc, respectively) by these couds is also detected. We conclude that only a small fraction of the interstellar matter (ISM) toward the IR dust cloud is located at distance less than 100 pc, which may be all that is left of the remnant clouds of TWA; the remaining remnant cloud having dissipated in the last ~ 1 Myr. Such a short dissipation timescale may be due to an external perturbation or kinematic segregation that has a large stellar proper motion relative to the natal cloud.

preprint2009arXiv

Similarity between the Molecular Loops in the Galactic Center and the Solar Chromospheric Arch Filaments

We carried out two-dimensional magnetohydrodynamic simulations of the Galactic gas disk to show that the dense loop-like structures discovered by the Galactic center molecular cloud survey by NANTEN 4 m telescope can be formed by the buoyant rise of magnetic loops due to the Parker instability. At the initial state, we assumed a gravitationally stratified disk consisting of the cool layer ($T \sim 10^3$ K), warm layer ($T \sim 10^4$ K), and hot layer ($T \sim 10^5$ K). Simulation box is a local part of the disk containing the equatorial plane. The gravitational field is approximated by that of a point mass at the galactic center. The self-gravity, and the effects of the galactic rotation are ignored. Numerical results indicate that the length of the magnetic loops emerging from the disk is determined by the scale height of the hot layer ($\sim$ 100 pc at 1 kpc from the Galactic center). The loop length, velocity gradient along the loops and large velocity dispersions at their foot points are consistent with the NANTEN observations. We also show that the loops become top-heavy when the curvature of the loop is sufficiently small, so that the rising loop accumulates the overlying gas faster than sliding it down along the loop. This mechanism is similar to that of the formation of solar chromospheric arch filaments. The molecular loops emerge from the low temperature layer just like the dark filaments observed in the H$α$ image of the emerging flux region of the sun.

preprint2009arXiv

Star Formation Around Supergiant Shells in the LMC

We examine the recent star formation associated with four supergiant shells (SGSs) in the Large Magellanic Cloud (LMC): LMC 1, 4, 5, and 6, which have been shown to have simple expanding-shell structures. H II regions and OB associations are used to infer star formation in the last few Myr, while massive young stellar objects (YSOs) reveal the current ongoing star formation. Distributions of ionized, H I, and molecular components of the interstellar gas are compared with the sites of recent and current star formation to determine whether triggering has taken place. We find that a great majority of the current star formation has occurred in gravitationally unstable regions, and that evidence of triggered star formation is prevalent at both large and local scales.

preprint2009arXiv

Temperature and Density Distribution in the Molecular Gas Toward Westerlund 2: Further Evidence for Physical Association

Furukawa et al. 2009 reported the existence of a large mass of molecular gas associated with the super star cluster Westerlund 2 and the surrounding HII region RCW49, based on a strong morphological correspondence between NANTEN2 12CO(J=2-1) emission and Spitzer IRAC images of the HII region. We here present temperature and density distributions in the associated molecular gas at 3.5 pc resolution, as derived from an LVG analysis of the 12CO(J=2-1), 12CO(J=1-0) and 13CO(J=2-1) transitions. The kinetic temperature is as high as 60-150 K within a projected distance of 5-10 pc from Westerlund 2 and decreases to as low as 10 K away from the cluster. The high temperature provides robust verification that the molecular gas is indeed physically associated with the HII region, supporting Furukawa et al.'s conclusion. The derived temperature is also roughly consistent with theoretical calculations of photo dissociation regions (PDRs), while the low spatial resolution of the present study does not warrant a more detailed comparison with PDR models. We suggest that the molecular clouds presented here will serve as an ideal laboratory to test theories on PDRs in future higher resolution studies.

preprint2009arXiv

Warm and Dense Molecular Gas in the N159 Region: 12CO J=4-3 and 13CO J=3-2 Observations with NANTEN2 and ASTE

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

preprint2008arXiv

The Molecular Ridge Close to 30 Doradus in the Large Magellanic Cloud

With the ATNF Mopra telescope we are performing a survey in the 12CO(1-0) line to map the molecular gas in the Large Magellanic Cloud (LMC). For some regions we also obtained interferometric maps of the high density gas tracers HCO+ and HCN with the Australia Telescope Compact Array (ATCA). Here we discuss the properties of the elongated molecular complex that stretches about 2 kpc southward from 30 Doradus. Our data suggests that the complex, which we refer to as the ``molecular ridge,'' is not a coherent feature but consists of many smaller clumps that share the same formation history. Likely molecular cloud formation triggers are shocks and shearing forces that are present in the surrounding south-eastern HI overdensity region, a region influenced by strong ram pressure and tidal forces. The molecular ridge is at the western edge of the the overdensity region where a bifurcated velocity structure transitions into a single disk velocity component. We find that the 12CO(1-0) and HI emission peaks in the molecular ridge are typically near each other but never coincide. A likely explanation is the conversion of warmer, low-opacity HI to colder, high-opacity HI from which H2 subsequently forms. On smaller scales, we find that very dense molecular gas, as traced by interferometric HCO+ and HCN maps, is associated with star formation along shocked filaments and with rims of expanding shell-like structures, both created by feedback from massive stars.

preprint2008arXiv

The Second Survey of the Molecular Clouds in the Large Magellanic Cloud by NANTEN I: Catalog of Molecular Clouds

The second survey of the molecular clouds in 12CO (J = 1-0) was carried out in the Large Magellanic Cloud by NANTEN. The sensitivity of this survey is twice as high as that of the previous NANTEN survey, leading to a detection of molecular clouds with M_CO > 2 x 10^4 M_sun. We identified 272 molecular clouds, 230 of which are detected at three or more observed positions. We derived the physical properties, such as size, line width, virial mass, of the 164 GMCs which have an extent more than the beam size of NANTEN in both the major and minor axes. The CO luminosity and virial mass of the clouds show a good correlation of M_VIR propto L_CO^{1.1 +- 0.1} with a Spearman rank correlation of 0.8 suggesting that the clouds are in nearly virial equilibrium. Assuming the clouds are in virial equilibrium, we derived an X_CO-factor to be ~ 7 x 10^20 cm^-2 (K km s^-1)^-1. The mass spectrum of the clouds is fitted well by a power law of N_cloud(>M_CO) proportional to M_CO^{-0.75 +- 0.06} above the completeness limit of 5 x 10^4 M_sun. The slope of the mass spectrum becomes steeper if we fit only the massive clouds; e.g., N_cloud (>M_CO) is proportional to M_CO^{-1.2 +- 0.2} for M_CO > 3 x 10^5 M_sun.

preprint2006arXiv

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy's Evolution (SAGE) I: Overview and Initial Results

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

preprint2000arXiv

Molecular Clouds Around a Run-away O Star, zeta Oph

Molecular clouds around a run-away O star z Oph have been surveyed with NANTEN telescope and their streaming motion caused by z Oph has been detected. z Oph is the earliest member of the Sco OB2 association and is a runaway star rapidly moving with an HII region, S27. We detected 2 major filamentary cloud complexes; one including L156 (L156 complex) is lying across nearly the center of S27, and the other one (L204 complex) is located near the eastern edge of S27. Total masses of them traced by the 12CO are 520 Mo and 1110 Mo, respectively. Denser cloud cores detected in C18O are locally distributed on the near side of the L204 complex to z Oph, and lower density gas traced by 12CO spreads toward the opposite side. Both complexes have radial velocity shifts that are correlated with the gas density. These spatial and velocity structures can be interpreted as follows; (1) the L156 complex is stuck on the expanding Stromgrem sphere and has been accelerated, (2) the gas in L204 complex was compressed and has also been accelerated outward by z Oph, resulting the radial velocity shifts of diffuse low-density gas relative to the dense cores embedded in the clouds. These density and velocity structures indicate dynamical interaction between the HII region and the molecular clouds. The cloud complexes are divided into 7 clouds by intensity distributions. In order to investigate the acceleration mechanism, we calculated momentum and kinetic energy for each cloud. They range from 60 to 800 Mo km/s and from 0.9 to 21 times 10^45 erg, respectively. We examined the effects of the stellar wind and photo evaporation by UV field of z Oph and found that the stellar wind can hardly input the momentum during the crossing time of the movement of z Oph. UV radiation seems to be a more likely origin of the streaming gas motion.