Active galaxies contain a supermassive black hole at their center that grows by accreting matter from the surrounding galaxy. The accretion process in about the central 10 parsecs has not been directly resolved in previous observations because of the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using submillimeter interferometry. A dense inflow of molecular gas was evident on subparsec scales. We calculated 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. Our observations also reveal a dense gas disk surrounding the inflow that is gravitationally unstable, which drives the accretion into about the central 1 parsec.

We investigated dusty and dust-free gas dynamics for a radiation-driven sub-parsec-scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole 10⁷ M _⊙ and bolometric luminosity 10⁴⁴ erg s⁻¹ based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven “lotus-like” multi-shell outflow is launched from the inner part (r ≲ 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 (θ) from the disk plane, reflecting the nonspherical radiation field and nonuniform dust-free gas. Moreover, we found that the sublimation radius of θ ∼ 20°–60° varies on a timescale of several years. The “inflow-induced outflow” contributes to the obscuration of the nucleus in the sub-parsec region. The column density of the dust-free gas is N _H ≳ 10²² cm⁻² for r ≲ 0.04 pc. Gases near the disk plane (θ ≲ 30°) 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.

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 (≤10³ K) and thin, and (2) warm (∼10⁴ 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-β, the ratio between gas pressure and magnetic pressure, is low (β < 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.

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” (Wada et al. 2016), 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 midplane 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 ≃ 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., ≳85°) 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.

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 multiwavelength 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 / cm − 3 ≥ 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 pc. In this case, the degree of polarization may change over a timescale of approximately 10 yr.

A recent hydrodynamic model, the 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 have confirmed that this picture is consistent with multiwavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovibrational absorption lines (ΔJ = ± 1, v = 0 − 1, λ ∼ 4.7 μm) could probe the 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 θ _obs = 50°–80°. At the inclination angle θ _obs = 77°, we observe multi-velocity components: inflow (v _LOS = 30 km s⁻¹), systemic (v _LOS = 0 km s⁻¹), and outflows (v _LOS = −75, − 95, and −105 km s⁻¹). The inflow and outflow components (v _LOS = 30 and −95 km s⁻¹) are collisionally excited at the excitation temperatures 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 few parsec region of AGN tori, and the observations can probe the gas circulation inside the tori.