Researcher profile

Keiichi Wada

Keiichi Wada contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

29 published item(s)

preprint2026arXiv

The Origin of Spin-Alignment of Dark Matter Subhalos

Subhalo spin is essential for modeling galaxy formation and controlling systematic uncertainties in intrinsic alignment (IA) studies. However, the physical mechanisms governing subhalo spin acquisition within the tidal environments of host halos remain poorly understood. In this work, we investigate the alignment between subhalo and host halo spins using the high-resolution cosmological $N$-body simulation, Shin-Uchuu. We find that the spin alignment between subhalos and host halos becomes increasingly pronounced toward the central regions. Our analysis reveals that subhalos typically acquire spin in the same direction as their orbital angular momentum. Since the orbital angular momentum of most subhalos is aligned with the host halo spin, an overall alignment between subhalo and host spins emerges. When classified by orbital orientation, however, subhalo spins in the inner regions are found to be oriented perpendicularly or anti-parallel to the host spin for polar and retrograde orbits, respectively. These results provide strong evidence that subhalo spins are acquired through torques exerted by the tidal field of the host halo. Furthermore, we demonstrate that the mass ratio and the radial distance from the host center are the primary parameters governing subhalo spin alignment, while the dependence on the accretion redshift is found to be negligible.

preprint2022arXiv

Origin and evolution of dust-obscured galaxies in galaxy mergers

Dust Obscured Galaxies (DOGs), which are observationally characterized as faint in the optical and bright in the infrared, are the final stage of galaxy mergers and are essential objects in the evolution of galaxies and active galactic nuclei (AGNs). However, the relationship between torus-scale gas dynamics around AGNs and DOGs lifetime remain unclear. We obtained evolution of the spectral energy distributions (SEDs) of a galaxy merger system with AGN feedback, from post-processed pseudo-observations based on an N-body/Smoothed Particle Hydrodynamics (SPH) simulation. We focused on a late stage merger of two identical galaxies with a supermassive black hole (SMBH) of 10$^8$ M$_\odot$. We found that the infrared luminosity of the system reaches ultra- and hyper-luminous infrared galaxy classes (10$^{12}$ and 10$^{13}$ L$_\odot$, respectively). The DOGs phase corresponds to a state in which the AGNs are buried in dense gas and dust, with the infrared luminosity exceeding 3.3 $\times$ 10$^{12}$ L$_\odot$. We also identified the sub-categories of DOGs, namely bump and power-law DOGs from the SEDs and their evolution. The bump DOGs tend to evolve to power-law DOGs on several Myrs. We found that contribution from the hot dust around the nucleus in the infrared radiation is essential for identifying the system as a power-law DOG; the gas and dust distribute non-spherically around the nucleus, therefore, the observed properties of DOGs depend on the viewing angle. In our model, the lifetime of merger-driven DOGs is less than 4 Myrs, suggesting that the observed DOGs phase is a brief aspect of galaxy mergers.

preprint2022arXiv

Probing dynamics and thermal properties inside molecular tori with CO rovibrational absorption lines

