Search

We investigated the influence of the Eddington ratio on sub-parsec-scale outflows in active galactic nuclei (AGNs) with supermassive black holes (SMBHs) masses of $10^7$ M$_{\odot}$ using two-dimensional radiation hydrodynamics simulations. When the range of Eddington ratio, $\gamma_{\rm Edd} > 10^{-3}$, the radiation force exceeds the gas pressure, leading to stronger outflows and larger dust sublimation radius. Although the sub-parsec-scale outflows is a time-dependence phenomena, our simulations demonstrated that the radial distributions can be well explained by the steady solutions of the spherically symmetric stellar winds. The dynamic structure of sub-parsec-scale outflows is influenced by the dust sublimation radius and the critical radii determined by the dynamical equilibrium condition. Although significantly affecting the outflow velocity, the Eddington ratio exerts minimal effects on temperature and number density distribution. Furthermore, our analytical solutions highlight the importance of the dust sublimation scale as a crucial determinant of terminal velocity and column density in dusty outflows. Through comparisons of our numerical model with the obscuring fraction observed in nearby AGNs, we revealed insights into the Eddington ratio dependence and the tendency towards the large obscuring fraction of the dusty and dust-free gases. The analytical solutions are expected to facilitate an understanding of the dynamical structure and radiation structures along the line of sight and their viewing angles from observations of ionized outflows., Comment: 22 pages, 10 figures, accepted for publication in ApJ

We included the intrinsic line profiles of the strongest fluorescent lines in the X-ray radiative transfer code SKIRT to model the cold-gas structure and kinematics based on high-resolution line observations from XRISM/Resolve and Athena/X-IFU. The intrinsic line profiles of the Ka and Kb lines of Cr, Mn, Fe, Co, Ni, and Cu were implemented based on a multi-Lorentzian parameterisation and line energies are sampled from these Lorentzian components during the radiative transfer routine. In the optically thin regime, the SKIRT results match the intrinsic line profiles as measured in the laboratory. With a more complex 3D model that also includes kinematics, we find that the intrinsic line profiles are broadened and shifted to an extent that will be detectable with XRISM/Resolve; this model also demonstrates the importance of the intrinsic line shapes for constraining kinematics. We find that observed line profiles directly trace the cold-gas kinematics, without any additional radiative transfer effects. With the advent of the first XRISM/Resolve data, this update to the X-ray radiative transfer framework of SKIRT is timely and provides a unique tool for constraining the velocity structure of cold gas from X-ray microcalorimeter spectra., Comment: 6 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compact obscured nuclei (CONs) represent a significant phase of galactic nuclear growth. These galaxies hide growing SMBHs or unusual starbursts in their very opaque, extremely compact (r $<$ 100 pc) centres. They are found in approximately 30 % of the luminous and ultra-luminous infrared galaxy (LIRG and ULIRG) population. Here, we present high-resolution ALMA observations ($\sim$30 mas, $\sim$5 pc) of ground-state and vibrationally excited HCN towards ESO 320-G030 (IRAS 11506-3851). ESO 320-G030 is an isolated luminous infrared galaxy known to host a compact obscured nucleus and a kiloparsec-scale molecular wind. Our analysis of these high-resolution observations excludes the possibility of a starburst-driven wind, a mechanically or energy driven active galactic nucleus (AGN) wind, and exposes a molecular MDH wind. These results imply that the nuclear evolution of galaxies and the growth of SMBHs are similar to the growth of hot cores or protostars where gravitational collapse of the nuclear torus drives a MHD wind. These results mean galaxies are capable, in part, of regulating the evolution of their nuclei without feedback., Comment: 9 pages, 10 figures

Recent submillimeter observations have revealed signs of pc-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multi-wavelength observations in this paper series. The major axis position-velocity diagram (PVD) of CO(3-2) reproduces the observed faint emission at the systemic velocity, and our calculations confirm that this component originates from failed winds falling back to the disk plane. The minor-axis PVD of [CI]($^3P_1$-$^3P_0$), when created using only the gas with positive radial velocities, presents a sign of blue- and redshifted offset peaks similar to those in the observation, suggesting that the observed peaks indeed originate from the outflow, but that the model may lack outflows as strong as those in the Circinus galaxy. Similar to the observed HCN(3-2), the similar dense gas tracer HCO$^+$(3-2) can exhibit nuclear spectra with inverse P-Cygni profiles with $\sim$0.5 pc beams, but the line shape is azimuthally dependent. The corresponding continuum absorbers are inflowing clumps at 5-10 pc from the center. To detect significant absorption with a high probability, the inclination must be fairly edge-on ($\gtrsim$85$^\circ$), and the beam size must be small ($\lesssim$1 pc). These results suggest that HCN or HCO$^+$ and [CI] lines are effective for observing pc-scale inflows and outflows, respectively., Comment: 19 pages, 16 figures, 1 animation (see arXiv ancillary files), accepted for publication in ApJ

