Source author record

Yoshiaki Sofue

Yoshiaki Sofue 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

65works
4topics
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

65 published item(s)

preprint2026arXiv

Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds

The central molecular zone (CMZ), surrounding the Galactic centre, is the largest reservoir of dense molecular gas in the Galaxy. Despite its relative proximity, the 3D structure of the CMZ remains poorly constrained, primarily due to projection effects. We aim to constrain the line-of-sight location of two molecular clouds in the CMZ -- the 50 and 20 km/s clouds -- and to investigate their possible physical connection using stellar kinematics and photometry. This study serves as a pilot for future applications across the full CMZ. We estimated the line-of-sight position of the clouds by analysing stellar kinematics, stellar densities, and stellar populations towards the cloud regions and a control field. We find an absence of westward moving stars in the cloud regions, which indicates that they lie on the near side of the CMZ. This interpretation is supported by the stellar density distributions. The similar behaviour observed in the two clouds, as well as in the region between them (the ridge), suggests that they are located at comparable distances and are physically linked. We also identified an intermediate-age stellar population (2-7 Gyr) in both regions, consistent with that observed on the near side of the CMZ. We estimated the line-of-sight distances at which the clouds and the ridge become kinematically detectable (i.e. where the proper motion component parallel to the Galactic plane differs from that of the control field at the 3 sigma level) by converting their measured proper motions parallel to the Galactic plane using a theoretical model of the stellar distribution. We find that the 50 and 20 km/s clouds are located at $43\pm8$ pc and $56\pm11$ pc from Sgr A*, respectively, and that the ridge lies at $56\pm11$ pc; this supports the idea that the clouds are physically connected through the ridge.

preprint2025arXiv

Rotation and stability of the circumnuclear gas disk in the Galactic Center potential by the ALMA CMZ Exploration Survey (ACES)

We investigated the gravitational potential and mass distribution in the Galactic Center by examining the morphology and kinematics of the circumnuclear gaseous disk revealed by the molecular line data from the ALMA CMZ Exploration Survey (ACES). We obtain an estimate of the shape of the potential {within the central $\sim 20$ pc} to reproduce the observed properties of the circumnuclear gas disk (CND) by simulating the motion of test particles for various axial ratios and show that the potential is approximately spherical. We construct a rotation curve by applying the terminal velocity method to the position-velocity diagrams, and calculate the mass distribution in the Galactic Center. The distribution of mass density is found to be of cusp type, approximated by $ρ_{\rm mass} \sim 1.56\times 10^5(R/1 {\rm pc})^{-1.9}~M_{\odot} {\rm pc}^{-3}$, where $R$ is the distance from the nucleus. We discuss the tidal effect caused by the gravitational potential that produces the rotation curve and show that the gas disk is stable against self-gravitational contraction within a critical radius of $ R_{\rm T}\sim 14 ~(ρ_{\rm gas}/10^5 {\rm H_2~cm^{-3}})^{-1/2}~{\rm pc}$. This suggests suppression of star formation and a top-heavy IMF in the circmunuclear region.

preprint2022arXiv

Discovery of non-equilibrium ionization plasma associated with the North Polar Spur and Loop I

We investigated the detailed plasma condition of the North Polar Spur (NPS)/Loop I using archival $Suzaku$ data. In previous research collisional ionization equilibrium (CIE) have been assumed for X-ray plasma state, but we also assume non-equilibrium ionization (NEI) to check the plasma condition in more detail. We found that most of the plasma in the NPS/Loop I favors the state of NEI, and has the density-weighted ionization timescale of $n_e t\sim10^{11-12}$ s cm$^{-3}$ and the electron number density $n_e\sim$ a few $\times$ 10$^{-3}$ cm$^{-3}$. The plasma shock age, $t$, or the time elapsed after the shock front passed through the plasma, is estimated to be on the order of a few $\rm{Myr}$ for the NPS/Loop I, which puts a strict lower limit to the age of the whole NPS/Loop I structure. We found that NEI results in significantly higher temperature and lower emission measure than those currently derived under CIE assumption. The electron temperature under NEI is estimated to be as high as 0.5 keV toward the brightest X-ray NPS ridge at $Δθ=-20^\circ$, which decreases to 0.3 keV at $-10^\circ$, and again increases to $\sim 0.6$ keV towards the outer edge of Loop I at $Δθ\sim0^\circ$, about twice the currently estimated temperatures. Here, $Δθ$ is the angular distance from the outer edge of Loop I. We discuss the implication of introducing NEI for the research in plasma states in astrophysical phenomena.

preprint2022arXiv

Molecular Fraction in the Galactic Center: The Central Molecular and HI Zones

By mapping the molecular fraction of the Galactic Center (GC), we quantitatively address the question of how much molecular and central the CMZ (Central Molecular Zone) is. For this purpose we analyze the CO and HI-line archival data, and determine the column- (surface-) and volume-molecular fractions, $f_{\rm mol}^Σ$ and $f_{\rm mol}^ρ$, which are the ratio of column-mass density of H$_2$ projected on the sky to that of total gas (H$_2$ + HI) from the line intensities, and the ratio of volume-mass densities of \Htwo\ to total gas from the brightness temperature, respectively. It is shown that $f_{\rm mol}^Σ$ is as high as $ \sim 0.9-0.95$ in the CMZ, and $f_{\rm mol}^ρ$ is $0.93-0.98$ in the GC Arms I and II attaining the highest value of $\sim 0.98$ toward Sgr B2. The expanding molecular ring (EMR, or the parallelogram) has a slightly smaller $f_{\rm mol}^Σ$ as $\sim 0.9-0.93$. We define the CMZ as the region with $f_{\rm mol}^Σ\ge 0.8-0.9$ between the shoulders of plateau-like distribution of H$_2$ column density from $l=-1^\circ.1$ to $+1^\circ.8$ having Gaussian vertical distribution with a half thickness of $\pm 0^\circ.2$. The CMZ is embedded in the Central HI Zone (CHZ), which is defined as an HI disc between $l\sim -2^\circ$ and $+2^\circ.5$, $b=-0^\circ.5$ and $+0^\circ.5$. Based on the analysis, we discuss the origin of CMZ and interstellar physics such as the volume filling factors of molecular and HI gases inferred from the difference between $ f_{\rm mol}^Σ$ and $f_{\rm mol}^ρ$.

preprint2022arXiv

One-kpc Expanding Cylinder of HI-Gas around the Galactic Center

We report the discovery of an expanding cylinder of HI gas of radius 1 kpc and vertical extent 800 pc by analyzing the 21-cm line survey data from the literature. The cylinder is expanding at 150 km s$^{-1}$ and rotating at 100 km s$^{-1}$, and is interpreted as due to a high-velocity conical wind at $\sim 180$ km s$^{-1}$ from the Galactic Center. The total mass of the cylinder is estimated to be $\sim 8.5\times 10^5 M_\odot$ and kinetic energy $\sim 3\times 10^{53}$ ergs.

preprint2022arXiv

Three-Dimensional Structure of the Central Molecular Zone

A detailed comparison of HI and CO line cube data of the Galactic Center (GC) region from the archives is obtained. The central molecular zone (CMZ) is shown to be embedded in the HI disc (central HI zone, CHZ) of radius $\sim 320$ pc and vertical scale height $\sim 70$ pc. A radio continuum belt is shown to run parallel to molecular Arms I and II. The belt draws a double-infinity shape on the sky, connecting Sgr E ($l\sim -1^\circ.2)$, C, B1, B2 and Sgr D ($+1^\circ.2$), and is interpreted as a warping star-forming ring. The molecular Arms are closely associated with the HI arms on the longitude-velocity diagram (LVD), showing coherent rigid-body ridges. Due to the close relationship between HI and CO, the HI line absorption can be used to determine the Arms' position relative to Sgr A, B1, B2 and C.Combining the trigonometric data of proper motions of Sgr A$^*$ and maser sources of Sgr B2 as well as radial velocities, the 3D velocity vector of Sgr B2 is determined. From these analyses, the molecular Arm I with Sgr B2 is shown to be located in the near side of Sgr A$^*$, and Arm II with Sgr C in the other side, both composing a pair of symmetrical Arms around the GC. We present a possible 3D view of Sgr A through E and Arms I and II along with a parameter list.

preprint2021arXiv

Dark Supernova Remnants revealed by CO-line Bubbles in the W43 Molecular Complex along the 4-kpc Arm

