Explore the vast range of titles available.

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 the torus-scale gas dynamics around AGNs and the DOGs’ lifetime remains unclear. We obtained the evolution of the spectral energy distributions (SEDs) of a galaxy merger system with AGN feedback from postprocessed 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 108 M ⊙. We found that the infrared luminosity of the system reaches ultra- and hyperluminous infrared galaxy classes (1012 and 1013 L ⊙, respectively). The DOG phase corresponds to a state in which the AGNs are buried in dense gas and dust, with the infrared luminosity exceeding 3.3 × 1012 L ⊙. We also identified subcategories of DOGs, namely bump and power-law DOGs, from the SEDs and their evolution. The bump DOGs tend to evolve to power-law DOGs over 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 are distributed nonspherically 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 Myr, suggesting that the observed DOG phase is a brief aspect of galaxy mergers.

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 107 M ⊙ and bolometric luminosity 1044 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 ≲ 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 ≳ 1022 cm−2 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 investigated the influence of the Eddington ratio on sub-parsec-scale outflows in active galactic nuclei (AGNs) with supermassive black holes (SMBHs) masses of 107 M ⊙ using two-dimensional radiation hydrodynamics simulations. When the range of the Eddington ratio is γ 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 are a time-dependent phenomenon, 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 reveal insights into the Eddington ratio dependence and the tendency toward 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.

We performed N-body/smoothed particle hydrodynamics simulations of isolated spiral galaxies with various bulge-to-disk mass ratios (Mbulge/Mdisk) from 0.02 to 0.2 to investigate mass transport from galactic scales (10 kpc) down to circumnuclear disk scales (≲100 pc). Our analysis revealed these main findings: (1) Gravitational torque from stellar spiral arms causes gas accretion with ∼1 M⊙ yr−1 along the gas spiral arms from a few kiloparsecs to a few 100 pc scale. The density of accreting gas is a few 100 cm−3, comparable to the gas arms. The pressure gradient force is over an order of magnitude weaker than the stellar gravitational torque. (2) Gravitational torque from barred structure causes episodic gas clump accretion with ∼1 M⊙ yr−1 on timescales of 10 Myr from the kiloparsec to a few 100 pc scale. The densities of these clumps exceed 700 cm−3, and this accretion occurs along elliptical orbits with a delayed phase relative to the bar potential. (3) Episodic gas clumpy accretion is important for galactic center instability, confirmed by Mbulge/Mdisk = 0.02 but not by Mbulge/Mdisk = 0.1 and 0.2. This difference occurs because in the bulge-dominated potentials, bar instability is suppressed, and rapid gas clump accretion does not occur. These findings suggest that gas clump accretion events driven by bars could be a source of high-density gas to the galactic center of the spiral galaxy, potentially promoting temporary activity in the galactic center.

We evaluated the jet power and the density of ambient matter in 3C 84 by using the momentum balance along the jet axis and the transonic condition for the cocoons observed at two different scales (approximately 1 and 6 pc). For the inner cocoon, we precisely determined the ratio of jet power to ambient density Lj/na to be (0.3–0.7) × 1043 erg s−1 cm3. Similarly, for the outer cocoon, we found that this value is more than an order of magnitude larger at (0.9–3.7) × 1044 erg s−1 cm3. This indicates that the outer cocoon is formed by a powerful jet that propagates through an ambient density of 20–300 cm−3 with a jet power of 1045−46.5 erg s−1. On the other hand, the inner cocoon is formed by a weaker jet with a power of 1043−44 erg s−1, propagating through a relatively low-density environment of 6–20 cm−3. These results suggest that (1) with respect to the difference in na, it appears to support the hypothesis that the inner cocoon, recently formed about 10 yr ago, is expanding in the low-density environment created by the jet emitted about 25–50 yr ago; (2) to achieve the short-lived and high Lj that generated the outer cocoon, a large mass-accretion rate is required over a short period to activate the jet. These may imply an extreme accretion event driven by the tidal disruption events of massive stars and/or the disk instability.

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 ∼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 the magnetic bubbles for our model that corresponds to a necessary condition for outflow growth. We found that the north–south asymmetrical 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 nonsteady 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, e.g., NGC 1377.

Although the X-ray spectra of Seyfert 1 galaxies exhibit absorption lines of He-like iron and H-like iron at blueshifted velocities of approximately 500 km s−1, the physical origin of these absorption lines remains uncertain. In this study, we performed X-ray radiative transfer based on the subparsec-scale thermally driven outflows. The initial step involved calculating the photoionization equilibrium using the Cloudy code, which is based on three-dimensional radiative hydrodynamic simulations. Subsequently, X-ray radiative transfer was performed using the Monte Carlo Simulation for Astrophysics and Cosmology code. Our findings indicate that when the angle of inclination ranges within 55°–65°, the transmitted component of the X-ray spectrum displays absorption lines of He-like and H-like iron, exhibiting a blueshift of approximately 500 km s−1. The results suggest that the absorption lines are generated by a photoionized gas within 0.005 pc. Additionally, the results indicate that the scattered component of the X-ray spectrum exhibits emission lines originating from neutral iron fluorescence, He-like iron, and H-like iron. The emission lines are broadened by approximately 7000 km s−1 due to the Keplerian rotation. Furthermore, the model reproduced the H-like iron and H-like iron absorption lines in NGC 3783 observed by the Chandra High Energy Transmission Grating.

In this study, we examine photoionization outflows during the late stages of galaxy mergers, with a specific focus on the relation between the observed velocity of outflowing gas and the apparent effects of dust extinction. We used the N-body/smoothed particle hydrodynamics code ASURA for galaxy merger simulations. These simulations concentrated on identical galaxy mergers featuring supermassive black holes of 108 M ⊙ and gas fractions of 30% and 10%. From the simulation data, we derived velocity and velocity dispersion diagrams for the active galactic nuclei (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, 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 [O iii] λ5007 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.

Recent submillimeter observations have revealed signs of parsec-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 multiwavelength 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 [C i](3 P 1–3 P 0), when created using only the gas with positive radial velocities, presents a sign of blueshifted 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 ∼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 (≳85°), and the beam size must be small (≲1 pc). These results suggest that HCN or HCO+ and [C i] lines are effective for observing parsec-scale inflows and outflows, respectively.

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 logd/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.