We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or [Ne V]-weak ($L_{\rm [Ne\,V]}$/$L_{\rm 12,AGN} \leq -$3.4) among 30 AGNs in local U/LIRGs. For the [O IV]-weak AGNs, we estimate their covering fractions of Compton-thick (CT; $N_{\rm H} \geq 10^{24}$ cm$^{-2}$) material with an X-ray clumpy torus model to be $f^{\rm (spec)}_{\rm CT} = 0.55\pm0.19$ on average. This value is consistent with the fraction of CT AGNs ($f^{\rm (stat)}_{\rm CT} = 53\pm12$%) among the [O IV]-weak AGNs in U/LIRGs and much larger than that in Swift/BAT AGNs ($23\pm6$%). The fraction of [O IV]-weak AGNs increases from $27^{+13}_{-10}$% (early) to $66^{+10}_{-12}$% (late mergers). Similar results are obtained with the [Ne V] line. The [O IV] or [Ne V]-weak AGNs in late mergers show larger $N_{\rm H}$ and Eddington ratios ($\lambda_{\rm Edd}$) than those of the Swift/BAT AGNs, and the largest $N_{\rm H}$ is $\gtrsim$10$^{25}$ cm$^{-2}$ at ${\log}\lambda_{\rm Edd} \sim -$1, close to the effective Eddington limit for CT material. These suggest that (1) the circumnuclear material in buried AGNs is regulated by the radiation force from high-$\lambda_{\rm Edd}$ AGNs on the CT obscurers, and (2) their dense material with large $f^{\rm (spec)}_{\rm CT}$ ($\sim 0.5 \pm 0.1$) in U/LIRGs is a likely cause of a unique structure of buried AGNs, whose amount of material may be maintained through merger-induced supply from their host galaxies., Comment: 13 pages, 10 figures, and 1 table, accepted for publication in ApJ

In this study, we examine photoionization outflows during the late stages of galaxy mergers, with a specific focus on the relation between observed velocity of outflowing gas and the apparent effects of dust extinction. We used the N-body/smoothed particle hydrodynamics (SPH) code ASURA for galaxy merger simulations. These simulations concentrated on identical galaxy mergers featuring supermassive black holes (SMBHs) of 10$^8$ M$_\odot$ and gas fractions of 30\% and 10 \%. From the simulation data, we derived velocity and velocity dispersion diagrams for the AGN-driven ionized outflowing gas. Our findings show that high-velocity outflows with velocity dispersions of 500 km s$^{-1}$ or greater can be observed in the late stages of galactic mergers. Particularly, in buried AGNs, both the luminosity-weighted outflow velocity and velocity dispersion increase owing to the apparent effects of dust extinction. Owing to these effects, the velocity--velocity dispersion diagrams display a noticeable blue-shifted tilt in models with higher gas fractions. Crucially, this tilt is not influenced by the AGN luminosity but emerges from the observational impacts of dust extinction. Our results imply that the observed high-velocity \OIII outflow exceeding 1000 km s$^{-1}$ in buried AGNs may be linked to the dust extinction that occurs during the late stages of gas-rich galaxy mergers., Comment: 14 pages, accepted for publication in ApJ

The Imaging X-ray Polarimetry Explorer (IXPE) detected X-ray polarization in the nearest Seyfert 2 galaxy, the Circinus galaxy, for the first time. To reproduce the IXPE results, we computed the degree of polarization based on two types of radiative hydrodynamic simulations: a parsec-scale three-dimensional model and a sub-parsec-scale axisymmetric model with a higher spatial resolution. In a series of papers, we confirmed that these models naturally explain the multi-wavelength observations of the Circinus galaxy from radio to X-rays. We used a Monte Carlo Simulation for Astrophysics and Cosmology code to compute the linear polarization of continuum emission. We found that the degree of polarization based on the parsec-scale radiation-driven fountain model was smaller than that observed with the IXPE. The degree of polarization based on the sub-parsec-scale model depends on the hydrogen number density of the disk ($d$), and the degree of polarization obtained from our simulation is consistent with that observed with the IXPE in the case of $\log d/\mathrm{cm}^{-3} \geq 13$. We investigate where the photons are Compton scattered and imply that the origin of the X-ray polarization in the Circinus galaxy is the outflow inside $0.01 \ \mathrm{pc}$. In this case, the degree of polarization may change over a timescale of approximately ten years., Comment: 11 pages, 8 figures, accepted for publication in ApJ