Fine structure of the density distribution in giant molecular clouds (GMC) around W43 (G31+00+90 km/s at $\sim 5.5$ kpc) was analyzed using the FUGIN$^*$ CO-line survey at high-angular ($20''\sim 0.5$ pc) and velocity (1.3 km/s) resolutions ($^*$Four-receiver-system Unbiased Galactic Imaging survey with the Nobeyama 45-m telescope). The GMCs show highly turbulent structures, and the eddies are found to exhibit spherical bubble morphology appearing in narrow ranges of velocity channels. The bubbles are dark in radio continuum emission, unlike usual supernova remnants (SNR) or HII regions, and in infrared dust emission, unlike molecular bubbles around young stellar objects. The CO bubbles are interpreted as due to fully evolved buried SNRs in molecular clouds after rapid exhaustion of the released energy in dense molecular clouds. The CO bubbles may be a direct evidence for exciting and maintaining the turbulence in GMCs by SN origin. Search for CO bubbles as "dark SNRs" (dSNR) will have implication to estimate the supernova rate more accurately, and hence the star formation activity in the Milky Way.

preprint2021arXiv

Origin of Galactic Spurs: New Insight from Radio/X-ray All-sky Maps

In this study, we analyze giant Galactic spurs seen in both radio and X-ray all-sky maps to reveal their origins. We discuss two types of giant spurs: one is the brightest diffuse emission near the map's center, which is likely to be related to Fermi bubbles (NPSs/SPSs, north/south polar spurs, respectively), and the other is weaker spurs that coincide positionally with local spiral arms in our Galaxy (LAS, local arm spur). Our analysis finds that the X-ray emissions, not only from the NPS but from the SPS are closer to the Galactic center by ~5 deg compared with the corresponding radio emission. Furthermore, larger offsets of 10-20 deg are observed in the LASs; however, they are attributed to different physical origins. Moreover, the temperature of the X-ray emission is kT ~ 0.2 keV for the LAS, which is systematically lower than those of the NPS and SPS (kT ~ 0.3 keV) but consistent with the typical temperature of Galactic halo gas. We argue that the radio/X-ray offset and the slightly higher temperature of the NPS/SPS X-ray gas are due to the shock compression/heating of halo gas during a significant Galactic explosion in the past, whereas the enhanced X-ray emission from the LAS may be due to the weak condensation of halo gas in the arm potential or star formation activity without shock heating.

preprint2020arXiv

CO-Line and Radio Continuum Study of Elephant Trunks: The Pillars of Creation in M16

Molecular-line and radio continuum properties of the elephant trunks (ET, pillars of creation) in M16 are investigated by analyzing the 12CO(J=1-0), 13CO(J=1-0) and C18O(J=1-0)-line survey data with the Nobeyama 45-m telescope and the Galactic plane radio survey at 20 and 90 cm with the Very Large Array. The head clump of Pillar West I is found to be the brightest radio source in M16, showing thermal spectrum and property of a compact HII region with the nearest O5 star in NGC 6611 being the heating source. The radio pillars have cometary structure concave to the molecular trunk head, and the surface brightness distribution obeys a simple illumination law by a remote excitation source. The molecular density in the pillar head is estimated to be several $10^4$ H$_2$ cm$^{-3}$ and molecular mass $\sim 13-40 M\odot$. CO-line kinematics reveals random rotation of the clumps in the pillar tail at $\sim 1-2$ km/s, comparable to the velocity dispersion and estimated Alfven velocity. It is suggested that the random directions of velocity gradients would manifest a torsional magnetic oscillation of the clumps around the pillar axis.

preprint2020arXiv

CO-to-H$_2$ Conversion and Spectral Column Density in Molecular Clouds: The Vriability of $X_{\rm CO}$ Factor