A recent hydrodynamic model, "radiation-driven fountain model" (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers confirmed that this picture is consistent with multi-wavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovibrational absorption lines ($ΔJ=\pm1$, $v=0-1$, $λ\sim 4.7$ $\mathrm{μm}$) could probe physical properties of the inside tori. However, the origin of the CO absorption lines has been under debate. In this paper, we investigate the origin of the absorption lines and conditions for detecting them by performing line radiative transfer calculations based on the radiation-driven fountain model. We find that CO rovibrational absorption lines are detected at inclination angles $θ_\mathrm{obs} = 50-80$ $^{\circ}$. At the inclination angle $θ_\mathrm{obs} = 77$ $^{\circ}$, we observe multi-velocity components: inflow ($v_\mathrm{LOS}=30$ $\mathrm{kms^{-1}}$), systemic ($v_\mathrm{LOS}=0 \, \mathrm{kms^{-1}}$), and outflows ($v_\mathrm{LOS}=-75,\, -95,$ and $-105$ $\mathrm{kms^{-1}}$). The inflow and outflow components ($v_\mathrm{LOS}= 30$ and $-95$ $\mathrm{kms^{-1}}$) are collisionally excited at the excitation temperature of $186$ and $380$ K up to $J=12$ and $4$, respectively. The inflow and outflow components originate from the accreting gas on the equatorial plane at $1.5$ pc from the AGN center and the outflowing gas driven by AGN radiation pressure at $1.0$ pc, respectively. These results suggest that CO rovibrational absorption lines can provide us with the velocities and kinetic temperatures of the inflow and outflow in the inner a-few-pc regions of AGN tori, and the observations can probe the gas circulation inside the tori.

preprint2021arXiv

Spontaneous Formation of Outflows Powered by Rotating Magnetized Accretion Flows in a Galactic Center

We investigate how magnetically driven outflows are powered by a rotating, weakly magnetized accretion flow onto a supermassive black hole using axisymmetric magnetohydrodynamic simulations. Our proposed model focuses on the accretion dynamics on an intermediate scale between the Schwarzschild radius and the galactic scale, which is $\sim$1-100 pc. We demonstrate that a rotating disk formed on a parsec-scale acquires poloidal magnetic fields via accretion and this produces an asymmetric bipolar outflow at some point. The formation of the outflow was found to follow the growth of strongly magnetized regions around disk surfaces (magnetic bubbles). The bipolar outflow grew continuously inside the expanding bubbles. We theoretically derived the growth condition of magnetic bubbles for our model that corresponds to a necessary condition for outflow growth. We found that the north-south asymmetric structure of the bipolar outflow originates from the complex motions excited by accreting flows around the outer edge of the disk. The bipolar outflow comprises multiple mini-outflows and downflows (failed outflows). The mini-outflows emanate from the magnetic concentrations (magnetic patches). The magnetic patches exhibit inward drifting motions, thereby making the outflows unsteady. We demonstrate that the inward drift can be modeled using a simple magnetic patch model that considers magnetic angular momentum extraction. This study could be helpful for understanding how asymmetric and non-steady outflows with complex substructures are produced around supermassive black holes without the help of strong radiation from accretion disks or entrainment by radio jets such as molecular outflows in radio-quiet active galactic nuclei, NGC 1377.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): A multiwavelength view of WISE mid-infrared galaxies/active galactic nuclei

We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset. We have also performed (i) an X-ray spectral analysis, (ii) SED fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru HSC images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7,088 WISE W4 objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and/or probably obscured objects. We find that (i) 82% of the eFEDS-W4 sources are classified as X-ray AGN with $\log\,L_{\rm X} >$ 42 erg s$^{-1}$; (ii) 67% and 24% of the objects have $\log\,(L_{\rm IR}/L_{\odot}) > 12$ and 13, respectively; (iii) the relationship between $L_{\rm X}$ and the 6 $μ$m luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and $N_{\rm H}$ for the eFEDS-W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5% of the eFEDS-W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE W4 sample of AGNs at $z < 4$, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey.

preprint2021arXiv

Warm Absorbers in the Radiation-driven Fountain Model of Low-mass Active Galactic Nuclei

To investigate the origins of the warm absorbers in active galactic nuclei (AGNs), we study the ionization-state structure of the radiation-driven fountain model in a low-mass AGN (Wada et al. 2016) and calculate the predicted X-ray spectra, utilizing the spectral synthesis code Cloudy (Ferland et al. 2017). The spectra show many absorption and emission line features originated in the outflowing ionized gas. The O VIII 0.654 keV lines are produced mainly in the polar region much closer to the SMBH than the optical narrow line regions. The absorption measure distribution of the ionization parameter ($ξ$) at a low inclination spreads over 4 orders of magnitude in $ξ$, indicating multi-phase ionization structure of the outflow, as actually observed in many type-1 AGNs. We compare our simulated spectra with the high energy-resolution spectrum of the narrow line Seyfert 1 galaxy, NGC 4051. The model reproduces slowly outflowing (a few hundreds km s$^{-1}$) warm absorbers. However, the faster components with a few thousands km s$^{-1}$ observed in NGC 4051 are not reproduced. The simulation also underproduces the intensity and width of the O VIII 0.654 keV line. These results suggest that the ionized gas launched from sub-parsec or smaller regions inside the torus, which are not included in the current fountain model, must be important ingredients of the warm absorbers with a few thousands km s$^{-1}$. The model also consistently explains the Chandra/HETG spectrum of the Seyfert 2 galaxy, the Circinus galaxy.