We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non-LTE line radiative transfer module in the extended SKIRT code using standard benchmarks. We find excellent agreement between the SKIRT results, the published benchmark results, and results obtained using the ray-tracing non-LTE line radiative transfer code MAGRITTE, which validates our implementation. We apply the extended SKIRT code on a 3D hydrodynamic simulation of a dusty AGN torus model and generate multi-wavelength images with CO rotational-line spectra against the underlying dust continuum. We find that the low-J CO emission traces the geometrically thick molecular torus, whereas the higher-J CO lines originate from the gas with high kinetic temperature located in the innermost regions of the torus. Comparing the calculations with and without dust radiative transfer, we find that higher-J CO lines are slightly attenuated by the surrounding cold dust when seen edge-on. This shows that atomic and molecular lines can experience attenuation, an effect that is particularly important for transitions at mid- and near-infrared wavelengths. Therefore, our self-consistent dust and non-LTE line radiative transfer calculations can help interpret the observational data from Herschel, ALMA, and JWST., Comment: 12 pages, 9 figures, Accepted for publication in A&A

Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecular gas is evident on sub-parsec scales. We calculate that less than 3% of this inflow is accreted by the black hole, with the rest being ejected by multiphase outflows, providing feedback to the host galaxy. The observations also reveal a dense gas disk surrounding the inflow; the disk is gravitationally unstable which drives the accretion into the central ~1 parsec., Comment: First release on Nov 3, 2023 in Science. 32 pages in one column = Main (13 pages, 4 figures) + Supplement (19 pages, 9 figures + 2 tables). This is the accepted version after peer review

We investigated dusty and dust-free gas dynamics for a radiation-driven sub-pc scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole $10^7 M_\odot$ and bolometric luminosity $10^{44}$ erg s$^{-1}$ based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven ``lotus-like'' multi-shell outflow is launched from the inner part ($r \lesssim 0.04$ pc) of the geometrically thin disk, and it repeatedly and steadily produces shocks as mass accretion continues through the disk to the center. The shape of the dust sublimation radius is not spherical and depends on the angle ($\theta$) from the disk plane, reflecting the non-spherical radiation field and nonuniform dust-free gas. Moreover, we found that the sublimation radius of $\theta \sim 20$-$60$ deg varies on a timescale of several years. The ``inflow-induced outflow" contributes the obscuration of the nucleus in the sub-parsec region. The column density of the dust-free gas is $N_{\rm H} \gtrsim 10^{22}$ cm$^{-2}$ for $r \lesssim 0.04$ pc. Gases near the disk plane ($\theta \lesssim 30$ degree) can be the origin of the Compton-thick component, which was suggested by the recent X-ray observations of AGNs. The dusty outflow from the sub-parsec region can be also a source of material for the radiation-driven fountain for a larger scale., Comment: 13 pages, 9 figures, accepted for publication in ApJ