Analyzing the Galactic plane CO survey with the Nobeyama 45-m telescope, we compared the spectral column density (SCD) of H$_2$ calculated for $^{12}$CO line using the current conversion factor $X_{\rm CO}$ to that for $^{13}$CO line under LTE in M16 and W43 regions. Here, SCD is defined by $dN_{\rm H_2}/dv$ with $N_{\rm H_2}$ and $v$ being the column density and radial velocity, respectively. It is found that the $X_{\rm CO}$ method significantly under-estimates the H$_2$ density in a cloud or region, where SCD exceeds a critical value ($ \sim 3\times 10^{21}\ [{\rm H_2 \ cm^{-2} \ (km \ s^{-1})^{-1}}]$), but over-estimates in lower SCD regions. We point out that the actual CO-to-H$_2$ conversion factor varies with the H$_2$ column density or with the CO-line intensity: It increases in the inner and opaque parts of molecular clouds, whereas it decreases in the low-density envelopes. However, in so far as the current $X_{\rm CO}$ is used combined with the integrated $^{12}$CO intensity averaged over an entire cloud, it yields a consistent value with that calculated using the $^{13}$CO intensity by LTE. Based on the analysis, we propose a new CO-to-\Htwo conversion relation, $N_{\rm H_2}^* = \int X_{\rm CO}^*(T_{\rm B}) T_{\rm B} dv$, where $X_{\rm CO}^*=(T_{\rm B}/T_{\rm B}^*)^βX_{\rm CO}$ is the modified spectral conversion factor as a function of the brightness temperature, $T_{\rm B}$, of the ${12}$CO ($J=1-0$) line, and $β\sim 1-2$ and $T_{\rm B}^*=12-16$ K are empirical constants obtained by fitting to the observed data. The formula corrects for the over/under estimation of the column density at low/high-CO line intensities, and is applicable to molecular clouds with $T_{\rm B} \ge 1$ K (rms noise in the data) from envelope to cores at sub-parsec scales (resolution). (Full resolution copy available at http://www.ioa.s.u-tokyo.ac.jp/~sofue/news/2020_mn_Xco12co13_fugin.pdf)

preprint2020arXiv

Feedback between Sgr A and B : AGN-Starburst Connection in the Galactic Centre

Propagation of fast-mode magneto-hydrodynamic (MHD) compression waves is traced in the Galactic Centre. MHD waves produced by the active Galactic nucleus (Sgr A) focus on the molecular clouds such as Sgr B in the central molecular zone, which will trigger star formation, or possibly starburst. MHD waves newly excited by the starburst propagate backward, and focus on the nucleus (Sgr A), where implosive waves compress the nuclear gas to promote fueling the nucleus and may trigger nucleus activity. Echoing focusing of MHD waves between Sgr A (active galactic nucleus: AGN) and Sgr B (starburst) trigger each other at high efficiency by minimal energy requirement. It also solves the problem of angular momentum for AGN fueling, as the focusing waves do not require global gas flow.

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

Galactic Centre Threads as Nuclear MHD Waves

Propagation of fast-mode magnetohydrodynamic (MHD) compression waves is traced in the Galactic Center with a poloidal magnetic cylinder. MHD waves ejected from the nucleus are reflected and guided along the magnetic field, exhibiting vertically stretched fronts. The radio threads and non-thermal filaments are explained as due to tangential views of the waves driven by sporadic activity in Sgr A$^*$, or by multiple supernovae. In the latter case, the threads could be extremely deformed relics of old SNRs exploded in the nucleus.

preprint2020arXiv

Giant Elephant Trunks from Giant Molecular Clouds

We report the discovery of large elephant trunk (ET)-like objects, named giant elephant trunk (GET), of molecular gas in star forming complexes in the Scutum and Norma arms using the $^{12}$CO(J=1-0)-line survey data with the Nobeyama 45-m telescope. In comparison with the CO maps of ETs in M16 as derived from the same data, we discuss physical properties of the GETs. Their lengths are $\sim 20$ to 50 pc, an order of magnitude larger than ETs. GETs show a cometary structure coherently aligned parallel to the galactic plane, and emerge from bow-shaped concave surface of giant molecular clouds (GMC) facing the HII regions, and point down-stream of the gas flow in the spiral arms. The molecular masses of the head clumps are $\sim 10^3 -10^4 M_\odot$, about 3 to 4 times the virial masses, indicating that the clumps are gravitationally stable. Jeans masses calculated for the derived density and assumed kinetic temperature are commonly sub-solar. We suggest that the GET heads are possible birth sites of stellar clusters, similarly to ETs' globules, but at much greater scale. We discuss the origin of the GETs by Rayleigh-Taylor instability due to deceleration of GMCs by low density gas stagnated in the galactic shock waves as well as by pressure of the HII regions.

preprint2020arXiv

Gravitational Focusing of Low-Velocity Dark Matter on the Earth's Surface

We show that the Earth acts as a high-efficiency gravitational collector of low-velocity flow of dark matter (DM). The focal point appears on the Earth's surface, when the DM flow speed is about 17 km/s with respect to the geo-center. We discuss diurnal modulation of the local DM density influenced by the Earth's gravity. We also touch upon similar effects on galactic and solar system objects.

preprint2020arXiv

Rotation Curve of the Milky Way and the Dark Matter Density

We review the~current status of the~study of rotation curve (RC) of the~Milky Way, and~present a~unified RC from the~Galactic Center to the galacto-centric distance of about 100 kpc. The~RC is used to directly calculate the~distribution of the~surface mass density (SMD). We then propose a~method to derive the~distribution of dark matter (DM) density in the~in the~Milky Way using the~SMD distribution. The~best-fit dark halo profile yielded a local DM density of $ρ_\odot = 0.36\pm 0.02$ GeV/cc. We also review the~estimations of the~local DM density in the~last decade, and~show that the~value is converging to a~value at $ρ_\odot=0.39\pm 0.09$ GeV/cc.

preprint2019arXiv

Focusing MHD Wave as a Trigger of Star Formation

Propagation of fast-mode magnetohydrodynamic (MHD) waves in the interstellar space is simulated, and focusing MHD wave (FMW) model is proposed for triggered star formation (SF). Waves from an SF region are trapped by nearby molecular clouds and converge onto their focal points, causing implosive compression. Even an isolated cloud suffers from long-distance invasion of waves from remote sources. Echoing SF occurs inside a cloud as well as between clouds. Repetitive refocusing in a filamentary cloud suggests spatial periodicity in SF sites along the filament. The model is applied to the SF regions M16 and M17, where MHD waves produced by M16 are shown to converge onto the focal point of nearby GMC and trigger the SF in M17.

preprint2016arXiv

CO-to-H_2 Conversion Factor of Molecular Clouds using X-Ray Shadows

A new method to determine the CO-to-H$_2$ conversion factor $X_{\rm CO}$ using absorption of diffuse X-ray emission by local molecular clouds was developed. It was applied to the Ophiucus (G353+17) and Corona Australis (G359-18) clouds using CO-line and soft X-ray archival data. We obtained a value $X_{\rm CO} =1.85 \pm 0.45 \times 10^{20} {\rm H_2 cm^{-2}/(K~ km~ s^{-1})}$ as the average of least-$chi^2$ fitting results for R4 (0.7 keV) and R5 (0.8 keV) bands. Full resolution pdf available at http://www.ioa.s.u-tokyo.ac.jp/~sofue/htdocs/2016Xco/

preprint2016arXiv

Rotation and Mass in the Milky Way and Spiral Galaxies

[PASJ Review Paper] Rotation curves are the basic tool for deriving the distribution of mass in spiral galaxies. In this review, we describe various methods to measure rotation curves in the Milky Way and spiral galaxies. We then describe two major methods to calculate the mass distribution using the rotation curve. By the direct method, the mass is calculated from rotation velocities without employing mass models. By the decomposition method, the rotation curve is deconvolved into multiple mass components by model fitting assuming a black hole, bulge, exponential disk and dark halo. The decomposition is useful for statistical correlation analyses among the dynamical parameters of the mass components. We also review recent observations and derived results. ( Full resolution copy is available at URL: http://www.ioa.s.u-tokyo.ac.jp/~sofue/htdocs/PASJreview2016/ )

preprint2016arXiv

Three-Dimensional Aquila Rift: Magnetized HI Arch Anchored by Molecular Complex

Three dimensional structure of the Aquila Rift of magnetized neutral gas is investigated by analyzing HI and CO line data. The projected distance on the Galactic plane of the HI arch of the Rift is ~250 pc from the Sun. The HI arch emerges at l~ 30 deg, reaches to altitudes as high as ~ 500 pc above the plane at l~ 350deg, and returns to the disk at l~270 deg. The extent of arch at positive latitudes is ~1 kpc and radius is ~100 pc. The eastern root is associated with the giant molecular cloud complex, which is the main body of the optically defined Aquila Rift. The HI and molecular masses of the Rift are estimated to be M_HI~1.4x10^5Msun and M_H2 ~3x10^5 Msun. Gravitational energies to lift the gases to their heights are E_grav: HI ~1.4x10^51 and E_grav: H_2 ~0.3x10^51 erg, respectively. Magnetic field is aligned along the HI arch of the Rift, and the strength is measured to be B~10 microG using Faraday rotation measures of extragalactic radio sources. The magnetic energy is estimated to be E_mag~1.2x10^51 erg. A possible mechanism of formation of the Aquila Rift is proposed in terms of interstellar magnetic inflation by a sinusoidal Parker instability of wavelength of ~2.5 kpc and amplitude ~500 pc.

preprint2016arXiv

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: IV. 3D Molecular Fraction and Galactic-Scale HI-to-H$_2$ Transition

Three dimensional (3D) distribution of the volume-density molecular fraction, defined by fmol_rho=rho_H2/(rho_HI + rho_H2), is studied in the Milky Way Galaxy. The molecular front appears at galacto-centric distance of R \sim 8 kpc, where the phase transition from atomic to molecular hydrogen occurs suddenly with fmol_rho dropping from \sim 0.8 to 0.2 within a radial interval as narrow as \sim 0.5 kpc. The front in fmol_rho is much sharper than that for surface density molecular fraction. The front also appears in the vertical direction with a full width of the high-fmol_rho disk to be \sim 100 pc. The radial and vertical fmol_rho profiles, particularly the front behaviors, are fitted by theoretical curves calculated using the observed density profile and assumed radiation field and metallicity with exponential gradients. The molecular fraction was found to be enhanced along the Perseus and some other spiral arms. The fmol_rho arms imply that the molecular clouds are produced from HI gas in the spiral arms and are dissociated in the interarm region. We also show that there is a threshold HI density over which the gas is transformed into molecules.

preprint2015arXiv

Dark Halos of M31 and the Milky Way

Grand rotation curves (GRC) within ~400 kpc of M31 and the Milky Way were constructed by combining disk rotation velocities and radial velocities of satellite galaxies and globular clusters. The GRC for the Milky Way was revised using the most recent Solar rotation velocity. The derived GRCs were deconvolved into a de Vaucouleurs bulge, exponential disk, and a dark halo with the Navarro-Frenk-White (NFW) density profile by the least chi-squares fitting. Comparison of the best-fit parameters revealed similarity of the disks and bulges of the two galaxies, whereas the dark halo mass of M31 was found to be twice the Galaxy's. We show that the NFW model may be a realistic approximation to the observed dark halos in these two giant spirals.

preprint2015arXiv

JVLA S and X-band Polarimetry of the Merging Cluster Abell 2256

We report polarimetry results of a merging cluster of galaxies Abell 2256 with Karl G. Jansky Very Large Array (JVLA). We performed new observations with JVLA at S-band (2051-3947 MHz) and X-band (8051-9947 MHz) in the C array configuration, and detected significant polarized emissions from the radio relic, Source A, and Source B in this cluster. We calculated the total magnetic field strengths toward the radio relic using revised equipartition formula, which is 1.8-5.0 microG. With dispersions of Faraday rotation measure, magnetic-field strengths toward Sources A and B are estimated to be 0.63-1.26 microG and 0.11-0.21 microG, respectively. An extremely high degree of linear polarization, as high as ~ 35 %, about a half of the maximum polarization, was detected toward the radio relic, which indicates highly ordered magnetic lines of force over the beam sizes (~ 52 kpc).The fractional polarization of the radio relic decreases from ~ 35 % to ~ 20 % around 3 GHz as the frequency decreases and is nearly constant between 1.37 and 3 GHz. Both analyses with depolarization models and Faraday tomography suggest multiple depolarization components toward the radio relic and imply the existence of turbulent magnetic fields.

preprint2015arXiv

Metal Enrichment in the Fermi Bubbles as a Probe of Their Origin

The Fermi bubbles are gigantic gamma-ray structures in our Galaxy. The physical origin of the bubbles is still under debate. The leading scenarios can be divided into two categories. One is the nuclear star forming activity similar to extragalactic starburst galaxies and the other is the past active galactic nucleus (AGN) like activity of the Galactic center supermassive black hole. In this paper, we propose that metal abundance measurements will provide an important clue to probe their origin. Based on a simple spherically symmetric bubble model, we find that the generated metallicity and abundance pattern of the bubbles' gas strongly depend on assumed star formation or AGN activities. Star formation scenarios predict higher metallicities and abundance ratios of [O/Fe] and [Ne/Fe] than AGN scenarios do because of supernovae ejecta. Furthermore, the resultant abundance depends on the gamma-ray emission process because different mass injection histories are required for the different gamma-ray emission processes due to the acceleration and cooling time scales of non-thermal particles. Future X-ray missions such as ASTRO-H and Athena will give a clue to probe the origin of the bubbles through abundance measurements with their high energy resolution instruments.

preprint2015arXiv

Rotation Curve Decomposition for Size-Mass Relations of Bulge, Disk, and Dark Halo in Spiral Galaxies

Rotation curves of more than one hundred spiral galaxies were compiled from the literature, and deconvolved into bulge, disk, and dark halo using $χ^2$ fitting in order to determine their scale radii and masses. Correlation analyses were obtained of the fitting parameters for galaxies that satisfied selection and accuracy criteria. Size-mass relations indicate that the sizes and masses are positively correlated among different components in such a way that the larger or more massive is the dark halo, the larger or more massive are the disk and bulge. Empirical size-mass relations were obtained for bulge, disk and dark halo by the least-squares fitting. The disk-to-halo mass ratio was found to be systematically greater by a factor of three than that predicted by cosmological simulations combined with photometry. A preliminary mass function for dark halo was obtained, which is represented by the Schechter function followed by a power law.

preprint2015arXiv

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: III. The Total Neutral Gas Disk

We present newly obtained three-dimensional gaseous maps of the Milky Way Galaxy; HI, H$_2$ and total-gas (HI plus H$_2$) maps, which were derived from the HI and $^{12}$CO($J=1$--0) survey data and rotation curves based on the kinematic distance. The HI and H$_2$ face-on maps show that the HI disk is extended to the radius of 15--20 kpc and its outskirt is asymmetric to the Galactic center, while most of the H$_2$ gas is distributed inside the solar circle. The total gas mass within radius 30 kpc amounts to $8.0\times 10^9$ M$_\odot$, 89\% and 11\% of which are HI and H$_2$, {respectively}. The vertical slices show that the outer HI disk is strongly warped and the inner HI and H$_2$ disks are corrugated. The total gas map is advantageous to trace spiral structure from the inner to outer disk. Spiral structures such as the Norma-Cygnus, the Perseus, the Sagittarius-Carina, the Scutum-Crux, and the Orion arms are more clearly traced in the total gas map than ever. All the spiral arms are well explained with logarithmic spiral arms with pitch angle of $11\degree$ -- $15\degree$. The molecular fraction to the total gas is high near the Galactic center and decreases with the Galactocentric distance. The molecular fraction also locally enhanced at the spiral arms compared with the inter-arm regions.

preprint2014arXiv

Rotation Curve Anomaly and Galactic Warp in M51

We revisit the anomaly of rotation curve in the nearly face-on galaxy M51 that shows an apparently faster decrease of rotation velocity than the Keplerian law in the outer disk, further showing apparent counter rotation in the outermost HI disk. We interpret this anomaly as due to warping of the galactic disk, and determined the warping structure of M51's disk using the tilted-ring method, assuming that the intrinsic rotation curve is normal. It is shown that the disk is nearly flat in the inner disk at a constant inclination angle, but the disk suddenly bends at radius 7.5 kpc by about 27$°$. The inclination angle, then, decreases monotonically outward reaching a perfect face-on ring at 18 kpc, beyond which the disk is warped in the opposite sense to the inner disk, resulting in apparent counter rotation.

preprint2014arXiv

The North Polar Spur and Aquila Rift

Soft X-ray intensity at 0.89 keV along the North Polar Spur is shown to follow the extinction law due to the interstellar gas in the Aquila Rift by analyzing the ROSAT archival data, which proves that the NPS is located behind the rift. The Aquila-Serpens molecular clouds, where the X-ray optical depth exceeds unity, are shown to have a mean LSR velocity of v=7.33 +/- 1.94 km/s, corresponding to a kinematic distance of r=0.642 +/- 0.174 kpc. Assuming a shell structure, a lower limit of the distance to NPS is derived to be 1.01 +/- 0.25 kpc, with the shell center being located farther than 1.1 kpc. Based on the distance estimation, we argue that the NPS is a galactic halo object.

preprint2013arXiv

A Grand Rotation Curve and Dark Matter Halo in the Milky Way Galaxy

A grand rotation curve of the Milky Way Galaxy is constructed, which covers a wide range of radius from the Galactic Center to ~1 Mpc, and is deconvolved into bulge, disk and halo components by least-squares fitting. We determined the scale radii and masses of the bulge and disk to be M_b=(1.652 +/-0.083)x10^10 Msun, a_b=(0.522 +/- 0.037) kpc, M_d=(3.41 +/- 0.41)x10^10 Msun and a_d=(3.19 +/- 0.3) kpc. The dark halo was fitted by the Navaro-Frenk-White (NFW) density profile, rho=rho_0/[(R/h) (1+R/h)^2], and the fit yielded h=12.5 +/- 0.9 kpc and rho_0=(1.06 +/- 0.14)x10^-2 Msun pc^-3. The local dark matter density near the Sun at R_0=8 kpc is estimated to be rho_0_sun=(6.12 +/-0.80)x10^-3 Msun pc^-3 = 0.235 +/- 0.030 GeV cm^-3. The total mass inside the gravitational boundary of the Galaxy at R~ 385 kpc, a half distance to M31, is estimated to be M_h ~ (7.03 +/- 1.01)x10^11 Msun. This leads to the stellar baryon fraction of (M_b+M_d)/(M_b+M_d+M_h)=0.072 +/-0.018. Considering expected baryon fraction in the Local Group, we suggest that baryons in the form of hot gas are filling the dark halo with temperature of ~ 10^6 K and emission measure ~ 10^-5 pc cm^-6. Such hot halo gas may share a small fraction of the observed X-ray background emission.

preprint2013arXiv

Rotation Curve and Mass Distribution in the Galactic Center --- From Black Hole to Entire Galaxy ---

Analyzing high-resolution longitude-velocity (LV) diagrams of the Galactic Center observed with the Nobeyama 45-m telescope in the CO and CS line emissions, we obtain a central rotation curve of the Milky Way. We combine it with the data for the outer disk, and construct a logarithmic rotation curve of the entire Galaxy. The new rotation curve covers a wide range of radius from r ~ 1 pc to several hundred kpc without a gap of data points. It links, for the first time, the kinematical characteristics of the Galaxy from the central black hole to the bulge, disk and dark halo. Using this grand rotation curve, we calculate the radial distribution of surface mass density in the entire Galaxy, where the radius and derived mass densities vary over a dynamical range with several orders of magnitudes. We show that the galactic bulge is deconvolved into two components: the inner (core) and main bulges. Both the two bulge components are represented by exponential density profiles, but the de Vaucouleurs law was found to fail in representing the mass profile of the galactic bulge.

preprint2013arXiv

The Mass Distribution and Rotation Curve in the Galaxy

The mass distribution in the Galaxy is determined by dynamical and photometric methods. Rotation curves are the major tool for determining the dynamical mass distribution in the Milky Way and spiral galaxies. The photometric (statistical) method utilizes luminosity profiles from optical and infrared observations, and assumes empirical values of the mass-to-luminosity (M/L) ratio to convert the luminosity to mass. In this chapter the dynamical method is described in detail, and rotation curves and mass distribution in the Milky Way and nearby spiral galaxies are presented. The dynamical method is categorized into two methods: the decomposition method and direct method. The former fits the rotation curve by calculated curve assuming several mass components such as a bulge, disk and halo, and adjust the dynamical parameters of each component. Explanations are given of the mass profiles as the de Vaucouleurs law, exponential disk, and dark halo profiles inferred from numerical simulations. Another method is the direct method, with which the mass distribution can be directly calculated from the data of rotation velocities without employing any mass models. Some results from both methods are presented, and the Galactic structure is discussed in terms of the mass. Rotation curves and mass distributions in external galaxies are also discussed, and the fundamental mass structures are shown to be universal.

preprint2013arXiv

The Paα Luminosity Function of HII Regions in Nearby Galaxies from HST/NICMOS

The HII region luminosity function (LF) is an important tool for deriving the birthrates and mass distribution of OB associations, and is an excellent tracer of the newly formed massive stars and associations. To date, extensive work (predominantly in Hα) has been done from the ground, which is hindered by dust extinction and the severe blending of adjacent (spatially or in projection) HII regions. Reliably measuring the properties of HII regions requires a linear resolution <40 pc, but analyses satisfying this requirement have been done only in a handful of galaxies, so far. As the first space-based work using a galaxy sample, we have selected 12 galaxies from our HST NICMOS Paα survey and studied the luminosity function and size distribution of HII regions both in individual galaxies and cumulatively, using a virtually extinction-free tracer of the ionizing photon rate. The high angular resolution and low sensitivity to diffuse emission of NICMOS also offer an advantage over ground-based imaging by enabling a higher degree of de-blending of the HII regions. We do not confirm the broken power-law LFs found in ground-based studies. Instead, we find that the LFs, both individual and co-added, follow a single power law dN(L)/dlnL ~ L^-1, consistent with the mass function of star clusters in nearby galaxies, and in agreement with the results of the existing analyses with HST data. The individual and co-added size distribution of HII regions are both roughly consistent with dN(D)/dlnD ~ D^-3, but the power-law scaling is probably contaminated by blended regions or complexes.

preprint2011arXiv

Accuracy Diagrams for the Galactic Rotation Curve and Kinematical Distances

We revisit the methods to determine the Galactic rotation curve and kinematical distances from radial velocities and proper motions. We construct "accuracy diagrams" to show the distributions in the galactic plane of expected uncertainties in the derived quantities such as rotation velocities and kinematical distances. We discuss how to optimize the source selection for measurements of kinematical quantities based on the accuracy diagrams.

preprint2010arXiv

Deep CO Observations and the CO-to-H_2 Conversion Factor in DDO 154, a Low Metallicity Dwarf Irregular Galaxy

We present a deep spectroscopic search for CO emission in the dwarf irregular galaxy DDO154, which has an Oxygen abundance of only 1/20 the solar value. The observations were conducted in order to constrain the CO-to-$\mathrm{H_2}$ conversion factor at low metallicity. No CO was detected, however, despite being one of the sensitive observations done towards galaxies of this type. We succeed in putting a strong lower limit on the conversion factor, at least 10 times the Galactic value. Our result supports previous studies which argue for a high conversion factor at low metallicity.

preprint2010arXiv

The Primordial Origin Model of Magnetic Fields in Spiral Galaxies

We propose a primordial-origin model for the composite configurations of global magnetic fields in spiral galaxies. We show that uniform tilted magnetic field wound up into a rotating disk galaxy can evolve into composite magnetic configurations comprising bisymmetric spiral (S=BSS), axisymmetric spiral (A=ASS), plane-reversed spiral (PR), and/or ring (R) fields in the disk, and vertical (V) fields in the center. By MHD simulations we show that these composite galactic fields are indeed created from weak primordial uniform field, and that the different configurations can co-exist in the same galaxy. We show that spiral fields trigger the growth of two-armed gaseous arms. The centrally accumulated vertical fields are twisted and produce jet toward the halo. We find that the more vertical was the initial uniform field, the stronger is the formed magnetic field in the galactic disk.

preprint2009arXiv

Determining Star Formation Timescale and Pattern Speed in Nearby Spiral Galaxies

We present a revised method for simultaneous determination of the pattern speed and star formation timescale of spiral galaxies, its application, and results for CO and Ha images of nearby spiral galaxies. Out of 13 galaxies, we were able to derive the 2 parameters for 5 galaxies. We categorize them as "C" galaxies, and find (1) The corotation radius is close to the edge of the CO data, and is about half of the optical radius for 3 galaxies. (2) The star formation timescale is roughly consistent with the free-fall time of typical molecular clouds, which indicates that the gravitational instability is the dominant mechanism triggering star formation in spiral arms. (3) The timescale is found to be almost independent of surface density of molecular gas, metallicity, or spiral arm strengths. The number of "C" galaxies and the quality of CO data, however, are not enough to confirm these relationships. We also find that 2 other galaxies show no offsets between CO and Ha, although their arms are clearly traced, and categorize them as "N" galaxies. The presence of a bar could account for this feature, since these 2 galaxies are both barred. With one galaxy excluded from our analysis due to its poor rotation curve, offsets of the remaining 5 galaxies are found to be ambiguous. We categorize them as "A" galaxies. The possible reasons for this ambiguity are (1) the density wave is weaker, and/or (2) observational resolution and sensitivity are not enough to detect the spiral arms and their offsets clearly. The former is supported by our finding that the arm strengths of "A" galaxies are slightly weaker than that of "C" galaxies. [abridged]

preprint2009arXiv

NH3 in the Galactic Center is formed in Cool Conditions

It is an open question why the temperature of molecular gas in the Galactic center region is higher than that of dust. To address this problem, we made simultaneous observations in the NH_3 (J,K) = (1,1), (2,2), and (3,3) lines of the central molecular zone (CMZ) using the Kagoshima 6 m telescope. The ortho-to-para ratio of NH_3 molecules in the CMZ is 1.5--3.5 at most observed area. This ratio is higher than the statistical equilibrium value, and suggests that the formation temperature of NH_3 is 11--20 K. This temperature is similar to the dust temperature estimated from the submillimeter and infrared continuum. This result suggests that the NH_3 molecules in the CMZ were produced on dust grains with the currently observed temperature (11--20 K), and they were released into the gas phase by supernova shocks or collisions of dust particles. The discrepancy between warm molecular gas and cold dust can be explained by the transient heating of the interstellar media in the CMZ approximately 10^5 years ago when the NH_3 molecules were released from the dust.

preprint2008arXiv

A Chain of Dark Clouds in Projection Against the Galactic Center

In the J, H, and Ks bands survey of the the Galactic Center region over an area of 2deg x 5deg, we have found many dark clouds, among which a distinguished chain of dark clouds can be identified with a quiescent CO cloud. The distances of the clouds is estimated to be 3.2-4.2 kpc, corresponding to the Norma arm by our new method to determine distance to dark clouds using the cumulative number of stars against J-Ks colors. Adopting these estimated distances, the size is about 70 pc in length and the total mass of the cloud is 6x10^4 M_solar. Three compact HII regions harbor in the cloud, indicating that star forming activities are going on at the cores of the quiescent CO cloud on the spiral arm.

preprint2008arXiv

Pseudo Rotation Curve connecting the Galaxy, Dark Halo, and Local Group

We construct a Galacto-Local Group rotation curve, combining the Galactic rotation curve with a diagram, where galacto-centric radial velocities of outer globular clusters and member galaxies of the Local Group are plotted against their galacto-centric distances. The high-velocity ends of this pseudo rotation curve within a radius R ~ 150 kpc are well traced by a rotation curve calculated for the NFW (Navaro, Frenk, White) and Burkert dark halo models. The NFW model indicates that the Galaxy's mass within 385 kpc, half distance to M31, is ~ 4 x 10^{11}M_{sun}. High-velocity ends of the pseudo rotation curve for the entire Local Group up to 1.5 Mpc indicate isothermal nature with a terminal velocity of ~ 200 km/s. In order for the Local Group to be gravitationally bound, an order of magnitude larger mass than those of the Galaxy and M31 is required. This fact indicates that the Local Group contains dark matter of mass ~ 5 x 10^{12}M_{sun} filling the space between the Galaxy and M31. The mass density of the Galactic dark halo becomes equal to that of the Local Group's dark matter at R ~ 100 kpc, beyond which the intracluster dark matter dominates. If we define the Galaxy's radius at this distance, the enclosed Galactic mass is ~ 3 x 10^{11}M_{sun}. We comment on the barionic matter fraction of galaxies in the Local Group.

preprint2007arXiv

A Revised Estimate of the CO (J=1-0) Emission from the Host Galaxy of GRB 030329 Using the Nobeyama Millimeter Array

A sensitive observation of the CO (J=1-0) molecular line emission in the host galaxy of GRB 030329 (z =0.1685) has been performed using the Nobeyama Millimeter Array in order to detect molecular gas and hidden star formation. No sign of CO emission was detected, which invalidates our previous report on the presence of molecular gas. The 3sigma upperlimit on the CO line luminosity (L'_CO) of the host galaxy is 6.9 x 10^8 K km s^-1 pc^2. The lowerlimit of the host galaxy's metallicity is estimated to be 12+log(O/H) ~ 7.9, which yields a CO line luminosity to H_2 conversion factor of alpha_CO = 40 Msun (K km s^-1 pc^2)^-1. Assuming this alpha_CO factor, the 3sigma upperlimit on the molecular gas mass of the host galaxy is 2.8 x 10^10 Msun. Based on the Schmidt law, the 3sigma upperlimit of the total star formation rate (SFR) of the host galaxy is estimated to be 38 Msun yr^-1. These results independently confirm the inferences of previous observations in the optical, submillimeter, and X-ray band, which regard this host galaxy as a compact dwarf galaxy, and not a massive, aggressively star forming galaxy.

preprint2007arXiv

VLBI Observations of Water Masers in Onsala 1: Massive Binary Star Forming Site?

We present proper motions of water masers toward the Onsala 1 star forming region, observed with the Japanese VLBI network at three epochs spanning 290 days. We found that there are two water maser clusters (WMC1 and WMC2) separated from each other by 1".6 (2900 AU at a distance of 1.8 kpc). The proper motion measurement reveals that WMC1 is associated with a bipolar outflow elongated in the east-west direction with an expansion velocity of 69+-11 km/s. WMC1 and WMC2 are associated with two 345 GHz continuum dust emission sources, and located 2" (3600 AU) east from the core of an ultracompact HII region traced by 8.4 GHz radio continuum emission. This indicates that star formation activity of Onsala 1 could move from the west side of ultracompact HII region to the east side of two young stellar objects associated with the water masers. We also find that WMC1 and UC HII region could be gravitationally bound. Their relative velocity along the line of sight is 3 km/s, and total mass is 37 Mo. Onsala 1 seems to harbor a binary star at different evolutionary stage.

preprint2006arXiv

ASTE $^{12}$CO($J=$3--2) Survey of Elliptical Galaxies

We report $^{12}$CO($J=$3--2) observations of 15 nearby elliptical galaxies carried out with the ASTE telescope. Thirteen were selected without regard to the presence of other tracers of cold interstellar matter. CO emission was detected from three of the galaxies, two of which are undetected by IRAS at 100 microns, suggesting that cold ISM may be present in more ellipticals than previously thought. The molecular gas masses range from $2.2 \times 10^6$ to $4.3 \times 10^8$ $M_\odot$. The ratio of the CO(3--2) and (1--0) lines, $R_{31}$, has a lower value for elliptical galaxies than for spiral galaxies except for NGC 855, for which the value is close to the mean for spirals. The molecular gas in NGC 855 has a mean density in the range 300 -- 1000 cm$^{-3}$ adopting a temperature range of 15 -- 100 K.

preprint2006arXiv

ASTE observations of nearby galaxies: A tight correlation between CO(J=3-2) emission and Halpha

Star formation rates (SFRs) obtained via extinction corrected H alpha are compared to dense gas as traced by CO(J=3-2) emission at the centers of nearby galaxies, observed with the ASTE telescope. It is found that, although many of the observed positions are dusty and therefore heavily absorbed at H alpha, the SFR shows a striking correlation with dense gas in the form of the Schmidt law with an index 1.0. The correlation is also compared between gas traced by CO(J=1-0) and application of H alpha extinction correction. We find that dense gas produces a far better correlation with SFR in view of surface density values.

preprint2006arXiv

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: II. The Molecular Gas Disk

We created a three-dimensional distribution map of molecular gas throughout the Milky Way Galaxy using the latest $^{12}$CO($J=$1--0) survey data cube and rotation curve based on the kinematic distance. The radial distribution of the molecular gas shows a central peak and a second peak around 0.5 $R_0$. The thickness of the molecular disk slightly increases from 48 pc to 160 pc with the galactocentric distance within a radius range of 0--11 kpc. We were able to trace the Outer, the Perseus, the Sagittarius-Carina, the Scutum-Crux, and the Norma arms as logarithmic spiral arms with pitch angles of $11\arcdeg - 15\arcdeg$. Considering that the pitch angles of the spiral arms are within this range, the Norma and the Outer arms {seem to} be identified as the same spiral arm. We could also trace a midplane displacement, whose amplitude is nearly constant inside a 10 kpc radius and increases beyond this radius. The ridges of midplane displacement form leading spiral arms.

preprint2006arXiv

Threshold between Spontaneous and Cloud-Collisional Star Formation

Based on simple physical and geometric assumptions, we have calculated the mean surface molecular density of spiral galaxies at the threshold between star formation induced by cloud-cloud collision and spontaneous gravitational collapse. The calculated threshold is approximately $\log Σ_\mathrm{crit} \sim 2.5$, where $ Σ\quad \mathrm{M_{\solar}}\cdot \mathrm{pc}^{-2}$ is the observed surface mass density of an assumed flat gas disk. Above this limit, the rate of molecular cloud collisions dominates over spontaneous molecular cloud collapse. This model may explain the apparent discontinuity in the Schmidt law found recently at $2 \lesssim \log Σ\lesssim 3$.

preprint2005arXiv

Magnetic-Reconnection and Current-Sheet Model for the Radio Arc and Threads in the Galactic Center

We propose a new mechanism to explain the radio Arc and threads in the Galactic Center by current sheets produced by local magnetic shears due to interaction of a moving cloud and the vertical field based on three-dimensional magneto-hydrodynamical simulations. Magnetic reconnetion and acceleration of cosmic-ray electrons in the current sheet will result in high contrast of radio emissivity inside and outside the Arc and threads.

preprint2005arXiv

The Schmidt Law at High Molecular Densities

We have combined Halpha and recent high resolution CO(J=1-0) data to consider the quantitative relation between gas mass and star formation rate, or the so-called Schmidt law in nearby spiral galaxies at regions of high molecular density. The relation between gas quantity and star formation rate has not been previously studied for high density regions, but using high resolution CO data obtained at the NMA(Nobeyama Millimeter Array), we have found that the Schmidt law is valid at densities as high as $10^3 \mathrm{M_\odot} \mathrm{pc}^{-2}$ for the sample spiral galaxies, which is an order of magnitude denser than what has been known to be the maximum density at which the empirical law holds for non-starburst galaxies. Furthermore, we obtain a Schmidt law index of $N=1.33\pm0.09$ and roughly constant star formation efficiency over the entire disk, even within the several hundred parsecs of the nucleus. These results imply that the physics of star formation does not change in the central regions of spiral galaxies. Comparisons with starburst galaxies are also given. We find a possible discontinuity in the Schmidt law between normal and starburst galaxies.

preprint2005arXiv

Virgo High-Resolution CO Survey :V. Circumnuclear Elliptical Ring in NGC 4569

We present high-resolution ($1\farcs8$ -- $4\farcs5$) CO data of the Virgo spiral galaxy NGC 4569, obtained using the Nobeyama Millimeter Array. We found that the molecular gas is highly concentrated in the circumnuclear region with two off-center peaks. A CO image with the highest angular resolution of $2\farcs0 \times 1\farcs8$ shows that six blobs likely form a circumnuclear elliptical ring (CER) with a semi-major axis radius of 720 pc. The CER shows a strongly twisted velocity field, and the position--velocity (PV) diagram shows significant forbidden velocity components. These kinetic features are understood as being non-circular motion due to a bar-potential. We found that the CER coincides with the H$α$ bright central core and that the mass ratio of the molecular gas to the dynamical mass is about 18%. These results support a gaseous inflow scenario induced by a weak bar potential and self-gravity of the gas.

preprint2004arXiv

Offsets between H-alpha and CO arms of a spiral galaxy NGC 4254: A New Method for Determining the Pattern Speed of Spiral Galaxies

We examined offsets between HII regions and molecular clouds belonging to spiral arms of a late type spiral galaxy NGC 4254 (M99). We used a high resolution CO(1-0) image obtained by Nobeyama Millimeter Array (NMA) and an H-alpha image. We derived angular offsets (theta) in the galactic disk, and found that these offsets show a linear dependence on the angular rotation velocity of gas (Omega_G). This linear relation can be expressed by an equation: theta =(Omega_G - Omega_P) * t_{H-alpha}, where Omega_P and t_{H-alpha} are constant. Here, Omega_P corresponds to the pattern speed of spiral arms and t_{H-alpha} is interpreted as the timescale between the peak compression of the molecular gas in spiral arms and the peak of massive star formation. We may thus determine Omega_P and t_{H-alpha} simultaneously by fitting a line to our theta - Omega_G plot, if we assumed they are constant. From our plot, we obtained t_{H-alpha} =4.8 (+/- 1.2) Myr and Omega_P = 26 (+10/-6) km/s/kpc, which are consistent with previous studies. We suggest that this theta - Omega_G plot can be a new tool to determine the pattern speed and the typical timescale needed for star formations.

preprint2004arXiv

Virgo High-Resolution CO Survey IV. Spiral-Driven Gas Dynamics in the Non-Barred Seyfert Galaxy NGC 4501

We report on high-resolution interferometer observations of the 12CO(J=1-0) emission in the central 5 kpc region of the Seyfert 2 galaxy NGC 4501. The observations were made using the Nobeyama Millimeter Array during a long-term CO line survey of Virgo spirals. The major features are: (1) a nuclear concentration with a radius of r~5'' (390 pc), which is resolved into double peaks, and (2) spiral arms which extend out from the nuclear region. The nuclear component has a mass of 1.3x10^8 Mo, which corresponds to ~3.5% of the dynamical mass, and shows a slight non-circular motion. The double peaks are separated by ~4''.7. (370 pc), and located on the root of optical spiral arms in a HST image. The gas arms are associated with the spiral dust lanes, and are linked to the central double peaks. The non-circular motions along the molecular arms indicate the fact that the gas is driven by the density wave, rather than the stochastic processes. We calculated the gas cloud orbits in a stellar spiral potential, and explained the observed CO spiral arms and non-circular motions. We suggest that the central gas condensation arises from spiral-driven gas transfer. We estimated and compared the effect of two possible mechanisms of angular-momentum transfer: galactic shock, and gravitational torques. We discuss that the galactic shock is dominant.

preprint2003arXiv

Faint 6.7um Galaxies and their Contributions to the Stellar Mass Density in the Universe

We discuss the nature of faint 6.7um galaxies detected with the mid-infrared camera ISOCAM on board the Infrared Space Observatory (ISO). The 23 hour integration on the Hawaii Deep Field SSA13 has provided a sample of 65 sources down to 6uJy at 6.7um. For 57 sources, optical or near-infrared counterparts were found with a statistical method. All four Chandra sources, three SCUBA sources, and one VLA/FIRST source in this field were detected at 6.7um with high significance. Using their optical to mid-infrared colors, we divided the 6.7um sample into three categories: low redshift galaxies with past histories of rapid star formation, high redshift ancestors of these, and other star forming galaxies. Rapidly star forming systems at high redshifts dominate the faintest end. Spectroscopically calibrated photometric redshifts were derived from fits to a limited set of template SEDs. They show a high redshift tail in their distribution with faint (<30uJy) galaxies at z>1. The 6.7um galaxies tend to have brighter K magnitudes and redder I-K colors than the blue dwarf population at intermediate redshifts. Stellar masses of the 6.7um galaxies were estimated from their rest-frame near-infrared luminosities. Massive galaxies (M_star~10e11M_sun) were found in the redshift range of z=0.2-3. Epoch dependent stellar mass functions indicate a decline of massive galaxies' comoving space densities with redshift. Even with such a decrease, the contributions of the 6.7um galaxies to the stellar mass density in the universe are found to be comparable to those expected from UV bright galaxies detected in deep optical surveys.

preprint2003arXiv

Galactic Center Shells and a Recurrent Starburst Model

By applying filtering techniques to remove straight filaments in the 20-cm VLA radio image of the Galactic Center Arc region, we have shown that numerous concentric radio shells of radii 5 to 20 pc are surrounding the Pistol and Sickle region, which we call Galactic Center Shells (GCS). Each shell has thermal energy of the order of $10^{49-50}$ erg. Several CO-line shells are associated, whose kinetic energies are of the order of $10^{49-50}$ erg. Summing up the energies of recognized GCSs, the total energy amounts to $\sim 10^{51}$ erg. The GCSs show an excellent correlation with the FIR shells observed at 16--26 microns with the MSX. We propose a model in which GCSs were produced by recurrent and/or intermittent starbursts in the Pistol area during the last million yr. The most recent burst occurred some $10^5$ years ago, producing an inner round-shaped shell (GCS I); earlier ones a million years ago produced outer shells (GCS II and III), which are more deformed by interactions with the surrounding ISM and Sgr A halo. We argue that recurrent starbursts had also occurred in the past, which produced larger scale hyper shell structures as well. A burst some million years ago produced the Galactic Center Lobe, and a much stronger one 15 million years ago produced the North Polar Spur.

preprint2003arXiv

The Virgo High-Resolution CO Survey. III. NGC 4254

We present high-angular-resolution (1".5 - 5") interferometer observations of the \co emission in the central region of the SA(s)c galaxy NGC 4254. The observations were obtained using the Nobeyama Millimeter-wave Array (NMA) during the course of a long-term CO line survey of Virgo spirals. We present the spectra, channel maps, integrated intensity distributions, velocity fields, position--velocity diagrams, and compare the data with various optical images. The rotation velocity is already finite at the nucleus, or at least it rises steeply to 80 km/s within the central 1", indicating the existence of a massive core of 10^8 Msun within 1" (80 pc) radius. The CO intensity maps show that the inner disk has well-developed multiple spiral arms, winding out from a bar-shaped elongated molecular complex. In addition to the bisymmetric spiral arms, an asymmetric tightly wound arm with high molecular gas density is found to wind out from the molecular bar. The molecular spiral arms, particularly the tightly wound arm, well traces optical dark lanes, and are associated with H$α$ arms having many H {\sc ii} regions. The inner asymmetric spiral structures can be explained by ram-pressure distortion of inter-arm low density regions of the inner disk by the intra-cluster gas wind, and is indeed reproduced by a hydrodynamical simulation.

preprint2003arXiv

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: I. The HI Disk

We derived the three-dimensional distribution of HI gas in the Milky Way Galaxy using the latest HI survey data cubes and rotation curves. The distance of the HI gas was determined by the kinematic distance using a rotation curve. We solved the near-far problem in the inner Galaxy by a fitting method which involves introducing a model of vertical HI distribution. In our resultant maps we could trace three prominent arms: the Sagittarius-Carina arm, the Perseus arm, and the Outer arm. These three arms were found to be logarithmic spiral arms. The pitch angles of the Sagittarius-Carina, Perseus, and Outer arms were estimated to be about $11\arcdeg$, $18\arcdeg$, and $7\arcdeg$, respectively. The Sun is located in a region rich in HI gas between the Sagittarius-Carina arm and the Perseus arm. The HI disk shows large and asymmetric warping in the outer disk: the HI disk goes up to about 1.5 kpc above the Galactic plane in the northern hemisphere, and down to about 1 kpc in the southern hemisphere, which means asymmetric warping. The inner HI disk is also found to be tilting. The radius of the HI disk is about 17 kpc and the HI mass within this radius is estimated to be $2.5 \times 10^9$ M$_\odot$, which corresponds to 1.5 % of the dynamical mass predicted from the rotation curve. We also found that the HI outskirt is largely swelling in the fourth quadrant, and hence the Galaxy is significantly lopsided. The scale-height of the HI layer increases with the radius, and is correlated with the HI volume density at the centroid of the HI layer.

preprint2002arXiv

Ram-Ppressure Effects on Dense Molecular Arms in the Central Regions of Spiral Galaxies by Intracluster Medium

We investigated ram-pressure effects by an intracluster wind on an inner disk of spiral galaxies by a hydrodynamical simulation. Even if the wind is mild and not strong enough to strip the gas disk, the ram pressure disturbs orbits of the inter-arm gas significantly. This results in asymmetric dense molecular arms in the inner few kpc region of a galaxy. This mechanism would explain the asymmetric CO gas distributions in the central regions often observed in Virgo spirals.

preprint2002arXiv

Spiral Molecular Front in Galaxies: Quick Transition from Atomic to Molecular Hydrogen in Spiral Arms

We derived a two-dimensional map of the molecular fraction, fmol, (the ratio of molecular-gas density to that of total gas) in the spiral galaxy M51, and examined the behavior of molecular fronts (MF), where MF represents the place where fmol changes drastically from nearly zero to unity and vice versa. We show that the MF phenomenon occurs not only radial, but also in the azimuth direction through the spiral arms, and fmol changes rapidly in the arm-to-inter-arm transition regions. The existence of the azimuthal MF indicates that the atomic gas (HI) is quicklytransformed to molecular gas (H_2) during the passage through spiral arms. We performed a numerical simulation of MF based on an HI-to-H_2 phase transition theory, and reproduced the observations. We estimated a azimuthal scale length of the transition to be less than 200 pc, corresponding to a time scale of about 2 Myr for H_2 gas formation. The azimuthal width of a molecular arm is estimated to be at most 2.5 kpc, where the gas can remain in molecular phase for about 25 Myr.

preprint1999arXiv

Dark Bulge, Exponential Disk, and Massive Halo in the Large Magellanic Cloud

The rotation curve of the Large Magellanic Cloud, which we have derived from high-resolution HI position-velocity diagrams observed by Kim et al (1998), shows a steep central rise and flat rotation with a gradual rise toward the edge. Using the rotation curve, we have calculated the distribution of surface mass density, and show that the LMC has a dark compact bulge, an exponential disk, and a massive halo. The bulge is 1.2 kpc away from the center of the stellar bar, and is not associated with an optical counterpart. This indicates that the "dark bulge" has a large fraction of dark matter, with an anomalously high mass-to-luminosity (M/L) ratio. On the contrary, the stellar bar has a smaller M/L ratio compared to the surrounding regions. (http://www.ioa.s.u-tokyo.ac.jp/~sofue/99lmc.htm)

preprint1998arXiv

Is M82 a Disk-Truncated Bulge by a Close Encounter with M81?

The rotation curve of the small-mass starburst galaxy M82 has a steep nuclear rise, peaking at 200 pc radius, and then, declines in a Keplerian fashion. This rotation curve mimics that for a central bulge of spiral galaxies with high concentration of stellar mass. The declining rotation indicates that its extended disk mass is missing. In order to explain this peculiar rotation characteristics, we propose a hypothesis that M82 is a surviving central bulge of a much larger disk galaxy, whose outer disk was truncated during a close encounter with M81. We simulate a tidal truncation of the disk of a companion galaxy by a tidal penetration through its more massive parent galaxy. The model can well reproduce the observed peculiar feature of M82.

preprint1998arXiv

Nuclear-to-Disk Rotation Curves of Galaxies in the H_alpha and [NII] Emission Lines

We have obtained optical CCD spectroscopy along the major axes of 22 nearby spiral galaxies of Sb and Sc types in order to analyze their rotation curves. By subtracting the stellar continuum emission, we have obtained position velocity (PV) diagrams of the H alpha and [NII] lines. We point out that the H alphaline is often superposed by a broad stellar absorption feature (Balmer wind) in the nuclear regions, and, therefore, the [NII] line is a better tracer of kinematics in the central a few hundred pc regions. By applying the envelope-tracing technique to the H alpha and [NII] PV diagrams, we have derived nucleus-to-disk rotation curves of the observed galaxies. The rotation curves rise steeply within the central a few hundred parsecs, indicating rapidly rotating nuclear disk and mass concentration near the nucleus. [For more rotation curves : http://www.ioa.s.u-tokyo.ac.jp/]

preprint1997arXiv

On the Keplerian Rotation Curves of Galaxies

We propose a criterion for examining whether or not the uncertainty of the outer rotation curves is sufficiently small to distinguish a Keplerian rotation curve from flat ones. We have applied this criterion to both Keplerian and non-Keplerian rotation curves so far obtained, and investigated their relative fraction. We also studied the minimum extent of the dark halos indicated by non-Keplerian rotation curves using the criterion. We have found that one cannot rule out the possibility that a significant fraction of rotation curves become Keplerian within 10-times the disk scale length. If the Keplerian rotation curves so far observed trace the mass truncation, several galaxies may have rather small halos, the extent of which is not larger than twice that of the optical disk.

preprint1996arXiv

CO vs HI Rotations of Galaxies

Using the fact that molecular gas is distributed in the inner disk while HI in the outer disk, we have derived most-completely sampled rotation curves by combining HI and CO position-velocity diagrams. We show that inner rotation curves have a steep increase within the central few hundred-pc region. The inner behavior of rotation can be fitted by the Miyamoto-Nagai potential, only if a central mass component with scale radius of 100 -- 200 pc and a mass of $\sim 3-5\times 10^{9}\Msun$ is assumed in addition to the usual bulge. This implies that the central $\sim 100$-pc region has a stellar mass density as high as several hundred $\Msun{\rm pc}^{-3}$, orders-of-magnitudes higher than that for a single bulge and disk (several $\Msun {\rm pc}^{-3}$). Such a concentric multi-component structure within a bulge (``bulge in bulge'') may have crucial implication for the dynamical evolution of the central regions of galaxies. Since the molecular gas distribution is closely correlated with optical disk, CO line can be used for the line width-luminosity (Tully Fisher) relation in addition to (instead of) the HI line. We show that the CO-line width-luminosity relation can be a powerful routine method to derive distances of cosmologically-remote galaxies and the Hubble parameter for a significantly large volume of the universe. By applying the CO-line width-luminosity relation to CO-line data by large-aperture mm-wave telescope, we have derived distances to galaxies at $cz \sim 20,000$. In order to measure the inclination and magnitudes, we performed optical imaging observations of thus CO-detected galaxies using the CFHT 3.6-m telescope at high resolution. The radio and optical data have been combined to derive the distances to the galaxies. We further derived the Hubble ratio for these galaxies by applying the CO-line

preprint1996arXiv

Mass of the Galaxy Inferred from Outer Rotation Curve

Using an outer rotation curve of the Galaxy, we explore the galactic constants and the mass of the Galaxy. We show that Θ_0 of 200 km/s is more favorable than the IAU standard value of 220 km/s, and also show that if Θ_0 is smaller than 207 km/s the rotation curve beyond 2R_0 is declining in Keplerian fashion. In the case of Θ_0= 200 km/s and R_0= 7.6 kpc, the total mass and the extent of the Galaxy inferred from the rotation curve are 2.0+/-0.3x 10^{11} M_\odot and 15 kpc, respectively. These results may significantly change the previous view of the Galaxy, that its outer region is dominated by a massive dark halo extending out to several tens of kpc.

preprint1996arXiv

Nuclear-to-Outer Rotation Curves of Galaxies in the CO and HI Lines

We have derived rotation curves of nearby galaxies, which are almost completely sampled from the inner to outer regions. We used high-resolution position-velocity diagrams in the CO line along the major axes for the inner regions, and HI-line data for the outer regions. We combined the CO and HI data to obtain total rotation curves from the nuclear to outer regions. The rotation curves show a steep nuclear rise within a radius smaller than the resolutions of CO-line observations, indicating a compact massive concentration near the nucleus. We show that the nuclear steep rise is general and universal, and a rigid body-like gentle rise is exceptional. Keywords: Galaxies: kinematics - Galaxies: rotation - Galaxies: structure - Galaxies: general - ISM: CO emission - ISM: HI gas

preprint1996arXiv

Radio Continuum and Molecular Gas in the Galactic Center

Nonthermal radio emission in the galactic center reveals a number of vertical structures across the galactic plane, which are attributed to poloidal magnetic field and/or energetic outflow. Thermal radio emission comprises star forming regions distributed in a thin, dense thermal gas disk. The thermal region is associated with dense molecular gas disk, in which the majority of gas is concentrated in a rotating molecular ring. Outflow structures like the radio lobe is associated with rotating molecular gas at high speed, consistent with a twisted magnetic cylinder driven by accretion of a rotating gas disk.

preprint1995arXiv

Effects of Starbursts on the Structure of Young Galaxies

We have investigated the effects of primeval starbursts on the galactic structure, showing that the dynamical structure of young galaxies is changed by relaxation of a stellar system formed by a starburst-driven outflowing gas (superwind). When a superwind interacts with halo gas, an expanding dense supershell of shocked gas is formed. Using a similarity solution for the expanding shell and a condition for the gravitational instability of a gaseous shell, we show that an expanding gaseous shell with a mass of several $10^{10} M_\odot$ and a radius of several kpc becomes gravitationally unstable to form stars. A stellar shell is thus formed, relaxes and will evolve into a fat stellar system. \par In order to investigate the fate of the stellar shell and its dynamical influence on the host galaxy, we performed three dimensional $N$-body simulations of a stellar shell plus disk system. The evolution of the stellar shell was computed under two types of environment: an external flattened potential and a live stellar disk. We find that the relaxation process and resultant structure are significantly affected by the external disk potential, and also depend on the velocity dispersion of the initial shell stars and the mass of the shell. The final galactic structure ranges from a thick disk and a compact bulge with a high density core to a thin disk with a diffuse bulge. Isodensity maps of the computed disk plus bulge system resemble isophotes of S0 galaxies. We suggest that the effects of primeval starbursts would explain some observational properties of S0 galaxies.