preprint2019arXiv

Mass accretion toward black holes in the final phase of galaxy mergers

We studied the final phases of galactic mergers, focusing on interactions between supermassive black holes (SMBHs) and the interstellar medium in a central sub-kpc region, using an N-body/hydrodynamics code. This numerical experiment aims to understand the fate of the gas supplied by mergers of two or more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the mass accretion rate to one SMBH exceeds the Eddington accretion rate when the distance between two black holes (BHs) rapidly decreases. However, this rapid accretion phase does not last for more than $10^7$ yrs, and it drops to $\sim$ 10% of the Eddington rate in the quasi-steady accretion phase. The rapid accretion is caused by the angular momentum transfer from the gas to the stellar component, and the moderate accretion in the quasi-steady phase is caused {by turbulent viscosity and gravitational torque in the disk. The second merger event enhances the mass accretion to the BHs; however, this phase takes place on a similar timescale to the first merger event. We also found that the AGN feedback and the mass accretion to BHs can coexist in the central region of merged galaxies, if the amount of feedback energy is given as $(2 \times 10^{-4} - 2 \times 10^{-3} )\dot{M} c^2$, where $\dot{M}$ is the accretion rate to $r= 1$ pc. The accretion rate is suppressed by $\sim$ 1/50 in the quasi-steady accretion phase for $0.02 \dot{M} c^2$. The fraction of the gas that finally falls to each BH is approximately 5-7% of the supplied total gas mass ($10^8 M_\odot$), and 15-20% of the gas forms a circumnuclear gas inside 100 pc. This remnant gas heavily obscures the luminous phase of the active galactic nuclei (AGN) during merger events, and the moderate AGN feedback does not alter this property.

preprint2019arXiv

Obscuring fraction of active galactic nuclei implied by supernova and radiative feedbacks

We study the obscuring structure of circumnuclear disks (CNDs) by considering supernova (SN) feedbacks from nuclear starburst and the effect of anisotropic radiative pressure from AGNs. We suppose that the mass accretion onto a central supermassive black hole (SMBH) is triggered by SN-driven turbulence within CNDs, and we explore how the structures of CNDs depend on the BH mass ($M_{\rm BH}$) and AGN luminosity ($L_{\rm AGN}$). We find that the obscuring fraction ($f_{\rm obs}$) peaks at $\sim10\%$ of the Eddington luminosity ($L_{\rm Edd}$), and its maximal value is $f_\mathrm{obs} \sim 0.6$ for less massive SMBHs (e.g., $M_{\rm BH} < 10^{8}M_{\odot}$). This is because the scale height of CNDs is determined by the SN-driven accretion for a smaller $L_{\rm AGN}$, while the dusty molecular gas in CNDs is blown away by the radiation pressure from AGNs beyond the critical luminosity. On the other hand, for massive SMBHs (e.g., $M_{\rm BH} > 10^{8}M_{\odot}$), $f_{\rm obs}$ is always smaller than $0.2$, and it is almost independent of $L_{\rm AGN}$ because the scale height of CNDs is mainly controlled by the maximal star-formation efficiency ($C_{\rm *, max}$) in CNDs. By comparison with the obscuring fractions suggested from the mid-infrared observations of nearby AGNs, the SN plus radiative feedback model with $C_{\rm *, max} = 10^{-7}\, {\rm yr}^{-1}$ well reproduces the observations for $M_\mathrm{BH} = 10^8 M_\odot$. We also find that the intense starburst or the existence of dust-free absorbers inside CNDs are necessary, to explain X-ray observations.

preprint2019arXiv

Planet Formation around Super Massive Black Holes in the Active Galactic Nuclei

As a natural consequence of the elementary processes of dust growth, we discovered that a new class of planets can be formed around supermassive black holes (SMBHs). We investigated a growth path from sub-micron sized icy dust monomers to Earth-sized bodies outside the &#34;snow line&#34;, located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the &#34;radial drift barrier&#34; does not prevent the formation of planetesimals. In the early phase of the evolution, low collision velocity between dust particles promotes sticking; therefore, the internal density of the dust aggregates decreases with growth. When the porous aggregate&#39;s size reaches 0.1--1 cm, the collisional compression becomes effective, and the decrease in internal density stops. Once 10--100 m sized aggregates are formed, they are decoupled from gas turbulence, and the aggregate layer becomes gravitationally unstable, leading to the formation of planets by the fragmentation of the layer, with ten times the mass of the earth. The growth time scale depends on the turbulent strength of the circumnuclear disk and the black hole mass $M_{BH}$, and it is comparable to the AGN&#39;s lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs.

preprint2016arXiv

Multi-phase Nature of a Radiation-Driven Fountain with Nuclear Starburst in a Low-mass Active Galactic Nucleus

The structures and dynamics of molecular, atomic, and ionized gases are studied around a low-luminosity active galactic nucleus (AGN) with a small ($2\times 10^6 M_\odot$) black hole using 3D radiation hydrodynamic simulations. We studied, for the first time, the non-equilibrium chemistry for the X-ray dominated region in the &#34;radiation-driven fountain&#34; (Wada 2012) with supernova feedback. A double hollow cone structure is naturally formed without postulating a thick &#34;torus&#34; around a central source. The cone is occupied with an inhomogeneous, diffuse ionized gas and surrounded by a geometrically thick ($h/r \gtrsim 1$) atomic gas. Dense molecular gases are distributed near the equatorial plane, and energy feedback from supernovae enhances their scale height. Molecular hydrogen exists in a hot phase ( > 1000 K) as well as in a cold ( < 100 K), dense ( >$10^3$ cm$^{-3}$) phase. The velocity dispersion of H$_2$ in the vertical direction is comparable to the rotational velocity, which is consistent with near infrared observations of nearby Seyfert galaxies. Using 3D radiation transfer calculations for the dust emission, we find polar emission in the mid-infrared band (12$μm$), which is associated with bipolar outflows, as suggested in recent interferometric observations of nearby AGNs. If the viewing angle for the nucleus is larger than 75 deg, the spectral energy distribution (~ 2 -- 60 $μm$) of this model is consistent with that of the Circinus galaxy. The multi-phase interstellar medium observed in optical/infrared and X-ray observations is also discussed.

preprint2015arXiv

Radiation Hydrodynamic Simulations of Line-Driven Disk Winds for Ultra Fast Outflows

Using two-dimensional radiation hydrodynamic simulations, we investigate origin of the ultra fast outflows (UFOs) that are often observed in luminous active galactic nuclei (AGNs). We found that the radiation force due to the spectral lines generates strong winds (line-driven disk winds) that are launched from the inner region of accretion disks (~30 Schwarzschild radii). A wide range of black hole masses ($M_{\rm BH}$) and Eddington ratios ($\varepsilon$) was investigated to study conditions for causing the line-driven winds. For $M_{\rm BH} = 10^6-10^9 M_\odot$ and $\varepsilon = 0.1-0.7$, funnel-shaped disk winds appear, in which dense matter is accelerated outward with an opening angle of 70-80 deg and with 10% of the light speed. If we observe the wind along its direction, the velocity, the column density, and the ionization state are consistent with those of the observed UFOs. As long as the obscuration by the torus does not affect the observations of X-ray bands, the UFOs could be statistically observed in about 13-28% of the luminous AGNs, which is not inconsistent with the observed ratio (~40%). We also found that the results are insensitive to the X-ray luminosity and the density of the disk surface. Thus, we can conclude that the UFOs could exist in any luminous AGNs, such as narrow-line Seyfert 1s (NLS1s) and quasars with $\varepsilon > 0.1$, in which fast line-driven winds are associated.

preprint2012arXiv

Analysis of H2O Masers in Sharpless 269 using VERA Archival data --- Effect of maser structures on astrometric accuracy

Astrometry using H2O maser sources in star forming regions is expected to be a powerful tool to study the structures and dynamics of our Galaxy. Honma et al. (2007) (hereafter H2007) claimed that the annual parallax of S269 is determined within an error of 0.008 milliarcsec (mas), concluding that S269 is located at 5.3 kpc +- 0.2 kpc from the sun, and R= 13.1 kpc. They claimed that the rotational velocity of S269 is equal to that of the sun within a 3% error. This small error, however, is hardly understood when taking into account the results of other observations and theoretical studies of galactic dynamics. We here reanalyzed the VERA archival data using the self-calibration method (hybrid mapping), and found that clusters of maser features of S269 are distributed in much wider area than that investigated in H2007. We confirmed that, if we make a narrow region image without considering the presence of multiple maser spots, and only the phase calibration is applied, we can reproduce the same maser structures in H2007. The distribution extent of maser spots in the feature differs 0.2 mas from east to west between our results and H2007. Moreover, we found that change of relative positions of maser spots in the cluster reaches 0.1 mas or larger between observational epochs. This suggests that if one simply assumes the time-dependent, widely distributed maser sources as a stable single point source, it could cause errors of up to 0.1 mas in the annual parallax of S269. Taking into account the internal motions of maser spot clusters, the proper motion of S269 cannot be determined precisely. We estimated that the peculiar motion of S269 with respect to a Galactic circular rotation is ~20 km/s. These results imply that the observed kinematics of maser emissions in S269 cannot give a strong constraint on dynamics of the outer part of the Galaxy, in contrast to the claim by H2007.

preprint2012arXiv

Dynamics of Non-Steady Spiral Arms in Disk Galaxies

In order to understand the physical mechanisms underlying non-steady stellar spiral arms in disk galaxies, we analyzed the growing and damping phases of their spiral arms using three-dimensional $N$-body simulations. We confirmed that the spiral arms are formed due to a swing amplification mechanism that reinforces density enhancement as a seeded wake. In the damping phase, the Coriolis force exerted on a portion of the arm surpasses the gravitational force that acts to shrink the portion. Consequently, the stars in the portion escape from the arm, and subsequently they form a new arm at a different location. The time-dependent nature of the spiral arms are originated in the continual repetition of this non-linear phenomenon. Since a spiral arm does not rigidly rotate, but follows the galactic differential rotation, the stars in the arm rotate at almost the same rate as the arm. In other words, every single position in the arm can be regarded as the co-rotation point. Due to interaction with their host arms, the energy and angular momentum of the stars change, thereby causing the radial migration of the stars. During this process, the kinetic energy of random motion (random energy) of the stars does not significantly increase, and the disk remains dynamically cold. Owing to this low degree of disk heating, the short-lived spiral arms can recurrently develop over many rotational periods. The resultant structure of the spiral arms in the $N$-body simulations is consistent with some observational nature of spiral galaxies. We conclude that the formation and structure of spiral arms in isolated disk galaxies can be reasonably understood by non-linear interactions between a spiral arm and its constituent stars.

preprint2012arXiv

Radiation-driven Fountain and Origin of Torus around Active Galactic Nuclei

We propose a plausible mechanism to explain the formation of the so-called &#34;obscuring tori&#34; around active galactic nuclei (AGNs) based on three-dimensional hydrodynamic simulations including radiative feedback from the central source. The X-ray heating and radiation pressure on the gas are explicitly calculated using a ray-tracing method. This radiation feedback drives a &#34;fountain&#34;, that is, a vertical circulation of gas in the central a few to tens parsecs. Interaction between the non-steady outflows and inflows causes the formation of a geometrically thick torus with internal turbulent motion. As a result, the AGN is obscured for a wide range of solid angles. In a quasi-steady state, the opening angles for the column density toward a black hole < 10^23 cm^-2 are approximately +-30 deg and +-50 deg for AGNs with 10% and 1% Eddington luminosity, respectively. Mass inflows through the torus coexist with the outflow and internal turbulent motion, and the average mass accretion rate to the central parsec region is 2x10^-4 ~ 10^-3, M_sun/yr this is about ten times smaller than accretion rate required to maintain the AGN luminosity. This implies that relatively luminous AGN activity is intrinsically intermittent or that there are other mechanisms, such as stellar energy feedback, that enhance the mass accretion to the center.

preprint2011arXiv

Interplay between Stellar Spirals and the ISM in Galactic Disks

We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic simulations in a `live&#39; stellar disk. We focus especially on spiral structures excited in a isolated galactic disk without a stellar bar. Using high-resolution, 3-dimensional N-body/SPH simulations, we found that the spiral features of the gas in galactic disks are formed by essentially different mechanisms from the galactic shock in stellar density waves. The stellar spiral arms and the interstellar matter on average corotate in a galactic potential at any radii. Unlike the stream motions in the galactic shock, the interstellar matter flows into the local potential minima with irregular motions. The flows converge to form dense gas clouds/filaments near the bottom of the stellar spirals, whose global structures resemble dust-lanes seen in late-type spiral galaxies. The stellar arms are non-steady; they are wound and stretched by the galactic shear, and thus local densities of the arm change on a time scale of ~ 100 Myrs, due to bifurcating or merging with other arms. This makes the gas spirals associated with the stellar arms non-steady. The association of dense gas clouds are eventually dissolved into inter-arm regions with non-cirucular motions. Star clusters are formed from the cold, dense gases, whose ages are less than ~30 Myrs, and they are roughly associated with the background stellar arms without a clear spatial offset between gas spiral arms and distribution of young stars.

preprint2011arXiv

Nuclear Star-Forming Ring of the Milky Way: Simulations

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way&#39;s nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past.

preprint2011arXiv

Origin of multiple nuclei in ultraluminous infrared galaxies

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

preprint2011arXiv

Shock-induced star cluster formation in colliding galaxies

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

preprint2011arXiv

The Structure and Dynamics of an AGN Torus: CO Line Predictions for ALMA from 3D Hydrodynamical Simulations with X-ray Driven Chemistry

Many efforts have been made to model the mass distribution and dynamical evolution of the circumnuclear gas in active galactic nuclei (AGNs). However, chemical evolution is not included in detail in three-dimensional (3-D) hydrodynamic simulations. The X-ray radiation from the AGN can drive the gas chemistry and affect the thermodynamics, as well as the excitation of the interstellar medium (ISM). Therefore, we estimate the effects (on chemical abundances and excitation) of X-ray irradiation by the AGN, for atomic and molecular gas in a 3-D hydrodynamic model of an AGN torus. We obtain the abundances of various species from an X-ray chemical model. A 3-D radiative transfer code estimates the level populations, which result in line intensity maps. Predictions for the CO J=1-0 to J=9-8 lines indicate that mid-J CO lines are excellent probes of density and dynamics in the central (<60 pc) region of the AGN, in contrast to the low-J CO lines. Analysis of the X_CO/αconversion factors shows that only the higher-J CO lines can be used for gas mass determination in AGN tori. The [C II] 158 um emission traces mostly the hot (T_k>1000m K) central (<60 pc) region of the AGN torus. The [C II] 158 um line will be useful for ALMA observations of high redshift (z>1) AGNs. The spatial scales (>0.25 pc) probed with our simulations match the size of the structures that ALMA will resolve in nearby (<45 Mpc at 0.01&#34;) galaxies.

preprint2010arXiv

On the Interpretation of the l-v Features in the Milky Way Galaxy

We model the gas dynamics of barred galaxies using a three-dimensional, high-resolution, $N$-body+hydrodynamical simulation and apply it to the Milky Way in an attempt to reproduce both the large-scale structure and the clumpy morphology observed in Galactic H\emissiontype{I} and CO $l-v$ diagrams. Owing to including the multi-phase interstellar medium, self-gravity, star-formation and supernovae feedback, the clumpy morphology, as well as the large-scale features, in observed $l-v$ diagrams are naturally reproduced. We identify in our $l-v$ diagrams with a number of not only large-scale peculiar features such as the &#39;3-kpc arm&#39;, &#39;135-km s$^{-1}$ arm&#39; and &#39;Connecting arm&#39; but also clumpy features such as `Bania clumps&#39;, and then link these features in a face-on view of our model. We give suggestions on the real structure of the Milky Way and on the fate of gas clumps in the central region.

preprint2009arXiv

Formation of High-redshift (z>6) Quasars Driven by Nuclear Starbursts

Based on the physical model of a supermassive black hole (SMBH) growth via gas accretion in a circumnuclear disk (CND) proposed by Kawakatu & Wada (2008), we describe the formation of high-$z$ ($z > 6$) quasars (QSOs) whose BH masses are M_{BH}> 10^{9} M_{\odot}. We derive the necessary conditions to form QSOs at z > 6 by only gas accretion: (i) A large mass supply with M_{sup} > 10^{10}M_{\odot} from host galaxies to CNDs, because the final BH mass is only 1-10% of the total supplied mass from QSO hosts. (ii) High star formation efficiency for a rapid BH growth. We also find that if the BH growth is limited by the Eddington accretion, the final BH mass is greatly suppressed. Thus, the super-Eddington growth is required for the QSO formation. The evolution of the QSO luminosity depends on the redshift z_{i} at which accretion onto a seed BH is initiated. In other words, the brighter QSOs at z >6 favor the late growth of SMBHs (i.e., z_{i}=10) rather than early growth (i.e., z_{i}=30). Moreover, we predict the observable properties and the evolution of QSOs at z >6. In a QSO phase, there should exist a stellar rich massive CND, whose gas mass is about 10% of the dynamical mass inside 0.1-1 kpc}. On the other hand, in a phase where the BH grows (i.e., a proto-QSO phase), the proto-QSO has a gas rich massive CNDs whose gas mass is comparable to the dynamical mass (abridged).