Multi-wavelength observations of active galactic nuclei (AGNs) often reveal various time scales of variability. Among these phenomena, "changing-look AGNs" are extreme cases where broad emission lines become faint/bright or even disappear/emerge between multi-epoch observations, providing crucial information about AGN internal structures. We here focus on "changing-state" AGNs, specifically investigating the transition of optical spectra over years to tens of years. Based on the axisymmetric radiation-hydrodynamical simulations (Paper I) for the gas dynamics within the dust-sublimation radius, we investigate the spectral properties of ionized gas exposed to the radiation from an AGN with a 10^7 Msun supermassive black hole. We find significant time-dependent variations in the Balmer emission lines by utilizing post-process pseudo-three-dimensional calculations and the spectral synthesis code CLOUDY. The equivalent width of Halpha and Hbeta changes by a factor of 3, or the emission lines even disappear during 30 years for the same viewing angle. The time-dependent behaviour arises primarily from gas dynamics, particularly the formation of non-steady, radiation-driven outflows within the innermost region of the disc (r <10^-3 pc). The intricate interplay between non-spherical radiation sources at the core of AGNs and the dynamic behavior of gas within the dust sublimation radius gives rise to radiation-driven outflows. This non-steady outflow potentially contributes to the observed variability in Balmer line emissions over multi-year timescales in certain AGNs., Comment: 7 pages, 5 figures, accepted by MNRAS

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 ($\Delta J=\pm1$, $v=0-1$, $\lambda \sim 4.7$ $\mathrm{\mu 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 $\theta_\mathrm{obs} = 50-80$ $^{\circ}$. At the inclination angle $\theta_\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., Comment: 19 pages, 17 figures, 3 tables, accepted by The Astrophysical Journal

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., Comment: 17 pages, 12 figures, 2 tables, accepted for publication in ApJ

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 ($\xi$) at a low inclination spreads over 4 orders of magnitude in $\xi$, 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., Comment: 13 pages, 13 figures, accepted for publication in ApJ

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., Comment: single column, 36 pages, 16 figures, accepted for publication in the Astrophysical Journal

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 $\mu$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 $\mu$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., Comment: 18 pages, 19 figures, and 3 tables, accepted to appear on A&A, Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

The narrow-line region (NLR) consists of gas clouds ionized by the strong radiation from the active galactic nucleus (AGN), distributed in the spatial scale of AGN host galaxies. The strong emission lines from the NLR are useful to diagnose physical and chemical properties of the interstellar medium in AGN host galaxies. However, the origin of the NLR is unclear; the gas clouds in NLRs may be originally in the host and photoionized by the AGN radiation, or they may be transferred from the nucleus with AGN-driven outflows. For studying the origin of the NLR, we systematically investigate the gas density and velocity dispersion of NLR gas clouds using a large spectroscopic data set taken from the Sloan Digital Sky Survey. The [S II] emission-line flux ratio and [O III] velocity width of 9,571 type-2 Seyfert galaxies and 110,041 star-forming galaxies suggest that the gas density and velocity dispersion of NLR clouds in Seyfert galaxies (ne ~ 194 cm-3 and sigma([O III]) ~ 147 km s-1) are systematically larger than those of clouds in H II regions of star-forming galaxies (ne ~ 29 cm-3 and sigma([O III]) ~ 58 km s-1). Interestingly, the electron density and velocity dispersion of NLR gas clouds are larger for Seyfert galaxies with a higher [O III]/Hbeta flux ratio, i.e., with a more active AGN. We also investigate the spatially-resolved kinematics of ionized gas clouds using the Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey data for 90 Seyfert galaxies and 801 star-forming galaxies. We find that the velocity dispersion of NLR gas clouds in Seyfert galaxies is larger than that in star-forming galaxies at a fixed stellar mass, at both central and off-central regions. These results suggest that gas clouds in NLRs come from the nucleus, probably through AGN outflows., Comment: 15 pages, 13 figures, 1 table. Accepted for publication in PASJ

The nuclear obscurer of Active Galactic Nuclei (AGN) is poorly understood in terms of its origin, geometry and dynamics. We investigate whether physically motivated geometries emerging from hydro-radiative simulations can be differentiated with X-ray reflection spectroscopy. For two new geometries, the radiative fountain model of Wada (2012) and a warped disk, we release spectral models produced with the ray tracing code XARS. We contrast these models with spectra of three nearby AGN taken by NuSTAR and Swift/BAT. Along heavily obscured sight-lines, the models present different 4-20keV continuum spectra. These can be differentiated by current observations. Spectral fits of the Circinus Galaxy favor the warped disk model over the radiative fountain, and clumpy or smooth torus models. The necessary reflector (NH>10^25/cm^2) suggests a hidden population of heavily Compton-thick AGN amongst local galaxies. X-ray reflection spectroscopy is a promising pathway to understand the nuclear obscurer in AGN., Comment: Accepted in A&A. X-ray spectral models can be downloaded from https://github.com/JohannesBuchner/xars/blob/master/doc/warpeddisk.rst (warped disk) and https://github.com/JohannesBuchner/xars/blob/master/doc/wada.rst (radiative fountain model). The XARS code is at https://github.com/JohannesBuchner/xars/

We studied the absorption features of CO lines against the continuum originating from the heated dust in the obscuring tori around active galactic nuclei (AGNs). We investigated the formation of absorption lines corresponding to the CO rotational transitions using three-dimensional non-LTE line transfer simulations considering the dust thermal emission. As in Papers I--III of this series, we performed post-processed radiative transfer calculations using the "radiation-driven fountain model" (wada2016}, which yields a geometrically thick obscuring structure around the nucleus. This model is consistent with the spectral energy distribution of the nearest type-2 Seyfert galaxy, the Circinus galaxy. We found that the continuum-subtracted channel maps of $J = 4-3$ and higher transitions show absorption regions along the disk mid-plane for an edge-on viewing angle. The spectra consist of multiple absorption and emission features, reflecting the internal inhomogeneous and turbulent structure of the torus. The deepest absorption feature is caused by the gas on the near-side of the torus between $r =10$ and 15 pc, which is located in front of the AGN-heated dust inside $r \simeq 5$ pc. We also found that a spatial resolution of 0.5--1.0 pc is necessary to resolve the absorption features. Moreover, the inclination angle must be close to the edge-on angle (i.e., $\sim 85^\circ$) to observe the absorption features. The findings of the present study imply that combining our radiation-hydrodynamic model with high-resolution observations of CO (7-6) by ALMA can provide new information about the internal structure of the molecular tori in nearby AGNs., Comment: 21 pages, 12 figures, accepted by Astrophysical Journal

We present ALMA ~0.02"-resolution observations of the nucleus of the nearby (~14 Mpc) type-2 AGN NGC 1068 at HCN/HCO+/HNC J=3-2 lines, as well as at their 13C isotopologue and vibrationally excited lines, to scrutinize the morphological/dynamical/chemical/physical properties of dense molecular gas in the putative dusty molecular torus around a mass-accreting supermassive black hole. We confirm almost east-west-oriented dense molecular gas emission both morphologically and dynamically, which we regard as coming from the torus. Bright emission is compact (<3 pc), and low-surface-brightness emission extends out to 5-7 pc. These dense molecular gas properties are not symmetric between the eastern and western torus. The HCN J=3-2 emission is stronger than the HCO+ J=3-2 emission within the ~7 pc torus region, with an estimated dense molecular mass of (0.4-1.0)x10^6Msun. We interpret that HCN abundance is enhanced in the torus. We detect signatures of outflowing dense molecular gas and a vibrationally excited HCN J=3-2 line. Finally, we find that in the innermost (<1 pc) part of the torus, the dense molecular line rotation velocity, relative to the systemic velocity, is the opposite of that in the outer (>2 pc) part, in both the eastern and western torus. We prefer a scenario of counter-rotating dense molecular gas with innermost almost-Keplerian-rotation and outer slowly rotating (far below Keplerian) components. Our high-spatial-resolution dense molecular line data reveal that torus properties of NGC 1068 are much more complicated than the simple axi-symmetrically rotating torus picture in the classical AGN unification paradigm., Comment: 45 pages, 20 Figures. Accepted for publication in ApJ

We studied the magnetic field structures and dynamics of magnetized multiphase gas on parsec scales around supermassive black holes by using global 3D magnetohydrodynamics (MHD) simulations. We considered the effect of radiative cooling and X-ray heating due to active galactic nuclei (AGNs). The gas disk consists of a multiphase gas with (1) cold ($\leq 10^3$ K) and thin, and (2) warm ($\sim 10^4$ K) and thick components with a wide range of number densities. The turbulent magnetic energy at maximum is comparable to the thermal and turbulent kinetic energies in the turbulent motion. We confirmed that the turbulent velocity of the warm gas in the ambient cold gas is caused by magnetoconvective instability. The turbulent magnetic field due to magnetorotational instability (MRI) is developed in the disk, but the mean toroidal magnetic field dominates and supports in a quasi-steady state, where the plasma-$\beta$, the ratio between gas pressure and magnetic pressure, is low ($\beta < 1$). As often seen in adiabatic MHD simulations of rotating disks, the direction of the mean toroidal field periodically reverses with time even in multiphase gas structures. The direction reversal is caused by magnetic flux vertically escaping from the disk and by the combination of the MRI and the Parker instability., Comment: 21 pages, 1 table, 16 figures, accepted for publication in ApJ. Movies are available at https://astrophysics.jp/MHD_torus/

In Wada, Tsukamoto, and Kokubo (2019), we proposed for the first time that a new class of planets, "blanets" (i.e., black hole planets), can be formed around supermassive black holes (SMBHs) in the galactic center. Here, we investigate the dust coagulation processes and physical conditions of the blanet formation outside the snow line ($r_{snow} \sim $ several parsecs) in more detail, especially considering the effect of the radial advection of the dust aggregates. We found that a dimensionless parameter $\alpha = v_t^2/c_s^2$, where $v_t$ is the turbulent velocity and $c_s$ is the sound velocity, describing the turbulent viscosity should be smaller than 0.04 in the circumnuclear disk, to prevent the destruction of the aggregates due to collisions. The formation timescale of blanets $\tau_{GI}$ at $r_{snow}$ is, $\tau_{GI} \simeq$ 70-80 Myr for $\alpha = 0.01-0.04$ and $M_{BH} = 10^6 M_\odot$. The mass of the blanets ranges from $\sim 20 M_E$ to $3000 M_E$ in $r < 4$ pc for $\alpha = 0.02$ ($M_E$ is the Earth mass), which is in contrast with $ 4 M_E-6 M_E$ for the case without the radial advection. Our results suggest that \textit{blanets} could be formed around relatively low-luminosity AGNs ($L_{bol} \sim 10^{42}$ erg s$^{-1}$) during their lifetime ($\lesssim 10^8$ yr)., Comment: 17 pages, 7 figures, accepted by ApJ

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., Comment: 22 pages, 16 figures, accepted by ApJ

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., Comment: 17 pages, 11 figures, Accepted for publication in ApJ

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 "snow line", located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the "radial drift barrier" 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'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's lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs., Comment: 15 pages, 4 figures, Accepted for publication in ApJ

It is widely believed that in active galactic nuclei (AGNs) a supermassive black hole with an accretion disk is surrounded by an optically and geometrically thick torus at sub-parsec scale. However, it is not clear how is the mass supply toward the central engine caused and how it is related with the internal structures of the tori. The magnetic field in the tori may contribute to the accretion process via the magneto-rotational instability (MRI). Using global three dimensional magnetohydrodynamic (MHD) simulations taking the effects of X-ray heating and radiative cooling into account studied the numerical resolution for azimuthal direction for MRI driving. We found that a strongly magnetized disk consisted of a cold ($< 10^3$ K) and warm ($10^4$ K) gas is developed in about 30 rotational periods. We also found in high resolution model that the mean azimuthal magnetic fields reverse their direction quasi-periodically. We confirmed that the typical wave length of the MRI should be resolved with a least 20 azimuthal grid cells., Comment: 9 pages, 6 figures, Accepted for publication in Galaxies as part of the Special Issue The Power of Faraday Tomography

We used the Atacama Large Millimeter/Submillimeter Array (ALMA) to map the CO(3-2) and [CI](1-0) lines, as well as their underlying continuum emission, from the central $\sim 200$ pc region of the Circinus galaxy that hosts the nearest type 2 Seyfert-class active galactic nucleus (AGN), with a spatial resolution of $\sim 6-15$ pc. The lines and continuum-emitting regions consist of a circumnuclear disk (CND; 74 pc $\times$ 34 pc) and spiral arms. The distribution of the continuum emission revealed a temperature-dependent dust geometry and possibly polar dust elongation in the torus region. The molecular mass of the CND is $M_{\rm H2} \sim 3 \times 10^6~M_\odot$ with a beam-averaged H$_2$ column density of $\sim 5 \times 10^{23}$ cm$^{-2}$ toward the AGN position, which contributes significantly to the nuclear obscuration. The [CI](1-0)/CO(3-2) ratio at the AGN position is unusually high, suggesting an X-ray dominated region-type chemistry. We decomposed the observed velocity fields into rotational and dispersion components, and revealed multi-phase dynamic nature in the $r \lesssim 10$ pc torus region, i.e., the diffuse atomic gas is more spatially extended along the vertical direction of the disk than the dense molecular gas. Through comparisons with our model predictions based on the radiation-driven fountain scheme, we indicate that atomic outflows are the driver of the geometrical thickness of the atomic disk. This supports the validity of the radiation-driven fountain scheme in the vicinity of this AGN, which would explain the long-lasting mystery, the physical origin of the AGN torus., Comment: 26 pages, 23 figures (incl. 2 appendix figures), 7 tables; Accepted for publication in ApJ

We investigate the properties of the ionized gas irradiated by an active galactic nucleus (AGN) based on our "radiation-driven fountain" model for the nearest type-2 Seyfert galaxy, the Circinus galaxy (Wada et al. 2016). We conducted "quasi-three dimensional" spectral analysis using the spectral synthesis code CLOUDY and obtained the surface brightness distributions of lines, such as H$\alpha$, H$\beta$, [OIII], [NII], and [SII] for the central 16-parsec region. The ionized regions observed based on these lines show a conical morphology around the rotation axis, even if we do not phenomenologically postulate the presence of an optically thick "torus". This region also shows non-uniform internal structures, reflecting the inhomogeneous structure of fountain flows. Using ionization diagnostic diagrams, we investigated the spectral properties of the ionized gas. The diagrams based on the line ratios of [NII]/H$\alpha$ and [SII]/H$\alpha$ show that most regions of the cone have the same properties as those in the narrow line regions (NLRs) in AGNs, whereas using [OI]/H$\alpha$, the central 10-pc regions are rather LINER-like. The gas density, temperature, and ionizing parameter in regions identified as "NLR" are typically $n \sim 300-1500$ cm$^{-3}$, $T \sim 1-3 \times 10^4 $ K, and $ U \sim 0.01$, respectively. The morphology and [OIII] intensity are similar to the base of the observed [OIII] cone in the Circinus galaxy, implying some physical connections with the origin of the $\sim100$ parsec scale NLR., Comment: 10 pages, 7 figures, accepted by ApJ

We present the results of our ALMA Cycle 4 high-spatial-resolution (0.04-0.07") observations, at HCN J=3-2 and HCO+ J=3-2 lines, of the nucleus of NGC 1068, the nearby prototypical type 2 active galactic nucleus (AGN). Our previous ALMA observations identified the compact emission of these lines at the putative location of the torus around a mass-accreting supermassive black hole. We now report that we have detected the rotation of this compact emission, with the eastern and western sides being redshifted and blueshifted, respectively. Unlike the previously reported CO J=6-5 emission, both the morphological and dynamical alignments of the HCN J=3-2 and HCO+ J=3-2 emission are roughly aligned along the east-west direction (i.e., the expected torus direction), suggesting that these molecular lines are better probes of a rotating dense molecular gas component in the torus. The western part of the torus exhibits larger velocity dispersion and stronger emission in the HCN J=3-2 and HCO+ J=3-2 lines than the eastern part, revealing a highly inhomogeneous molecular torus. The dense molecular gas in the torus and that of the host galaxy at 0.5-2.0" from the AGN along the torus direction are found to be counter-rotating, suggesting an external process happened in the past at the NGC 1068 nucleus., Comment: 7 pages, 3 figures, accepted for publication in ApJL. Figures 1 and 2 are degraded. High resolution version is available at http://www2.nao.ac.jp/~masaimanishi/N1068/N1068bMemuH.pdf

In this study, we investigate the line emissions from cold molecular gas based on our previous "radiation-driven fountain model" (Wada 2016), which reliably explains the spectral energy distribution of the nearest type 2 Seyfert galaxy, the Circinus galaxy. Using a snapshot of the best-fit radiation-hydrodynamic model for the central r < 16 pc, in which non-equilibrium X-ray-dominated region chemistry is solved, we conduct post-processed, non-local thermodynamic equilibrium radiation transfer simulations for the CO lines. We obtain a spectral line energy distribution with a peak around J=6, and its distribution suggests that the lines are not thermalized. However, for a given line-of-sight, the optical depth distribution is highly non-uniform between $\tau \ll 1$ and $\tau \gg 1$. The CO-to-H2 conversion factor (X_CO), which can be directly obtained from the results, is not a constant and depends strongly on the integrated intensity, and it differs from the fiducial value for local objects. X_CO exhibits a large dispersion of more than one order of magnitude, reflecting the non-uniform internal structure of a "torus." We also found that the physical conditions differ between grid cells on a scale of a few parsecs along the observed lines of sight; therefore, a specific observed line ratio does not necessarily represent a single physical state of the ISM., Comment: 13 pages, 12 figures, accepted by ApJ

On both observational and theoretical grounds, the disk of our Galaxy should be accreting cool gas with temperature ~<10^5 K via the halo at a rate ~1 M_sun/yr. At least some of this accretion is mediated by high velocity clouds (HVCs), observed to be traveling in the halo with velocities of a few 100 km/s and occasionally impacting the disk at such velocities, especially in the outer regions of the Galaxy. We address the possibility of particle acceleration in shocks triggered by such HVC accretion events, and the detectability of consequent non-thermal emission in the radio to gamma-ray bands and high-energy neutrinos. For plausible shock velocities ~300 km/s and magnetic field strengths ~0.3-10 muG, electrons and protons may be accelerated up to ~1-10 TeV and ~30-10^3 TeV, respectively, in sufficiently strong adiabatic shocks during their lifetime of ~10^6 yr. The resultant pion decay and inverse Compton gamma-rays may be the origin of some unidentified Galactic GeV-TeV sources, particularly the "dark" source HESS J1503-582 that is spatially coincident with anomalous HI structure known as "forbidden-velocity wings". Correlation of their locations with star-forming regions may be weak, absent or even opposite. Non-thermal radio and X-ray emission from primary and/or secondary electrons may be detectable with deeper observations. The contribution of HVC accretion to Galactic cosmic rays is subdominant, but could be non-negligible in the outer Galaxy. As the thermal emission induced by HVC accretion is likely difficult to detect, observations of such phenomena may offer a unique perspective on probing gas accretion onto the Milky Way and other galaxies., Comment: ApJ, in press after minor revision; 16 pages, 6 figures

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 "radiation-driven fountain" (Wada 2012) with supernova feedback. A double hollow cone structure is naturally formed without postulating a thick "torus" 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$\mu 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 $\mu 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., Comment: 9 pages, 5 figures. Accepted for ApJL. A movie file for Fig.5b can be downloaded from http://astrophysics.jp/Circinus/

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., Comment: 14 pages, 11 figures, accepted for publication in PASJ

Active galactic nuclei (AGNs) are believed to be obscured by an optical thick "torus" that covers a large fraction of solid angles for the nuclei. However, the physical origin of the tori and the differences in the tori among AGNs are not clear. In a previous paper based on three-dimensional radiation-hydorodynamic calculations, we proposed a physics-based mechanism for the obscuration, called "radiation-driven fountains," in which the circulation of the gas driven by central radiation naturally forms a thick disk that partially obscures the nuclear emission. Here, we expand this mechanism and conduct a series of simulations to explore how obscuration depends on the properties of AGNs. We found that the obscuring fraction f_obs for a given column density toward the AGNs changes depending on both the AGN luminosity and the black hole mass. In particular, f_obs for N_H \geq 10^22 cm^-2 increases from ~0.2 to ~0.6 as a function of the X-ray luminosity L_X in the 10^{42-44} ergs/s range, but f_obs becomes small (~0.4) above a luminosity (~10^{45} ergs/s). The behaviors of f_obs can be understood by a simple analytic model and provide insight into the redshift evolution of the obscuration. The simulations also show that for a given L_AGN, f_obs is always smaller (~0.2-0.3) for a larger column density (N_H \geq 10^23 cm^-2). We also found cases that more than 70% of the solid angles can be covered by the fountain flows., Comment: 17 pages, 12 figures, appeared in ApJ. corrected typos in Fig. 7 and Fig. 8

The disk wind, which is powered by the radiation force due to spectral lines (line force), is studied for broad absorption line (BAL) quasars. We investigate the structure of the disk wind based on the non-hydrodynamic method and compare with wind properties inferred from X-ray observations of BAL quasars. In this paper, we apply the stellar wind theory to the initial condition (the mass outflow rate at the base of the wind). We found the funnel-shaped winds with a half opening angle of 50^{circ} for the case of epsilon=0.3-0.9 and M_{BH}=10^{7-8.5}M_odot, where epsilon is the Eddington ratio and M_{BH} is the black hole mass. Thus, the absorption features are observed for an observer of which a viewing angle is around 50^{circ}. A probability of BAL quasars is 7-11%, which is roughly consistent the abundance ratio of BAL quasars, 10-15%. Here, the probability is estimated by the solid angle, that the absorbing features would be detected, divided by 4pi. In contrast, if the Eddington ratio is smaller than 0.01 or if the black hole is very massive, M_{BH} < 10^9M_{odot}, the disk wind is not launched due to the less effective line force. Then, the quasars are identified as non-BAL quasars independently of the observer's viewing angle., Comment: 10 pages, 6 figures. Accepted for publication in PASJ

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., Comment: 16 pages (in emulateapj format), incl. 13 figures and 1 table; Accepted for publication in The Astrophysical Journal

We propose a plausible mechanism to explain the formation of the so-called "obscuring tori" 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 "fountain", 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., Comment: 12 pages, 9 figures, accepted publication in ApJ

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., Comment: 33 pages, 11 figures, to be published in New Astronomy

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., Comment: 12 pages, 9 figures, accepted for publication in ApJ

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'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., Comment: Accepted for publication in the Astrophysical Journal Letters, vol. 734, 2011

We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic simulations in a `live' 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., Comment: 13 pages, 12 figures, accepted by ApJ. Higher resolution of ms.pdf is available at http://d.pr/Nvjk A targzipped Supplementary movies is available at http://d.pr/TV61

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/\alpha 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") galaxies., Comment: 35 pages, 13 figures, accepted for publication in ApJ

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., Comment: 4 pages, 3 figures, Proceedings of IAU Symposium 270, Computational Star Formation

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 '3-kpc arm', '135-km s$^{-1}$ arm' and 'Connecting arm' but also clumpy features such as `Bania clumps', 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., Comment: accepted to PASJ

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)., Comment: 12 pages, 10 pages, accepted by ApJ

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 \theta_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., Comment: 15 pages, 12 figures, accepted by Astrophysical Journal