preprint2009arXiv

Molecular Gas Disk Structures around AGNs

We present new high resolution numerical simulations of the ISM in a central R ~32 parsecs region around a supermassive black hole at a galactic center. Three-dimensional hydrodynamic modeling of the ISM (Wada & Norman 2002) with the nuclear starburst now includes tracking of the formation of molecular hydrogen out of the neutral hydrogen phase. In a quasi equilibrium state, mass fraction of H2 is about 0.4 (total H2 mass is ~1.5 10^6 Msun) of the total gas mass for the uniform far UV (FUV) with G_0 = 10. The gas forms an inhomogeneous disk, whose scale-height becomes larger in the outer region. H2 forms a thin nuclear disk in the inner ~ 5 pc, which is surrounded by molecular clouds swelled up toward h < 10 pc. The velocity field of the disk is highly turbulent in the torus region, whose velocity dispersion is ~ 20 km/s on average. Average supernova rate (SNR) of ~ 5 10^-5/yr is large enough to energize these structures. Gas column densities toward the nucleus larger than 10^22 cm^-2 are observed if the viewing angle is smaller than θ_v ~ 50 deg from the edge-on. However, the column densities are distributed over almost two orders of magnitude around the average for any given viewing angle due to the clumpy nature of the torus. For a stronger FUV (G_0 =100), the total H2 mass in an equillibrium is only slightly smaller (~ 0.35). Finally the morphology and kinematics of the circumnuclear molecular gas disks emerging from our models is similar to that revealed by recent near infrared observations using VLTI/Keck.

preprint2009arXiv

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

preprint2008arXiv

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

preprint2007arXiv

Density structure of the interstellar medium and the star formation rate in galactic disks

The probability distribution functions (PDF) of density of the ISM in galactic disks and global star formation rate are discussed. Three-dimensional hydrodynamic simulations show that the PDFs in globally stable, inhomogeneous ISM in galactic disks are well fitted by a single log-normal function over a wide density range. The dispersion of the log-normal PDF (LN-PDF) is larger for more gas-rich systems, and whereas the characteristic density of LN-PDF, for which the volume fraction becomes the maximum, does not significantly depend on the initial conditions. %At the characteristic density, thermal pressure and kinetic %pressure due to turbulent motion are comparable, and the material is %statistically stagnated in a turbulent flow. Supposing the galactic ISM is characterized by the LN-PDF, we give a global star formation rate (SFR) as a function of average gas density, a critical local density for star formation, and star formation efficiency. We find that the observed SFR is well-fitted by the theoretical SFR in a wide range of the global gas density ($10 - 10^4 M_\odot$ pc$^{-2}$). Star formation efficiency (SFE) for high density gas ($n > 10^3$ cm$^{-3}$) is SFE $= 0.001 - 0.01$ for normal spiral galaxies, and SFE $= 0.01 - 0.1$ for starburst galaxies. The LN-PDF and SFR proposed here could be applicable for modeling star formation on a kpc-scale in galaxies or numerical simulations of galaxy formation, in which the numerical resolution is not fine enough to describe the local star formation.

preprint2007arXiv

Origin of Warm High-Velocity Dense Gas in ULIRGs

Possible origins of the molecular absorption discovered in some ULIRGs are investigated, based on a 3-D hydrodynamic model of star-forming interstellar gas in a galactic central region. The blue-shifted, warm ($\sim 200-300$ K), dense ($>10^6$ cm$^{-3}$) molecular gas suggested by CO absorption in IRAS 08752+3915 could be caused by the innermost region of the inhomogeneous inter-stellar medium (ISM) around a supermassive black hole. The infrequent observations of the dense gas with absorption in ULIRGs and Seyfert 2 galaxies could simply suggest that the high-density regions occupy only a very small volume fraction of the obscuring material. This is naturally expected if the inhomogeneous structure of the ISM is caused by non-linear development of instabilities. The model predicts a turbulent velocity field in the obscuring material, therefore blue- and red-shifted gases should be observable with nearly the same probability for the large enough statistical samples.

preprint1998arXiv

The Effect of Self-gravity of Gas on Gas Fueling in Barred Galaxies with a Supermassive Black Hole

In our previous paper, we have shown that a gas disk in the nuclear region of a barred galaxy which contains a central supermassive black hole (SMBH) rapidly evolves into a nuclear gas ring by the effect of an additional inner Lindblad resonance caused by the SMBH. In this paper, we investigate the fate of the gas ring, involving self-gravity of gas, using two-dimensional hydrodynamical simulations. We find that the gas ring becomes gravitationally unstable for a gas surface density of gas above a critical value, and fragments into several gas clumps. Some denser clumps increase their mass via the accretion of the surrounding gas and collisions with other clumps, and finally a very massive gas clump (10^7 M_sun) is formed. Due to the torque from the most massive clump, a part of the gas in the ring loses its angular momentum and falls into the galactic center. As a result, a nuclear gas disk (50 pc) is formed around the SMBH. The accretion rate for $R<50$ pc attains about 1 M_sun/yr for 3.5*10^7 yr. At the final phase of the bar-driven fueling, self-gravity is crucial for the angular momentum transfer of the gas. This is a new mechanism for gas fueling to the vicinity of the SMBH.

preprint1995arXiv

Bar-driven Fueling to a Galactic Central Reagin in a Massive Gas Disk

We have found an effective fueling process to a central region of galaxies with weak bar-like distortion by two dimensional hydrodynamical simulations. Gravitational instability of an elongated gas ring at the inner Lindblad resonance (ILR), which has been reported as an effective fueling mechanism, are not needed for this fueling process. A massive gaseous disk in a central region of galaxies sensitively responds to the weakly distorted potential, and a large amount of gas can be fed into within $1/20$ of a core radius of the potential in several $10^7$ yr. The ILRs, the dissipative nature of the gas, and self-gravity of the gas are essential for triggering this effective fueling. The accumulation process has not been ever known: the gas accumulates to form a dense `linear&#39; structure inclined at about 45 degree with respect to the bar potential in a {\it leading sense}. We also found that a counter rotating gaseous core can be formed as a result of the fueling. The sense of the rotation of the core depends on a fraction of the gas mass to the background mass. Physical mechanism of the fueling process can be understood using a linear theory of gaseous orbits in a weak barred potential.

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.